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PROPAGATION OF SINGULARITIES IN MANY-BODY 
SCATTERING IN THE PRESENCE OF BOUND STATES 

ANDRAS VASY 



Abstract. In this paper we describe the propagation of singularities of tem- 

- , pered distributional solutions u £ S' of {H — X)u = 0, where if is a many-body 

^J ■ Hamiltonian H = A + V, A>0,V = J]^ Va, and A is not a threshold of H, 

under the assumption that the inter-particle (e.g. two-body) interactions Va 

are real- valued polyhomogonoous symbols of order —1 (e.g. Coulomb-type with 

the singularity at the origin removed). Here the term 'singularity' refers to a 

r\ , • microlocal description of the lack of decay at infinity. Thus, we prove that the 

^^ ' set of singularities of m is a union of maximally extended broken bicharacteris- 

^~N tics of H. These are curves in the characteristic variety oi H — \, which can be 

(-H , quite complicated due to the existence of bound states. We use this result to 

describe the wave front relation of the S-matrices. We also analyze Lagrangian 
properties of this relation, which shows that the relation is not 'too large' in 
terms of its dimension. 



> 

^^ ■ 1. Introduction 

m ■ 

^D ' In this paper we describe the propagation of singularities of generahzed eigen- 

^— ^ . functions of a many-body Hamiltonian H = A + V, V = ^^ Va, on K" under the 

assumption that the inter-particle interactions Va are real-valued polyhomogeneous 
^s . symbols of order —1 (e.g. Coulomb-type with the singularity at the origin removed). 

■^ I More precisely, we use the 'many-body scattering wave front set' WFsc(m) at infin- 

"^ , ity for tempered distributions u € iS'(M"), and prove that for u € 5'(K.") satisfying 

2 ' {H—X)u = 0, WFsc(w) is a union of maximally extended generalized broken bichar- 

acteristics of H, broken at the collision planes. Here WFscCu) provides a microlocal 
description of the lack of decay of u modulo 5(R"), similarly to how the usual wave 
front set describes distributions modulo C°° functions. 

The definition of generalized broken bicharacteristics is quite technical due to 
the presence of bound states in the subsystems. However, if these bound states are 
absent, our definition becomes completely analogous to Lebeau's definition [^ for 
the wave equation in domains with corners. Indeed, in this case the propagation 
result itself, which was proved in |3J], is a direct (C°°-type) analogue of Lebeau's 
result for the propagation of analytic singularities for solutions of the wave equation 
in domains with corners. 

If there are bound states in the subsystems, but either the set of thresholds 
is discrete, or _ff is a four-body Hamiltonian, the geometry of generalized broken 
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2 ANDRAS VASY 

bicharacteristics is not much more complicated than in Lebeau's setting. The gen- 
eral definition reflects that when particles collide, the total energy as well as the 
external momentum is preserved. The complication in the presence of bound states 
is that kinetic energy is not preserved, even asymptotically. In summary, our re- 
sults provide a connection between quantum and classical objects, just as Lebeau's 
results connect the wave equation and geometric optics. Note that these results 
only provide the answer as to where in phase space the singularities of generalized 
eigenfunctions may be located; they leave open the question of what these singu- 
larities are like, i.e. we do not have FlO-type results. Such results exist at least in 
certain 3-body settings p3| , M, see the remarks in the next section. In addition, 
we show that in our setting, if the set of thresholds is discrete, then the broken 
bicharacteristic relation has the correct Lagrangian geometry to give rise to such 
FIO results. 

We also prove the corresponding result in the 'limiting absorption principle' 
setting, namely that under certain assumptions on WFsc(/), i?(A±iO)/ are defined, 
and WFsc(-R(A + iO)f) is a subset of the image of WFsc(/) U i?_(A) under the 
forward broken bicharacteristic relation. Here i?_(A) is the outgoing 'radial set'. 
Such a result makes the 'radial- variable' propagation estimates that have been used 
in many-body scattering, especially as derived in the works of Gerard, Isozaki and 
Skibsted [^ 0, more precise. 

We use this result to analyze that the wave front relation of the scattering ma- 
trices (S- matrices). These connect the incoming and outgoing data of generalized 
eigenfunctions of H, so one expects that their singularities are described by consid- 
ering the limit points of generalized broken bicharacteristics 7 = 7(t) as t —> ±00. 
In fact, in addition to the propagation of singularities result, the only ingredient 
that is required for this analysis is a good approximation for the incoming Poisson 
operators with incoming state a near the incoming region, and similar results for 
the outgoing Poisson operators with, say, outgoing state /?. In general, one expects 
a WKB-type construction, essentially as in Hadamard's parametrix construction. 
Indeed, this is what Melrose and Zworski do in the geometric two-body type set- 
ting, [^ . In the Euclidean many-body setting this construction has been done by 
Skibsted |32] in the short-range and by Bommier m in the long-range setting, in 
the latter case by adopting the construction of Isozaki and Kitada fl^, at least 
imder the assumption that the energies of the states a, (3 are below the continuous 
spectrum of the corresponding subsystem Hamiltonians. Such a construction is 
unnecessary if Vc are Schwartz, for then the product decomposition is sufhciently 
accurate to give a good approximation for the Poisson operator. We thus obtain 
the following result. 

Theorem. Suppose that H is a many-body Hamiltonian, and A is not a threshold 
or (LP'-) eigenvalue of H. Suppose also that either a and (3 are channels such that 
the corresponding eigenvalues £„, e^, of the subsystem Hamiltonians H"^ , H^ , are 
in the discrete spectrum of H"^ and H^ respectively, or that Vc is Schwartz for all 
c. Then the wave front relation of the S-matrix S'^q(A), is given by the generalized 
broken bicharacteristic relation of H as stated precisely below in Theorem \2.^ . 

Special cases, which have already been analyzed, include the free-to- free S-matrix 
in three-body scattering p3. pi p^, or indeed in many-body scattering under the 
additional assumption that there are no bound states in any subsystem Q| . In these 
cases the wave front relation is given by the broken geodesic relation, broken at the 
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collision planes, on §"^^ at distance tt. In both cases, one can naturally extend 
the results to geometric many-body type problems on asymptotically Euclidean 
manifolds. 

Also, Bommier ||l[ and Skibsted ||3^ have shown that the kernels of the 2-cluster 
to free cluster and 2-cluster to 2-cluster S-matrices are smooth (except for the diag- 
onal singularity if the 2-clusters are the same), and previously Isozaki had showed 
this in the three-body setting Q. We remark that (under our polyhomogeneous 
assumption) the proofs of Bommier and Skibsted in fact show that the 2-cluster 
to same 2-cluster S-matrices are (non-classical, if the potentials are long-range) 
pseudo-differential operators, which differ from the identity operator by operators 
of order —1 if the potentials are short range. In our geometric normalization this 
means that they are Fourier integral operators associated to the geodesic flow on 
the sphere at infinity to distance tt (along the cluster). 

It may seem that the (generalized) broken bicharacteristic relation is rather large 
since a single bicharacteristic can be continued in many ways after it hits a collision 
plane. In fact, this is not the case, since many bicharacteristics do not hit these 
collision planes at all. A more precise statement would be the that the wave front 
relation of the S-matrices is given by a union of Lagrangian submanifolds, so at 
least it has the same structure as if the S-matrices were (sums of) Fourier integral 
operators. Indeed, we prove that if the set of thresholds of H is discrete, e.g. if 
there are no bound states in any subsystems, then the wave front relations of the 
S-matrices are given by finite unions of smooth Lagrangian relations. 

We remark that the results of this paper would remain valid if we assumed only 
that Va € S~P{X'^), p > 0, as customary. In fact, the proof of the propagation 
of singularities for generalized eigenfunctions remains essentially unchanged, and 
the only difference in the above Theorem is that the parametrix for the Poisson 
operators is not as explicit, cf. Q; instead, one needs to use the constructions 
of Isozaki-Kitada ||l^ (as presented by Skibsted and Bommier) directly. The rea- 
son for the polyhomogeneous assumption is that the proofs are somewhat nicer, 
especially in notation, and it is a particularly natural assumption to make in the 
compactification approach we adopt. 

Our main tool in proving the propagation of singularities results consists of 
microlocally positive commutator estimates, i.e. on the construction of operators 
which have a positive commutator with H in the part of phase space, say U, where 
we wish to conclude that a generalized eigenfunction u has no scattering wave front 
set. These commutators are usually negative in another region of phase space, 
namely backwards (or forwards, depending on the construction) along generalized 
broken bicharacteristics through U. We thus assume the absence of this region 
from WFgc(u), and conclude that the positive commutator region, U, is also miss- 
ing from WFgc(u). Such techniques have been used by Hormander, Melrose and 
Sjostrand |12| , p4| to show the propagation of singularities for hyperbolic equations 
(real principal type propagation) such as the wave equation, including in regions 
with smooth boundaries. Indeed, the best way to interpret our results is to say 
that if — A is hyperbolic at infinity. In two-body scattering the analogy with the 
wave equation in domains without boundary is rather complete; this was the ba- 
sis of Melrose's proof of propagation estimates for scattering theory for 'scattering 
metrics' in p3. In many-body scattering, the lack of commutativity of the appro- 
priate pseudo-differential algebra, even to top order, makes the estimates (and their 
proofs) more delicate. We remark that, as can be seen directly from the approach 
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we take, the wave front set estimates can be easily turned into microlocal estimates 
on the resolvent considered as an operator between weighted Sobolev spaces; wave 
front set statements are a particularly convenient way of describing propagation. 
Indeed, there is some freedom in the precise definition of the wave front set; see 



the remarks preceeding the statement of Theorem 2.2. The alternative definitions 
differ slightly, but agree for generalized eigenfunctions, and the corresponding prop- 
agation results are based on the same positive commutator estimates. Thus, the 
reader may find the explicit microlocal estimates of Section M particularly clear. 
However, piecing together these estimates to describe propagation is more cumber- 
some than making simple geometric statements based on wave front sets; for this 
reason we emphasize the latter. 

Positive commutator estimates have also played a major role in many-body scat- 
tering starting with the work of Mourre p^ , Perry, Sigal and Simon p^ , Froese 
and Herbst [g|, Jensen |Q, Gerard, Isozaki and Skibsted [g[ 0| and Wang [g8|. In 
particular, the Mourre estimate is one of them; it estimates i[H, w ■ D^ -{- D^ ■ w\. 
This and some other global positive commutator results have been used to prove the 
global results mentioned in the first paragraph about some of the S-matrices with 
initial state in a two-cluster. They also give the basis for the existence, unique- 
ness and equivalence statements in our definition of the S-matrix by asymptotic 
expansions; these statements are discussed in [^ |6[ Q in more detail. 

More delicate (and often time-dependent) commutator estimates have been used 
in the proof of asymptotic completeness. This completeness property of many-body 
Hamiltonians was proved by Sigal and Soffer, Graf, Derezihski and Yafaev under 



different assumptions on the potentials and by different techniques ||2^, |29|, 31, ^ 
|[ g, ^. In particular, Yafaev's paper |3|] shows quite explicitly the importance of 
the special structure of the Euclidean Hamiltonian. This structure enables him to 
obtain a (time-independent) positive commutator estimate, which would not follow 
from the indicial operator arguments of [^ Q and the present paper, and which 
is then used to prove asymptotic completeness. 

We briefly outline of the positive commutator proofs. They consist of two parts: 
first, the construction of a symbol (and an associated pseudo-differential operator) 
that we claim has a microlocally positive commutator with iJ, when localized in 
energy, and second, the proof that the commutator is indeed positive in the ap- 
propriate part of the phase space. The first part has much in common with the 
analysis of generalized broken bicharacteristics, since both are intimately connected 
to various Hamilton vector fields associated to the subsystems. The second part is 
essentially a microlocal version of the proof of the Mourre estimate that was ob- 
tained by Froese and Herbst ||^. Indeed, we could follow the full indicial operator 
version of this proof, as was done in |gj and |3J]. However, in the presence of 
bound states in the proper subsystems we would have to rely more heavily on the 
approximate product structure of the Hamiltonian at each cluster. Indeed, in the 
full geometric problem considered in |Q , the L^ eigenvalues and eigenfunctions of 
the indicial operator of H can vary as one moves along the collision planes, which 
makes even the description of the bicharacteristics more complicated. However, it 
turns out that in the Euclidean many-body setting, after an explicit calculation 
of the indicial operators, we can use the Mourre estimate explicitly in the normal 
(non-commutative) variables at each cluster, and use the standard Poisson bracket 
formula for the commutator in the tangential (commutative) variables. Since it 
eliminates the need to present arguments that are essentially simple (microlocal) 
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modifications of the Froese-Herbst proof of the Mourre estimate, we adopt the 
second approach. 

The structure of the paper is the foUowing. In the next section we set up the 
framework of many-body scattering, mostly foUowing the book of Derezihski and 
Gerard, g], and we state the precise results on the propagation of singularities. 
In Section p| we relate this to the Melrose's approach to scattering via compact- 
ification EJ. We also include in this section the outline of the resolved space 
construction and the definition of the algebra of many-body scattering differential 
operators from p4| . As we indicate in this section, from an algebraic point of 
view, it is the lack of commutativity of this algebra to 'top weight' (at infinity) 
that gives rise to the breaks in the generalized broken bicharacteristics along which 
singularities propagate. In Section we recall the definition and basic properties 
of many-body scattering pseudo-differential operators as well as the definition of 
the corresponding wave front set from ||34| , and we analyze the characteristic va- 
riety of many -body Hamiltonians. Microlocal elliptic regularity is stated here in 
Corollary ^^. In Section || we describe generalized broken bicharacteristics, and in 
Section |6| we explain the positive commutator argument that is the key to our prop- 
agation results. Propagation of singularities itself, stated in Theorem |2.2| , is proved 
in Section 0. Sections p|-[lO| turn this result i nto t heorems on the resolvent and the 
scattering matrices, stated in Theorems 2.4- 2.10 . Finally, in Appendix ^ we show 



that if the set of thresholds of H is discrete, or if i/ is a four-body Hamiltonian, 
generalized broken bicharacteristics are piecewise integral curves of the Hamilton 
vector fields, with only a finite number of breaks, and then in Appendix H we 
analyze the Lagrangian structure of the broken bicharacteristic relation. 

Most results of this paper were announced in |3^ ; this paper contains the detailed 
proofs. 

I am very grateful to Andrew Hassell, Richard Melrose and Maciej Zworski for 
numerous very fruitful discussions; in particular, I would like to thank Richard 
Melrose for his comments on this paper. I am grateful to Maciej Zworski for in- 
troducing me to the work of Gilles Lebeau |]l9J. If there are no bound states in 
any subsystems, many-body scattering is 'philosophically' and, to a certain extent, 
technically (e.g. the structure of generalized broken bicharacteristics) is very sim- 
ilar to the wave equation in domain with corners. Thus, Lebeau's paper played 
an important direct role in my paper p3 , and remained philosophically important 
while working on the present manuscript. My joint projects with Andrew Hassell, 
as well as our discussions in general, provided very valuable insights into the bro- 
ken bicharacteristic geometry, especially towards understanding their Lagrangian 
structure. I would also like to thank Rafe Mazzeo, Erik Skibsted and Jared Wunsch 
for helpful discussions, their encouragement and for their interest in this research. 

2. Notation and detailed statement of results 

Before we can state the precise definitions, we need to introduce some basic (and 
mostly standard) notation. We refer to g for a very detailed discussion of the setup 
and the basic results. We consider the Euclidean space M", and let g be the standard 
Euclidean metric on it. We assume also that we are given a (finite) family X of 
linear subspaces Xa, a G I, oi R" which is closed under intersections and includes 
the subspace Xi = {0} consisting of the origin, and the whole space Xq = K". Let 
X"" be the orthocomplement of Xa- We write ga and g" for the induced metrics on 
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Xa and X° respectively. We let tt" be the orthogonal projection to X°- , iTa to Xa- 
A many-body Hamiltonian is an operator of the form 

(2.1) if = A + ^(7r'')*K; 

ael 

here A is the positive Laplacian, Vo — 0, and the Va are real-valued functions in an 
appropriate class which we take here to be poly homogeneous symbols of order —1 
on the vector space Xa to simplify the problem: 

(2.2) Va e 5phg(X'^). 

In particular, smooth potentials Va which behave at infinity like the Coulomb po- 
tential are allowed. Since {'K°')*Va is bounded and self-adjoint and A is self-adjoint 
with domain _ff^(R") on L^ — L^(M"), H is also a self-adjoint operator on L^ with 
domain H^{W). We let R{\) = {H - \)-^ for A G C \ R be the resolvent of H. 

There is a natural partial ordering on / induced by the ordering of X"" by inclu- 
sion. (Though the ordering based on inclusion of the Xa would be sometimes more 
natural, here we use the conventional ordering, we simply write Xa C Xi, if the op- 
posite ordering is required.) Let /i = {1} (recall that Xi = {0}); 1 is the maximal 
element of /. A maximal element of /\ /i is called a 2-cluster; I2 denotes the set of 
2-clusters. In general, once Ik has been defined for fc = 1, . . . , m — 1, we let /m (the 
set of m-clusters) be the set of maximal elements of /^ = / \ Uj!!7//fc, if /„ is not 
empty. If I'^^ ~ {0} (so /^-i-i is empty), we call H an m-body Hamiltonian. For 
example, if / ^ {0, 1}, and for all a,b ^ {0, 1} with a ^ b we have XaC] Xb = {0}, 
then 7f is a 3-body Hamiltonian. The A^-cluster of an iV-body Hamiltonian is also 
called the free cluster, since it corresponds to the particles which are asymptotically 
free. 

Our goal is to study generalized eigenfunctions oi H, i.e. solutions u E 5'(M") of 
{H — X)u — 0. Since H ~ X is an elliptic partial differential operator with smooth 
coefficients, {H - X)u e C°°{R") implies that u € C°°(R"). Thus, the place where 
such u can have interesting behavior is at infinity. Analysis at infinity can be viewed 
either as analysis of uniform properties, or as that of properties in the appropriate 
compactification of K". We adopt the second point of view by compactifying R" 
as in Q. Thus, we let 

(2.3) X = Xo = S'l 

to be the radial compactification of M" (also called the geodesic compactification) 
to a closed hemisphere, i.e. a ball, and §"^^ = 9§". Recall from ||2^ that RC : 
R" -^ S" is given by 

(2.4) RC{w) = {l/{l + \w\''y/^,w/{l + \w\Y^^)e§lcW'+\ weM". 

Here we use the notation RC instead of SP, used in |E3[, to avoid confusion with 
the standard stereographic projection giving a one-point compactification of R". 
We write the coordinates on R^ — Xa® X°- as (wa, w°). We let 

(2.5) Xa = C1(RC(X„)), Ca=Xan dS^. 

Hence, Ca is a sphere of dimension Ua — \ where ria = dim Xq. We also let 

(2.6) C = {Ca: ae I}. 

Thus, Co = 9§" = §"^^, and a < 6 if and only if Cb C Ca- Since throughout this 
paper we work in the Euclidean setting, where the notation X, Xa, etc., has been 
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used for the (non-compact) vector spaces, we always use a bar, as in X, Xa, etc., to 
denote the corresponding compact spaces. In particular, even when talking about 
general compact manifolds with boundary in the following sections, recalling the 
results of |34|, we will write them as X. 

We note that if a is a 2-cluster then Ca n Cb = unless Ca C Cb, i.e. b < a. We 
also define the 'singular part' of Ca as the set 

(2.7) Ca,sing = Ub^a{Cb H Ca), 

and its 'regular part' as the set 

(2.8) C^ = Ca \ Llb^aCb = Ca\ Causing- 

For example, if a is a 2-cluster then Ca,sing — and C^ — Ca- 

We usually identify (the interior of) §" with M". A particularly useful boundary 
defining function of S" is given by a; g C°°(S" ) defined as a; = r^^ = |w|^^ (for 
r > 1, say, smoothed out near the origin); so S"^^ = 9§" is given by a; = 0, x > 
elsewhere, and dx ^ a.t §""1. We write S';^g(§+) and 5';^,g(M") interchangeably 
We also remark that 

(2.9) ^^hg(S+) = 2:-"C°° (§';). 

We recall that under RC, C°°(§"), the space of smooth functions on S!f. vanishing 
to infinite order at the boundary, corresponds to the space of Schwartz functions 
5(M"), and its dual, C-°°{E:l), to tempered distributions 5'(R"). We also have the 
following correspondence of weighted Sobolev spaces 

(2.10) H^/iK) = H^'^ = H^'\^'') = {w)-^H^{W) 

where {w) = (1 + |wp)^/^. 

Corresponding to each cluster a we introduce the cluster Hamiltonian iJ" as an 
operator on L^(X°) given by 

(2.11) iJ" = A + ^H, 

b<a 

A being the Laplacian of the induced metric on X"' . Thus, if iJ is a A''-body 
Hamiltonian and a is a fc-cluster, then H"^ is a [N -I- 1 — /c)-body Hamiltonian. 
The L^ eigenfunctions of if" (also called bound states) play an important role 
in many-body scattering; we remark that by a result of Froese and Herbst, B, 
speCpp(7f ) C (—00,0] (there are no positive eigenvalues). Moreover, speCpp(i/°) 
is bounded below since iJ° differs from A by a bounded operator. Note that 
X^ = {0}, H^ = 0, so the unique eigenvalue of H^ is 0. 

The eigenvalues of if can be used to define the set of thresholds of H^ . Namely, 
we let 

(2.12) A, = u<aspecpp(ii'') 
be the set of thresholds of H'^ , and we also let 

(2.13) a;, = A, U spccpp(ii'') = U6<a specpp(i7''). 

Thus, e Aq for a 7^ and Ka C (—00, 0]. It follows from the Mourre theory (see 
e.g. H, ^) that Aa is closed, countable, and speCpp(ii°) can only accumulate at 
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Aa- Moreover, -R(A), considered as an operator on weighted Sobolev spaces, has a 
hniit 

(2.14) i?(A ± ^0) : H^^{Sl) ^ H^+^^'' (81) 

for / > 1/2, /' < —1/2, from either half of the complex plane away from 

(2.15) A = AiUspeCpj,(i7). 

In addition, L^ eigenfunctions of iJ" with eigenvalues which are not thresholds are 
necessarily Schwartz functions on X" (in fact, they decay exponentially, see Q). 
We also label the eigenvalues of H"^, counted with multiplicities, by integers m, and 
we call the pairs a — {a,m) channels. We denote the eigenvalue of the channel a 
by Eq, write ipa for a corresponding normalized eigenfunction, and let Cq be the 
orthogonal projection to V'a in L^{X'^). 

The phase space in scattering theory is the cotangent bundle T*R". Again, it is 
convenient to consider its appropriate partial compactification, i.e. to consider it as 
a vector bundle over S" . Thus, consider the set of all one- forms on R" of the form 

n 

(2.16) ^ajdwj 

where aj e C°°(§") (we drop RC from the notation as usual). This is then the set 
of all smooth sections of a trivial vector bundle over §" , with basis dwi, . . . , dwn- 
Following Melrose's geometric approach to scattering theory, see [^, we consider 
this as the (dual) structure bundle, and call it the scattering cotangent bundle of §" , 
denoted by ^''T*^^. Note that T*W can be identified with M" x R" via the metric 
g; correspondingly ^'^T*E>^ is identified with §" x R", i.e. we simply compactified 
the base of the standard cotangent bundle. We remark that the construction of 
scy*gn jg completely natural and geometric, just like the following ones, see p3[ , 
or Section S for a summary. 

However, in many-body scattering ^'^T*S^ is not the natural place for microfo- 
cal analysis for the very same reason that introduces the compressed cotangent 
bundle in the study of the wave equation on bounded domains. We can see what 
causes trouble from both the dynamical and the quantum point of view. Regarding 
dynamics, the issue is that only the external part of the momentum is preserved 
in a collision, the internal part is not; while from the quantum point of view the 
problem is that there is only partial commutativity in the algebra of the associated 
pseudo-differential operators, even to top order. To rectify this, we replace the full 
bundle "'^T^, Sl=C'^x R" over C'^ C S"-^ by "'^T^, X^ = C^ x X„, i.e. we consider 

(2.17) ^''t*Sl = yJa'^'T^.^Xa. 

Over Cq, there is a natural projection tTq : ^^T^iEi^ -^ ^^T^,Xa corresponding to 
the pull-back of one-forms; in the trivialization given by the metric it is induced 
by the orthogonal projection to Xa in the fibers. By putting the tTq together, we 
obtain a projection tt : ^'^Tg„_i§" — > ^'^T*'§,\. We put the topology induced by tt on 
sc2i*gn q;^}^ig definition is analogous to that of the compressed cotangent bundle 



in the works of Melrose, Sjostrand |24| and Lebeau |19 on the wave equation in 
domains with smooth boundaries or corners, respectively. 

We also recall from |g3[ that the characteristic variety So (A) of A— A is simply the 
subset of ^^T^„-iWl where g — \ vanishes; g being the metric function. If Ai — {0}, 
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the compressed characteristic set oi H — X will be simply tt{Y.q(X)) C ^^T*§" . In 
general, all the bound states contribute to the characteristic variety. Thus, we let 

(2.18) Sb(A) - {Ch e ''T^.Xb : A - |6p G speCppij"} C ^'T^.Xb; 

note that |^bp is the kinetic energy of a particle in a bound state of H^. If Ca C 
Ch, there is also a natural projection TTba : ^^Tq, X}, -^ ^"^T^iXa (in the metric 
trivialization we can use the orthogonal projection Xi, -^ Xa as above), and then 
we define the characteristic set of iJ — A to be 

(2.19) ±{X) = UaSa(A), ±aiX} = UcoC„ TT^a (Sb(A)) H ^'T^^Xa, 

SO t{X) C ''''t*Sl. We let ttj, be the restriction of tti, : """T^.^Xb -> t{X) to S6(A). 

We next recall from |Q the definition of generalized broken bicharacteristics in 
case there are no bound states in any of the subsystems. In fact, in this case the 
word 'generalized' can be dropped; for the generalized broken bicharacteristics have 
a simple geometry as stated below. First, note that the rescaled Hamilton vector 
field of the metric function g, i.e. 

(2.20) 2{w)i, ■ dy, e V(T*]R") = V(]R" x K") 

extends to a smooth vector field, """^Hg e V(''''T*Sl) = V{Sl x E"), with S'^ consid- 
ered as the radial compactification of M"; in fact, ^"^Hg is tangent to the boundary 

Definition. Suppose Ai = {0}, and / = [a, f3] is an interval. We say that a 
continuous map 7 : / ^ ^(-^) is a broken bicharacteristic of i/ — A if there exists 
a finite set of points tj£l,a = to<ti<...< tk-i < tfc = /3 such that for each 
Jj l\{t-.t-+i) is the image of an integral curve of ^'^Hg in So(A) under n. If / is an 
interval (possibly M), we say that 7 : / — + S(A) is a broken bicharacteristic oi H — X, 
if the restriction of 7 to every compact subinterval of / is a broken bicharacteristic 
in the above sense. 

Here 7(1) C S(A) — 7r(So(A)) corresponds to the conservation of kinetic energy 
in collisions (since there are no bound states), and the use of the compressed space 
scy*§n shows that external momentum is conserved in the collisions. It turns out, 
see [Q, that 7 is essentially the lift of a broken geodesic on S" of length n (if 
/ = K, otherwise shorter), broken at the collision planes, i.e. at C. In particular, 
even if / = R, it has only a finite number of breaks, and in fact, there is a uniform 
bound on the number of such breaks (depending only on the geometry, i.e. on C, 
not on 7). 

The definitions are less explicit if Ai ^ {0}, but they essentially still state that 
the total energy and the external momentum are preserved in collisions. Thus, 
generalized broken bicharacteristics will be continuous maps 7 defined on intervals 
J, 7 : / — !■ S(A) with certain appropriate generalization of the integral curve condi- 
tion described above. In order to take the bound states into consideration, we also 
need to consider the rescaled Hamilton vector fields ^^H^ of the metric gb in the 
subsystem b. Thus, under the inclusion map 

(2.21) lb : ^^ra^Xfc -> ''T^,^Sl 

induced by the inclusion Xb ^^ M" in the fibers, {ib)*^'^Hg = ^'^Hg (i.e. the re- 
striction of the vector field ^^Hg to ^^T^ Xb, considered as a subset of ^^T^ S" ). 
Thus, we require that lower bounds on the Hamilton vector fields ^^H^ applied to 
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TT-invariant functions, i.e. to functions / G C°°(^'^Tg„_iS") such that /(^) — /(^') if 
7r(^) — 7r(^'), imply lower bounds on the derivatives of /^r along 7. Here /tt is the 
function induced by / on ^'^T*§" , so / = /^r o tt. 

Definition 2.1. A generalized broken bicharacteristic of iJ — A is a continuous 
map 7 : / — > S(A), where / C M is an interval, such that for all to E I and for 
each sign + and — the following holds. Let ^0 = 7(^0), suppose that ^0 G ^'^T^'Xa- 
Then for aU 7r-invariant functions / S C°°('"'T*„_,E,l), 

(2.22) D^iU o 7)(io) > infri/gVdo) : lo e 7r-i(^o), ^ C ^j. 

Here D± are the one-sided lower derivatives: if g is defined on an interval /, 
iD±g){to) = liminit^to±ig{t) - g{to))/{t - to). 

Although it is not apparent, this definition is equivalent to the previous one if 
Ai = {0}. Moreover, in four-body scattering, even if Ai ^ {0}, one can describe 
the generalized broken bicharacteristics piecewise as projections of integral curves 
of ^'^Hg. In general many-body scattering, the lack of conservation of kinetic energy 
makes such a description harder, but if Ai is discrete, we obtain a description that 



parallels the one above. More precisely, in Theorem 5.10 we prove the following. 
Suppose that Ai is discrete and 7 : R — )■ S(A) is a continuous curve. Then 7 is a 
generalized broken bicharacteristic of iJ — A if and only if there exist io < ii < ^2 < 
... <tk such that 7l[tj,tj+i], as weU as 7l(-oo,to] and 7l[tfc,+oo), are the projections 
of integral curves of the Hamilton vector field ^''Hg for some a. In addition, there 
is a uniform bound on k (independent of 7) , depending only on C and Ai . Similar 
results hold if the interval of definition, M, is replaced by any interval. 

As mentioned in the introduction, 'singularities' (i.e. lack of decay at infinity) of 
u € S' are described by the many-body scattering wave front set, WFsc(u), which 
was introduced in ||34| , and which describes u modulo Schwartz functions, similarly 
to how the usual wave front set describes distributions modulo smooth functions. 
Just as for the image of the bicharacteristics, ^'^T*E>'^ provides the natural setting 
in which WFsc is defined: WFsc(u) is a closed subset of scy*§n ^ xhe definition of 
WFsc(u) relies on the algebra of many-body scattering pseudo-differential opera- 
tors, also introduced in Q. There are several possible definitions of WFsc, all of 
which agree for generalized eigenfunctions of H, but the one given in pj| that is 
modelled on the fibred-cusp wave front set of Mazzeo and Melrose pO] enjoys many 
properties of the usual wave front set. A slightly different definition, for which the 



crucial property in the last line of ( 4.35 ) still holds, was discussed in [^; this is the 



only essential prope rty fo r the positive commutator estimates in this paper. Also, 



the discussion after ( 9.13 ) indicates why 'finite order' versions of the wave front set, 
i.e. versions in which we only require a fixed (though possibly high) number of the 
symbol estimates for ps.d.o.'s, somewhat akin to the constructions of |6|, would be 
helpful; nonetheless, the Mazzeo-Melrose definition appears to be the most natural 
one from the point of view of general microlocal analysis. 

We recall the precise definitions in Section Q here we also translate our results 
into statements on the S-matrices where the usual wave front set can be used. 
We remark that in the two-body setting, when ^^t*^^ = ^'^T*^^, WFgc is just 
the scattering wave front set WFgc introduced by Melrose, |2^, which in turn is 
closely related to the usual wave front set via the Fourier transform. Thus, for 
(10,0 e "'=Tg*„_i§!J:, considered as S"-i x R" = dS^ x R", (w,^ ^ WFsc(w) means 
that there exists G C°°{Sl) such that 0(w) ^ and T{^u) is C°° near ^. If we 
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employed the usual conic terminology instead of the compactified one, we would 
think of (^ as a conic cut-off function in the direction uj. Thus, WFgc at infinity is 
analogous to WF with the role of position and momentum reversed. The definition 
of WFsc(u) is more comphcated, but li u — ij;{H)v for some V" S C'^{^) (any 
other operator in ^^^' {^.,0) would do instead of '>p(H)), then the following is 
a sufficient condition for (w,^o) & ^'^T^i Xa, considered as C'^ x Xa, not to be in 
WFsc(u). Suppose that there exists </>°e C°°(§!^), 0(u;) 7^ 0, and p e C^(X„), 
p(Ca) ^ 0, and ii7:-)*p)T{(bu) e 5(R") = S{Xo). Then (w,Ca) ^ WFse(u). We 
also remark that we state all of the following results for the absolute wave front 
sets (i.e. we work modulo Schwartz functions), but they have complete analogues 
for the relative wave front sets (working modulo weighted Sobolev spaces); indeed, 
it is the latter that is used to prove the results on the former. 

Our main result is then the following theorem, in which we allow arbitrary thresh- 
olds, and which describes the relationship between WFsc(w) and generalized broken 
bicharacteristics, if, for example, {H — X)u — 0. Note that if {H — X)u = 0, then 
u = tp{H)u for tp G C^(M), xj} = 1 near A, so the above description of WFsc('«) is 
applicable. 

Theorem 2.2. Let u E 5'(R"), A (^ Ai. Then 

(2.23) WFsc{u)\WFsc{{H~X)u) 

is a union of maximally extended generalized broken bicharacteristics of H — X in 
±{X)\WFsc{{H-X)u). 

We remark that the statement of the theorem is empty at points ^0 G ^^T^, Xa 
at which '"^i?g(Co) = for some ^0 G '^b^i^o) and some b with Ca C Cb- Indeed, at 
such points the constant curve (7(i) — S,o for all t in some interval) is a generalized 
broken bicharacteristic. A simple calculation shows that the set of these points ^0 
is R+{X) Ui?_(A), where 

(2.24) 

R±W = {e e ""'T^'^Xa : 3b, Ca (lCb,X- T(e)2 e speCppiH"), ±t(C) > 0} 

are the incoming (+) and outgoing (— ) radial sets respectively, and r is the sc-dual 
variable of the boundary defining function x, so in terms of the Euclidean variables 

(2.25) r = -ff; 

\w\ 

see the next section for further details. Hence, the theorem permits singularities 
to emerge 'out of nowhere' at the radial sets. Although we do not prove that 
this indeed does happen, based on general principles, this appears fairly likely. 
Moreover, the optimality of Theorem 2^ if Ai = {0} follows from |^, Q, see the 



remarks about this in ||3J|; the amplitude of the reflected 'wave' is given (to top 
order) by the appropriate subsystem S-matrix. 

There is a similar result for WFsc(w), u — R{X + iO)/; namely that WFgc(u) \ 
WFsc(/) is the image of WFsc(/) U i?_(A) under forward propagation, if e.g. / € 
H"^'^, s > 1/2. The set R-{X) appears here since there can be maximally extended 
generalized broken bicharacteristics which are either not disjoint from i?_(A), or 
simply whose closure is not disjoint from _R_ (A). In particular, even if / is Schwartz, 
WFsc(u) is not necessarily a subset of i?_(A), rather a subset of its image under 
forward propagation. Indeed, by duality, this is exactly what gives rise to the 
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conditions on WFsc(/) under which u ~ R{X + iOjf can be defined. To make it 
easier to state these results, we make the following definition. 

Definition 2.3. Suppose K C S(A). The image $+(/\) of K under the forward 
broken bicharacteristic relation is defined as 

(2.26) 

^+{K) ={^0 G S(A) : 3 a generalized broken bicharacteristic 7 : (—00, to] -^ ^W 



s.t. 7(to) = eo, 7((-oo, to]) nK^d}}. 

The image $_(_ftr) of K under the backward broken bicharacteristic relation is 
defined similarly, with [to, +00) in place of {—00, to]. 

Note that ^+{K) = U^^k^+{{£,}) directly from the definition. The result on 
the boundary values of the resolvent is then: 

Tiieorem 2.4. Suppose that A ^ A, / e 5(R"), and let u = R{X + M)f. Then 
WF 5c(u) C <f>+(i?_(A)). Moreover, R{\ + iO) extends by continuity to v e 5'(K") 
with WFsc{v) n $_(i?+(A)) = 0, and for such v, 

(2.27) WFsc(i?(A + ^0)1-) c $+(WFsc(«)) U $+(i?_(A)). 

The scattering matrices S'/3q(A) of H with incoming channel a, outgoing channel 
(3 can be defined either via the wave operators, or via the asymptotic behavior 
of generalized eigenf unctions. It was shown in pq | that the two are the same, 
up to normalization (free motion is factored out in the wave operator definition); 
here we briefly recall the second definition. We first state it for short-range Vc {Vc 
polyhomogeneous of order —2 for all c). Thus, for A £ (e^, cx)) \ A and g G C^{C'^), 
there is a unique u S 5'(R") such that {H — X)u = 0, and u has the form 

(2.28) u = e-*^'^^^V-'i™'^»/2((7r'')>„)i;_ + R{X + iO)f, 

where v^ £ C°°(§!Ji), V-\c, = g, and / e i/sT'^^^^''(§+), e' > 0. The Poisson 
operator Pq_+(A) is the map 

(2.29) Pa.+ (A) : C^{C'^) ^ 5'(M") defined by Pa,+ (A)5 = u. 

The term R[X + iO)/ has distributional asymptotics of a similar form 'at the chan- 
nel (i\ i.e. of the form eV^-^'^''~r-'i"°^''/2((7r'')*V'/3)w/3,+, see |3§ for the precise 
definitions. 

Only minor modifications are necessary for Vc G ST^^iA-^'^)- Namely, write 

(2.30) Ia=Yl ^b, la = (r-a/a)|c„ € C°°{Ca), 

la is C°° near C'^ with simple vanishing at C^ (since b ^ a means Cg <^ Cb, hence 
Co n Cb C Ca, sing), so To/a is C°° there. Then the asymptotics in ( ^.2g ) must be 
replaced by 



(2.31) e-*V^^V-'^™'^"/V*^~"/2xA=7r((^a-)*^^)^_ 

plus lower order terms; see Section |^ for details. 

The scattering matrix S^aiX) maps g — v- {c^ to W/3,+ |c' • It is also given by the 
formula 

(2.32) ^0„(A) = i ((g - A)P^,„(A))*P„,+ (A), 

2i^X- 6/3 
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A > niax(eQ,e^), A ^ A. Here Pp-{X) is a microlocalized version of the outgoing 
Poisson operator, microlocalized near the outgoing region for /3, i.e. where r is near 
— i^A — e/3, see Q. In fact, we can simply take P^__(A) to be a microfocal (cut- 
off) parametrix for P^__(A). This formula is closely related to that of Isozaki and 
Kitada [0- 

A very good parametrix, Pq^-|_(A), for Pq,+ (A) in the region of phase space where 
T is close to \/X — ta has been constructed by Skibsted |Q in the short-range and by 
Bommier |^ in the long-range setting, under the assumption that Cq € spec^(iJ'^). 
If we instead assume that the Vc are all Schwartz, then the trivial (product type) 
construction gives the desired parametrix. Their constructions enable us to de- 
duce the structure of the S-matrices immediately from our propagation theorem, 



Theorem 2.4 via (2.32) and 



(2.33) P„,+ (A) = P„^+(A) - P(A + z0)(i7 - A)P„,+ (A). 

Since the parametrix (near the incoming or outgoing sets) is important for turning 
the results on the propagation of singularities to wave front set results, in all our 
results on the Poisson operators and S-matrices 5/3^ (A) in this paper we make the 
following assumption: 

(2.34) either e„ e speCrf(iJ") and e^ e specjH''), or K G SiX") for aU c. 

It is easy to describe the wave front set of Pa,+ {X)g, g G C^°°{C'^), near its 
'beginning point', i.e. near the {a, -|-)-incoming set. Namely, it is the union of 
integral curves of ^'^H^ (in ^'^T^,Xa) (which are in particular bicharacteristics of 
H — X, hence broken bicharacteristics), one integral curve for each ( S WF{g) C 
S*C'g^\ we denote these integral curves by 7a,-(C)- It is actually convenient to 
replace the parameter t of the integral curve by s, the arclength parameter of its 
projection to Ca- The relationship between these two is that if we write s — S{t), 
then S solves the ODE dS/dt = 2(A - e^ - T(7(i))2)i/2, xhe reparameterized 
integral curves are then given by 

(2.35) 
Ta = \JX- eaCos(s- sq), {va, fJ-a) = v^A - £„ sin(s - So) exp((s - SQy''Hg){C) 

where s G {sq, sq -I- tt). This defines ja.-iC) up to replacing t hy t — ti for any fixed 
ti £ M, so we are abusing the notation slightly. Due to (2.33D, Theorem 2.4 describes 



WFsc(Pa,+ (A)5) elsewhere. A similar result also applies for P/3,_(A); in this case 
one simply has to replace the range of the arclength parameter by s e (sq — tt, sq). 
We denote the corresponding integral curves by 7/3^+ (C). 

Definition 2.5. The forward broken bicharacteristic relation with initial channel 
a is defined to be the relation TZa,+ C S(A) x S*C'^ given by 

(2.36) 7^„,+ -{(C,C)eS(A)x5*C^ $-({0) n 7a,-(C) y^ 0}- 

The backward broken bicharacteristic relation with initial channel /3, denoted by 
Tiff,- is defined similarly, with $+ in place of <&_ and 7/3.+ in place of 7q._. Finally, 
the forward broken bicharacteristic relation with initial channel a, final channel (3, 
T^pa C S*Cl X S*C'^ is defined as the composite relation of TZa,+ and 'R-^^_'- 

(2.37) 

nfj,,^{iCX')^S*C',xS*C',: 3CeS(A),(C,e)e7^„.+ , (C^0e7^^.-}. 
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Note that (C, ^ T^a.+ thus means that there exists a generaUzed broken bichar- 
acteristic 7 : R —> S(A) and io G M such that 7|(_oo,to] is given by 7q,_(C), and 
£, & 7{[to, +00)). Thus, Tla.+ should be thought of as the relation induced by $_ 
'at channel a' as time goes to —00. 

una Bx Aisa relation, K C A, hy n{K) we mean {£, ^ B : 3C G K, (^, () e 
7^}. Similarly, if C/ C B, by 7^-l(C/) we mean {C e ^ : 3^ e L/, (C, C) € 7^}. We 
can n{K) the image of K under 7^. Thus, if if C S'*C;, 

7^„.+(A0 - le e S(A) : 3C e if, $-({a) n 7a,-(C) y^ 0}, 

andif C/C t{\), then 

7^a,V(c/) - {C e 5*c; : ae e [/, $-({a) n 70,- (0 ^ 0}- 

This definition, ( 2.33| ) and Theorem 2.4 immediately prove the following propo- 
sition. 



Proposition 2.6. Suppose that H is a many-body Hamiltonian, A ^ A, and (2.34) 
holds. Suppose also that g £ C^{C'^). Then 

(2.38) WF5c(Pa.+ (A).g) \ i?+(A) C $+(i?_(A)). 

In addition, P^ _|_(A) extends by continuity from C'^{C'^) to distributions g G 
C-°°(C^) with (R+iX) X WF(5)) n 7^„,+ = (i.e. 7^„,+ (WF(5)) n i?+(A) = 0;. i/ 
5 is such a distribution, then 

(2.39) WF5c(Pa.+ (A)g) \ i?+(A) c $+(i?_(A)) U 7^„,+ (WF(.9)). 

One of the main features of ( ^.3^ ) is that in general WFsciPa,+ iX)g) cannot be 
expected to be contained in the radial sets; one also has to include the image of 
the outgoing radial set under forward propagation in the statement. As a corollary. 



(2.32) shows that in general 6*^0 (A) does not map smooth incoming data to smooth 
outgoing data. However, if /3 is a two-cluster channel, every generalized broken 
bicharacteristic 7 such that for some tg G M, 7|(-oo,to] i^ given by 7^,+ (C), C ^ S*Cb: 
is actually equal to 7/3_+(^) for all times, and, b being a 2-cluster, 7^ _|_(C) never 
intersects the radial sets, and as i — > ±00, -ffj^+{(){t) goes to i?;p(A). Thus, if /3 is a 
2-cluster, a is any cluster, Sf3a{X) maps smooth functions to smooth functions. On 
the other hand, if /3 is the free channel 0, then the absence of positive thresholds 
gives a similar conclusion. 



Corollary 2.7. Suppose that H is a many-body Hamiltonian, A ^ A, and (2.34) 



holds. Suppose g e C^{C'^). Then WF(5'0„(A)g) C 7^^_^_(i?_(A)). Thus, if P is 
the free channel or it is a two-cluster channel, then Sf3a{X)g is C°°. 

Our theorem on the wave front relation of the S-matrix is then the following. 



Theorem 2.8. Suppose that H is a many-body Hamiltonian, A ^ A, and ( 2.34 ) 
holds. Then S'/5q(A) extends by continuity from C^ {C"^) to distributions g G C~°°(C'^) 
with 7?.Q.-|-(WF(g)) n i?-|-(A) = 0. If g is such a distribution, then 

(2.40) WF(5^,(A)g) C UJ^liR^iX)) U 7^;5„(WF(g)). 

If Ai = {0}, then maximally extended generalized broken bicharacteristics are 
essentially the lift of generalized broken geodesies on §"~^ of length tt, so in this 
case we recover the following result of Q . 



MANY-BODY SCATTERING 15 

Corollary 2.9. // no subsystem of H has bound states and A > 0, then the wave 
front relation of Sqo{X) is given by the broken geodesic relation on E>"^^, broken at 
C, at distance n. 

It may seem that our results are too weak in the sense that some broken bicharac- 
teristics can be continued in many ways when they hit a collision plane. However, it 
should be noted that not every broken bicharacteristic hits a collision plane; indeed, 
if the dimension of a collision plane increases, more broken bicharacteristics will hit 
it, but each will generate a lower dimensional family of continuations. We can make 
this precise in terms of a Lagrangian characterization of the wave front relation of 
the S-matrices. To do so, recall first that the wave front relation of operators 
mapping C^iC'^) to C-°°{C'f,) is often understood as a subset of r*C^ x T*C'^ \ 
corresponding to the wave front set of its Schwartz kernel, with the sign of the 
second component, the one in T*C'g^, switched. For example, in this sense, the wave 
front relation of an FIO, say P, is given by a homogeneous canonical relation, i.e. 
by a Lagrangian submanifold A of r*C^ x T*C'^, Lagrangian with respect to the 
usual twisted symplectic form ujh — uja, that is conic with respect to the diagonal 
R+ action in T*C'ij x T*C'^. Here cua denotes the canonical symplectic form on 
T*Ca, etc. The wave front set mapping properties of such an operator are that Pu 
is defined if (0 x WF(m)) n A = 0, and if Pu is thus defined, 

WF(Pu) cA(WF(m)) U {C e T*C^ \ : (^, 0) G A} 

= {e e T*C', \ : (e, 0) e A or 3^' E WF(u), (C, O e A}. 

Note that WF(u) is a conic subset of T*C'^, so it is better to regard it as a subset 
of 5'*C^ (quotient out by the M"*" action on T*C'^ \ 0); indeed, this is what we have 
done. But every element ( of 5'*C^ (thought of as the quotient bundle) has several 
representatives ^ in T*C'^. Correspondingly, the wave front relation, understood 
as a relation connecting S*C'^ and S*C'^, i.e. as a subset of S*C[^ x S*C'^, only 
determines the corresponding conic relation on {T*Cl \ 0) x {T*C'^ \ 0) up to the 
rescaling of one factor with respect to the other. Thus, our results by themselves 
cannot pinpoint the wave front set of the Schwartz kernel of Spa (A) itself, only its 
image under the quotient map with respect to the M+ x R+ action. Note, however, 
that A n (T*C'i, X 0) and A n (0 x T*C^) show up in the mapping properties of P, 
namely whether it can be applied to all distributions, and whether it maps smooth 
functions to smooth functions. Thus, by Theorem |2.§| we would expect, if 5'^a(A) 
were an FIO, that the (twisted) wave front set of its kernel lies in 

(2.42) (7^^l_(i?_(A)) X 0) U (0 X 7^-V(i?+(A)) U Aq, 

where Aq C {T*C'^ \ 0) x {T*C'^ \ 0) is conic Lagrangian and projects to 7?./3q under 
the quotient map in both factors. While, as explained, our results as stated do not 



prove that the wave front set of the kernel of ^^^(A) lies in (2.42), we can prove its 
'quotiented' version. 



Theorem 2.10. Suppose that H is a many-body Hamiltonian, A ^ A, and (|2.34) 



holds. Suppose also that Ai is discrete. The wave front relation of Si3a{X) is 
a subset of the projection of a finite union of conic Lagrangian submanifolds of 
T*{C'g^ X C'li) \ 0, given by the broken bicharacteristic relation of H. Some of the 
Lagrangians may lie in r*C^ x or x T*Clj where denotes the zero section. 
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To see how this theorem is proved, recall that under our discreteness assump- 
tion the generalized broken bicharacteristics are actually broken bicharacteristics 
with a finite number of breaks (and a uniform bound on this finite number). Cor- 
respondingly, we can associate to each generalized broken bicharacteristic one of 
a finite number of collision patterns which describe from which clusters in which 
sequence the bicharacteristics reflect, and along which collision plane at which ki- 
netic energy they travel meanwhile. We prove in Appendix H that each of these 
patterns corresponds to a (typically not closed) conic Lagrangian submanifold of 
r*(C^ X C^) \ 0, the projection of the union of which gives the wave front relation 
(which is, however, closed), thereby proving this theorem. 

We remark that if there are no bound states in the proper subsystems, the 
proof of the result is essentially the same as in the three-body setting ^, ^; 
Hassell's proof employing Jacobi vector fields on the sphere is particularly easy 
(and transparent!) to adapt. 

The statement of this theorem shows that, under the discreteness assumption, 
the wave front relation of the S-matrices has the correct structure to be the wave 
front relation of a Fourier integral operator (FIO). However, the construction of 
such an FIO would require a much better understanding of the symplectic geometry 
of, and related analysis on, manifolds with corners that are equipped with certain 
boundary fibrations (corresponding to the resolved space [^;C] and the fibration 
given by the blow-down map which are discussed in the next section). In relatively 
simple settings this has been discussed in |g^, ^^ , but much progress is needed for 
dealing with the general case. 

3. Scattering geometry and analysis 

Next, we recall from |2^ Melrose's definition of the Lie algebra of 'scattering 
vector fields' VsdX), defined for every manifold with boundary X. Thus, 

(3.1) VsciX) ^ xVbiX) 

where Vh{X) is the set of smooth vector fields on X which are tangent to dX. If 
{x,yi, . . . ,2/„_i) are coordinates on X where x is a boundary defining function, 
then locally a basis of Vsc(-?) is given by 

(3.2) x^dx, xdy., j = 1,... ,n-l. 

Correspondingly, there is a vector bundle ^'^TX over X, called the scattering tangent 
bundle of X, such that Vsc(^) is the set of all smooth sections of ^'^TX: 

(3.3) Vsc(^) ^{ve C°°{X;'"'TX) : \fpeX, Vpe "%X}. 

The dual bundle of ^''TX (called the scattering cotangent bundle) is denoted by 
^'^T*X. Thus, covectors v E ^'^T*X, p near dX, can be written as 

,„ ., dx dv 

(3.4) y = r — +^j.--^. 

x^ X 

Hence, we have local coordinates (a;,y,r, /x) on ^'^T*X near dX. The scattering 
density bundle ^'^UX is the density bundle associated to ^'^T*X, so locally near 
dX it is spanned by a;~"~^ dxdy over C°°{X). Finally, Diffsc(-^) is the algebra of 
differential operators generated by the vector fields in Vsc(^); Diff^(-?) stands for 
scattering differential operators of order (at most) m. 
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An example is provided by the radial compactification of Euclidean space, X = 
§". We can use 'inverse' polar coordinates on R" to induce local coordinates on 
§" near 9S" as above. Thus, we let x = r~^ = \w\~^ (for r > 1, say, smoothed out 
near the origin), as in the introduction, write w — x^^to, lu G S"""'^, \w\ > 1, and 
let yj, j = 1, . . . ,n — 1, be local coordinates on S"^^. For example, one can take 
the yj to be n — 1 of the Wk/\w\, k = 1, . . . ,n. Then x g C°°(S!f:) is a boundary 
defining function of S" , and x and the yj give local coordinates near (9S!f: — §"~^. 

To establish the relationship between the scattering structure and the Euclidean 
scattering theory, we introduce local coordinates on X near p S dX as above, 
and use these to identify the coordinate neighborhood U oi p with a coordinate 
patch U' on the closed upper hemisphere S" (which is just a closed ball) near 
its boundary. Such an identification preserves the scattering structure since this 
structure is completely natural. We further identify §" with R" via the radial 
compactification RC as in ( |2.4| ). The constant coefficent vector fields d^j on R" 
lift under RC to give a basis of '''^TS" . Thus, V £ Vsc(S+) can be expressed as 
(ignoring the lifting in the notation) 

n 

(3.5) F = ^a,a»,, a,eC°°iSl). 

i=i 

As mentioned above, aj S C°°(§") is equivalent to requiring that RC* a^ is a 
classical (i.e. one-step polyhomogeneous) symbol of order on R". This description 
also shows that the positive Euclidean Laplacian, A, is an element of Diffgj,(§"), 
and that ''"^ilS" is spanned by the pull-back of the standard Euclidean density \dw\. 
In terms of the 'inverted' polar coordinates on R", covectors f • dw = ^ ■ ^j dwj 



take the form (3^) with 



(3.6) r = -^^-y^,r' + \^,\'^\e. 

\w\ 

Here /i is the orthogonal projection of ^ to the tangent space of the unit sphere 
S"~^ at y G §"~^, and |/i| denotes the length of a covector on §"^^ with respect to 
the standard metric h on the unit sphere. 

We next show that polyhomogeneous symbols on AT", pulled back to R" by tt", 
are smooth on the blown-up space [§" ; Co] ■ Recall that the blow-up process is 
simply an invariant way of introducing polar coordinates about a submanifold. A 
full description appears in |2^ and a more concise one in [pSl Appendix A] , but we 
give a brief summary here. Near Ca, S" can be identified with the inward-pointing 
normal bundle N'^Ca of Ca- Here N~^Ca is the image of the inward-pointing 
tangent bundle, T+S^, in the quotient bundle NCa = Tc^Sl/TCa, W G TqE^ is 
inward pointing if Wx{q) > 0. The blow up [§";Ca] of Ca in §" is then locally 
given by the blow up of N^Ca at the zero section, which amounts to introducing 
polar coordinates in its fibers. The front face of the blow up (the lift of Ca) is then 
identified with the inward pointing sphere bundle, S~^NCa, which is the quotient 
of N^Ca \ by the natural R"*" action in its fibers (0 denotes the zero section). 
We denote the blow-down map (which is C°°) by /3[§![;Ca] : [Sl;Ca] -^ §+. Now 
S~^NCa is a hemisphere bundle over Ca, which can be identified with the radial 
compactification of the normal bundle of Ca in S"'~^ whose fibers can in turn be 
identified with A". 
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To see this in more concrete terms, we proceed by finding local coordinates on 
[X; Ca] explicitly. Near Ca in S" we have \wa\ > c|w"| for some c > 0. Hence, near 
any point p G Ca one of the coordinate functions {wa)j which we take to be {wa)m 
for the sake of definiteness, satisfies |(w;a)m| > c'\{wa)j\, \{wa)m\ > c'|zi;°| for some 
c' > 0. Taking into account the coordinate form of RC we see that 

(3.7) 

x=\wa\ , Vj = -: 7- (j = 1, . . . , TO - 1), Zj = p (.7 = 1,... ,n-m) 

give coordinates on S" near p where we think of the j/j as coordinates on the 
unit sphere Ca- In these coordinates Ca is defined by a; = 0, z = 0. Sometimes 
we write these coordinates as Xa, Va, Za to emphasize the collision plane. Note 
that here the notation regarding y and z is reversed as compared to p4| . Wc call 
these the standard coordinates around Ca', we will always use these in this paper. 
Correspondingly, we have coordinates 

(3.8) X, yj {j = 1,... ,m-l), Zj = Zj/x (j = 1, . . . , n - m), 
i.e. 

(3.9) x^\wa\-\ ijj^^-^ (j = l,...,TO-l), Zj^iw'^), ij = l,...,n-m) 

near the interior of the front face ff of the blow-up [X;Ca], i.e. near the interior 
of ff = f3[X;Ca]*Ca- Similarly, one can easily write down local coordinates near 
the corner ffn(3[X; Ca]*dX, see Bj, Section 2]. As a result of such calculations, we 
deduce the following lemma and its corollary. 

Lemma 3.1. fpj. Lemma 2.5] j Suppose that X = §" and let j3 — I3[X] Ca] he the 
blow-down map. Then the pull-back /3*(RC~ )*7r'' o/ tt" : M" — > X"" extends to a 
C°° map, which we also denote by -k"" , 

(3.10) TT'':[X-Ca]^X\ 

Moreover, if x"" is a boundary defining function on X° (e.g. x"" — |w"|~^ for |w°| > 
1), then Pqx = (7r°)*a;° is a defining function for the lift of dX to [X;Ca], i.e. for 
13* dX. 

Corollary 3.2. f|34[ Corollary 2.6]; Suppose that X ^S^, f € S^f^^iX''). Then 

(3.11) (7rrfep,'xC^^{[X;Ca]). 

Here, following the previous lemma, we regard tt" as the map in (3.10|), and pgx is 



the defining function of /3[X]Ca]*dX , i.e. of the lift of dX, and the subscript phg 
refers to classical (one-step polyhomogeneous) symbols. 

This corollary shows that for a Euclidean many-body Hamiltonian, H = A -\- 
Sa^o, Va becomes a smooth function on the compact resolved space [§";Ca]. 
Thus, to understand H, we need to blow up all the Ca- The iterative construction 
was carried out in detail in [Q Section 2]; we refer to the discussion given there 
for details. However, we remind the reader that the Ca are blown up in the order 
of inclusion {opposite to the usual order on the clusters). That is, one starts with 
the blow-up of 2-clusters (which are disjoint); 3-clusters become disjoint upon this 
blow-up. One proceeds to blow-up the 3-clusters; 4-clusters become disjoint now. 
One proceeds this way, finally blowing up the N — 1-clusters. (The blow-up of the 
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A^-cluster is a diffeoniorphism, hence can be neglected.) We thus obtain a manifold 
with corners which is denoted by [S";C], and the blow-down map (which is C°°) is 
written as 

(3.12) (3sc = P[Sl;C]:[Sl;C]^Sl. 

The algebra of many-body differential operators is then defined as 

(3.13) Diffsc(§';,C) -C°°([S!^;C])®c-(S")Diffsc(S!^). 



That is, similarly to (|3J), P G Diff gj. (§*?:, C) means that 



(3.14) P= J2 ""^S, aaeC°°([§+;C]) 

|a|<m 

where we again ignored the pull-back by RC in the notation. 

The coordinates on ^'^T*X induced by the standard coordinates around Ca are 

.o 1 r\ dXa , dya dZa 

(3.15) Ta^Y^^^°■^ ^^°-' ; 

we included the subscript a to emphasize the element C = Cq of C around which 
the local coordinates are centered. Again, the roles of /i and v have been switched 
from l34|. In terms of the splitting M" = Xq = Xa © X"', and the corresponding 
splitting ^ = ($a,C°) of the dual coordinates, this gives 

(3.16) ^^ — Va, and, at Ca, Ta = r, ^a the orthogonal projection of ^a to TyCa- 

One of the main differences between Diffgc(§",C) and Diffsc(§") is that the 
former is not commutative to 'top weight'. That is, while for P E Diffg"(§"), 

Q e Diff™'(S!f.), we have [P,Q] G a; Diff™+'"'-^ (§!;:), this is replaced by [P,Q] e 
pco Diff™+'"'^i(S!^,C) for P e Diffg',(§'^,C), Q € Diff^J,'(§';,C) with pco a defining 
function for the lift of Cq. Thus, a vanishing factor (such as x above) is only 
present at the lift of the free face Co, i.e. there is no gain of a weight factor at 
the front faces ff. From an algebraic point of view, it is this non-commutativity 
that causes the breaks in the propagation of singularities by necessitating the use 
of the compressed bundle for wave front sets; one can only microlocalize in the 
commutative (tangential) variables, i.e. at Ca (where Za = 0) only in Xa, ya (i-e. 

Wa), and Ta, ^J,a (i-S. ^a)- 

For a closed embedded submanifold C of dX, the relative scattering tangent 
bundle ""^TiCiX) of C in X is the subbundle of ^'=TcX consisting of w G """TpX, 
p (z C, for which there exists 

(3.17) VeVsc{X;C)ciVsciX) 
with Vp = V (equality understood in ^'^TpX). Here 

(3.18) Vsc(^; C) = xVb{X; C) = x{V G Vb(X) : V is tangent to C} 

and tangency is defined using the (non-injective) inclusion map ^TX — > TX. 

Thus, in local coordinates (x, y, z) near p € C such that C is defined by a; = 0, 
z = 0, ^'^T{C; X) is spanned by x^dx and xdy- , j — 1, . . . ,m~l where n — to is the 
codimension on C in dX. Thus, in our Euclidean setting, X — §" , C = Ca = dXa, 
g the Euclidean metric, ^'^T{C; X) is naturally isomorphic to ^"^TcXa, i-e. it should 
be regarded as the bundle of scattering tangent vectors of the collision plane at 
infinity, spanned by 9(iu^) ., j — 1, . . . ,m, m — dimXa- 
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For C ^ Ca ^ C, the metric g defines the orthocomplement {^'^T(C\X))^ of 
'■''T{C; X) in """"TcX. The subbundle of """"T^X consisting of covectors that anni- 
hilate ("'=T(C;X))-L is denoted by '''=r*(C;X); it is called the relative scattering 
cotangent bundle of C in X. In our Euclidean setting, ^'^T*{Ca\X) is naturally 
isomorphic to ^^T^ Xa and is spanned by d{wa)j, j — I, . . . , to; so we simply write 

(3.19) ''T*{Ca;X)^''T^^Xa. 

In our standard local coordinates {x,y,z) near p £ C, C — Ca, C defined by 
a: = 0, z = 0, and xdz give an orthonormal basis of {^'^T{C;X))-^. Note that a 
basis of ^'^T{C; X) is given by x^d^ and xdy^ , while a basis of ^'^T*{C; X) is given 
by x~'^ dx, x~^ dyj. A covector in ^'^T*X can be written in these local coordinates 
as 

(3.20) ra^+^^.^ + ^^.^, 

•^a "^CL -^a 

Thus, local coordinates on ^'^TX-^X are given by (y, z, t, jjl, v), while on ^'^T*{Ca', X) 

by {y,r,v) = {ya,Ta,^J.a)■ 

Now if C = Ca, Cb e C with Ca C Ci,, we can further adjust our coordinates so 
that Cf, is defined by a; = 0, z" = 0, for some splitting z — (z' , z"). In fact, we can 
use the decomposition 

(3.21) R" = Xa e XI} e x''^ Xb" = X" n Xb, x"'' = x" n x^ 

the sum is of course still orthogonal with respect to g. Then the splitting z — (z', z") 

ra 



corresponds to the splitting X'^ = X§ ® X°^ in terms of the coordinates (3.7) 



Following ( |3.21| ), we write these coordinates as Za — [zab, ^a) 

With the corresponding splitting of the dual variable, v = [v' ,v") = {h'ab,i^a)j 
we obtain a well-defined projection 

(3.22) nta.'''T^'ACb;X)^"'T*{Ca;X), 

(3.23) ^TbaiO,ya,Ta,^J.a,l^ab) = {Va, Ta, ^J.a) ■ 

In our Euclidean setting this is just the obvious projection 

(3.24) TTba : ^'^Tgj^^Xb — > ^'^Tg-^^Xa 

under the inclusion Xa C Xb- We write tt for the collection of these maps. 

4. Indicial operators and the characteristic set 

We start by recalling from ||3J] the definition of the many-body pseudo-differential 
calculus via the quantization of symbols. Thus, the non-polyhomogeneous space 
^Scc(^+'^) i^ ^be following. Suppose that, with the notation of p2| , 

(4.1) ae^-"''([§;;C]x§!^). 

That is, identifying int(S") and int([S" ;C]) with R" as usual (via RC~^), suppose 
that a e C°°(M2, x R^) has the following property. For every P e Diffb(§'J:), acting 

on the second factor of §" (i.e. in the ^ variable), and Q e Diffj^ ([§";C]), acting 
on the first factor of §" (i.e. in the w variable), fc, k' e N, 

(4.2) PQaepZ"'p'9L°^iK^K) 
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where poo and pQ are defining functions of the first and second factors of §!f respec- 
tively. Let A = qL{o-) denote the left quantization of a: 

(4.3) Au{w) = (27r)-" / e*('"-"'')-«a(w, O^iw') dw' d^, 

understood as an oscillatory integral. Then A e ^^g"';, (§+, C). We could have equally 
well used other (right, Weyl, etc.) quantizations as well. 

In addition, the polyhomogeneous class ^gj,' (§",C) is given by the quantization 
of symbols 

Since differential operators ^ aa{w)D°' are just the left quantization of the symbols 
a{w,S,) — '^aa{w)£,°', it follows immediately that 

(4.5) Diff^,(§'^,C)cvI/g^,(S!^,C). 

It was shown in pj] that ^scc(§",C) and 'I'sc(§"7C) are ^-algebras with re- 
spect to composition and taking adjoints. Moreover, \1/^'^(S",C) is bounded from 
Hlf{^X) to Hl-'^-^+^Wl)] this follows from the inclusion *™^S!f:,C) C ^™(M") 
where the latter is Hormander's 'uniform' calculus in M", see [|2[ Section 18.1]. 

In the standard pseudo-differential calculus on compacts manifolds without bound- 
ary Fredholm properties are captured by the principal symbol. In our setting, 
Hormander's principal symbol map on ^™ (R") restricts to a principal symbol map 

(4.6) asc,™ : *"/(§';, C) ^ 5r('^r*[§;;C]), 

'S'h(^'^r*[§" ;C]) denoting the space of smooth symbols which are homogeneous of 
degree m on s<=T* [§!;:; C] , the pull-back of ^''T*E>\ by the blow-down map. How- 
ever, the invertibility of the principal symbol is not sufficient for an operator to be 
Fredholm, since the behavior near infinity must be taken into account as well. The 
additional piece of information that is necessary is the invertibility of the indicial 
operators whose definition we briefly recall below. First, however, we remind the 
reader what the structure of the lift of Ca to [X\C] is. 
For Ca&C let 

(4.7) Ca = {aeC: CbCCa), 

(4.8) C^^ICbGC: Ca^Ci,}. 

We carry out the blow-up [§" ;C] by first blowing up Ca- Since all elements of Ca 
are p-submanifolds (i.e. product submanifolds) of Cq, the lift (3^^\Ca]*Ca of Ca to 
[S";Ca] is naturally diffeomorphic to 

(4.9) Ca^[Ca\Ca]. 

Thus, over C^, the regular part of Ca, Ca can be identified with Ca- The front face of 
the new blow-up, i.e. of the blow up of /3[§" ;Ca]*Ca in [§" ;Ca] is thus a hemisphere 
(i.e. ball) bundle over Ca, namely S^NCa- We write the bundle projection, which 
is just the restriction of the new blow-down map to the front face, S~^NCa as 

(4.10) Pa : S+NCa -^ Ca. 

In our Euclidean setting, these fibers can be naturally identified with X°- via the 



projection tt" (extended as in Lemma 3.1) 
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Every remaining blow up in [S" ; C] concerns submanifolds that are either disjoint 
from this new front face or are the Hft of elements of C". The former do not 
affect the structure near the new front face, S^NCa — /?[§+; Ca; Ca]*Co, while the 
latter, which are given by the lifts of elements of C, correspond to blow ups that 
can be performed in the fibers of S'^NCa- Note that the lift of Cb E C°, meets 
the new front face only at its boundary since all Cb are subsets of S"^^ — 9§" . 
In particular, the hft /3g(.C'a of Ca to [§!j:;C] fibers over Ca and the fibers are 
diffeomorphic to a hemisphere (i.e. ball) with certain boundary submanifolds blown 
up. More specifically, the intersection of /3[§" ;Co; Ca]*Cb, Cb S C", with the front 
face S~^NCa is the image of T/3[§"; Ca]*Cb under the quotients defining the spherical 
normal bundle; (3g^Ca is obtained by blowing these up in S^NCa- Hence, the fiber 
of /3g,Ca over p e Ca is given by [S+NgCalTgC] where q = /3[§!^;Ca](p) G Ca- In 
particular, in our Euclidean setting, the fibers of Pg^Ca over Ca can be naturally 
identified with [X'^;C''] via Tr*". 

We now define "'=T*(C'a; §!f:) denote the pull-back of '''=r*(Ca;§!^) by the blow- 
down map P[Ca;Ca]- 

(4.11) '''T*{Ca;Sl) ^ f3[Ca;Car'T*{Ca;Sl). 

Thus, ^'^T*{Ca'-,S^) is, via the metric g, naturally diffeomorphic to Ca x Xa, where 
the metric is used to identify Xa with its dual. If Ca C Cb then nba lifts to a map 

(4.12) iTba : ''T;^cM'cSCb;X) ^ ^^r*(C„;X). 

For each operator A E 'i^^ {X,C), the Ca-indicial operator of A, denoted by 
Aa,i, will be a collection of operators, one for each C, E ^'^T*{Ca\X)^ acting on 
functions on the fiber pa^{p) of Pa- Recall that the interior of these fibers can be 
naturally identified with X"-. So suppose that u S C°°(p^^(p)), i.e. u is a Schwartz 
function on the the fiber above p which is a compactification of X"; we need to 
define Aa{C}u. For this purpose choose / G C°°(§";R) such that d{f/x), evaluated 
at P[Ca;Ca]ip), is_ equal to C- Then let A = e-'^/'^x-'Ae'^^/'' e *™/(X,C), and 
choose u' e C°°{[X;C]) such that M'lp-i(p) = u. Then 

(4.13) i,,(C)^i-(i«')lp-(p), 

which is independent of all the choices we made. 

In our Euclidean setting we can simply take / to be the pull-back of a function 
in C°°{Xa) (which we write, abusing the notation, as / <E C°°{Xa))- Thus, at least 
in the regular part C^ of Ca, the a-indicial operator of A arises by thinking of A 
as a scattering pseudo-differential operator on Xa with values in operators on X°-, 
and finding its scattering symbol at C'a C Ca = dXa- In fact, Aa,i{£,) is a smooth 
function of ^ G ^'^T*{Ca;X), so the above description also determines Aa^i at dCa- 
We often simplify (and thereby abuse) the notation and drop the index I, i.e. we 
write Aa = Aa.i, when the value of I is understood. 

In the case of Euclidean many-body scattering, Aad is a function on (3'^^'^T^ Xa 
with values in operators on S{X°')] here 

(4.14) Pa = l3[Ca;Ca] : Ca = [C, ; C J ^ Ca 
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is the blow-down map. Thus, Pa is simply the restriction of l3[Xa]Ca] to the lift 
Ca^f3[Xa;Ca]*Ca. Then 

(4.15) Aaj £ C°°(/3r^^g*^^Xa,^P™/(X^r)). 
We write the space on the right as 

(4.16) Aa,i e ^';:l^,^{f3rT*{Ca;X),Ca). 



An explicit formula for the indicial operator of A, given by (4.3), at Ca is 
(4.17) 



Note that the oscillatory testing definition that we adopted automatically shows 
that the indicial operator map is multiplicative. 

The principal symbol and the indicial operators together indeed describe com- 
pactness properties of ps.d.o's. In particular, if A e *g"' (§",C), asc,m{A) never 
vanishes, and for every a, Aa^i is invertible in ^^' , (/3*^'^T*(C'a; X),Ca), then there 

is a two-sided parametrix P for A such that PA - Id, AP - Id e *g^°°'°°(S!J:,C). 
For A € VE*^' (S",C) self-adjoint, m > 0, such that asc,miA) never vanishes (but 
no assumptions made on the indicial operators of A) , A £ C \ M, we conclude that 
{A - A)-i e *g™^°(S!^,C). The invertibihty of the indicial operators of A - A is 

automatic since (A — A)^ g = Aafi — A, and then an iterative argument, starting at 
a = 0, using the self-adjointness of A, shows the required invertibihty. In addition, 
an argument using almost analytic extensions and the Cauchy representation also 
shows that certain functions of self-adjoint ps.d.o's are ps.d.o's themselves. In par- 
ticular, for A g Vl/g"' (S" ,C) self-adjoint, m > 0, such that crsc,m{A) never vanishes, 
(f> £ C^(R), we conclude that (f>{A) £ *g;^'"(§!J:,C) (see Q, Propositions 4.7-4.8], 
and see |^ and ||ll| for a general discussion of functions of ps.d.o.'s). 

Now it is particularly easy to identify the indicial operators: for each ^ — 
{y,T,fi) £ ^'^T*{Ca]X), they are operators on X"-, identified as the fiber of the 
front face Pa^ip) over p. A simple calculation, see p% Sections 4 and 11] for more 
details, shows that the indicial operators of H are given by 

(4.18) i?a,o(0-^a,o((p,0)) + |eP, ^={y,T,fi)£^-T*{Ca;X), 

(4.19) Ha,o{P,0) = Az + Vip,Z) = H'' 

where Z are Euclidean coordinates on the interior of p~^{p), i.e. on X°, and Az is 
the Euclidean Laplacian. Note that under these natural identifications (that p~^{p) 
is regarded as X°- for all p) the a-indicial operators of H only depend on |^p, i.e. 
on the metric function of ga- 

Equations ( 4.18| )-( 4.19| ) show that Ha,o{p,0) is uniformly bounded below, so for 



any tp GC^(R) the set 

(4.20) UaCm£''T*{Ca;X): ^{Ha{0) ^ 0}) 

is compact. 

Note that iJa, o(0-^ is invertible in *sc(X°,C°) if and only if A- l^aP is not in 
the spectrum of PP^. Indeed, by the HVZ theorem, this means that A — |^aP is not 
an i^ eigenvalue of H"^ , and it is smaller than all of its thresholds , A — | ^a P < inf A^ , 
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i.e. smaller than the eigenvalues of subsystems of -ff". By an iterative argument 
corresponding to Ca C Cb, A — |^aP < inf Aa guarantees the invertibility of the 
indicial operators of Ha.oiO ~ ^^ ^^ Ha.oiO — A is Fredholm, and hence invertible if 
and only if it has no L^ eigenf unctions. Thus, we define the characteristic variety 
S(A) C ^'^T*X of iJ — A to be the union image of these sets under tt: 

(4.21) ±{X)^Ua±aiX)Cf*X, 

(4.22) 
Sa(A) ={Ca e '''T^,iCa;X) : A - l^al' > <J for some a e speCppiJ^ Ca C Ct} 

U {(a e ''T^,{Ca;X) ■■ A- iCaP e specppiin- 

If there are no bound states in the subsystems, a compactness argument as in [|^ 
allows one to show that the free Laplacian governs the propagation of singularities 
except that breaks occur at the collision planes Ca where the usual law of reflection 
is satisfied. This indicates that one should think of the characteristic variety S(A) 
as the union of projection of the characteristic sets corresponding to the bound 
states of the various subsystems. Thus, let 

(4.23) I]fc(A) = {Cb e '^T^'.Xt : A ~ |6|' e spec^^ H"}. 

For a with C'a C Cb, the projection of Sh(A) under nba is 

(4.24) 

TTbai^bW) = {Ca S '''T*^Xa I 3^ G '''T*^Xb with TTfeJCfc) = ^a and A - 1^6^ G SpCC^^ H''} 

= {Ca e '"T^^Xa : A - l^aP > c7 for some a G spcc^^ H''}. 
Thus, 

(4.25) S,(A) =Ucoc„7r,a(E,(A))n^^T5,(a;X). 

In fact, letting vr;, : ^"^T^ Xb — > ^'^T*X be the projection, we thus conclude that 

(4.26) i](A)=U7rb(Efc(A)). 

We write ttj, for the restriction of ivb to ^^(A). 

We proceed to recall the basic topological properties of ^'^T*X from |3J, Sec- 
tion 5]. We put the topology induced by tt : '"'T*^X -> ^'^t*X on '"'t*X. Thus, 
C C ^'^T*X is closed if and only if tt^^{C) is closed, so if / is continuous on 
^'^T*^X and 7r-invariant, then /tt is continuous on ^'^T*X. In particular, there 
are always continuous functions separating points on ^'^T*X: ii p{£,) 7^ p{£.') (here 
p : ^'^T*X -^ dX stands for projection to the base), one can use the pull-back of an 
appropriate function on dX, and if p(0 = p{^') & C'a, then ^,^' € ^'^T^,Xa arc of 

the form (j7a,Ca) and {y'a,Ca)j ^a 7^ fa J ^'^d the function ^ f-» ^a is well-defined and 
TT- invariant on a neighborhood of p(f ), so multiplying it by the pull-back of a cutoff 
on dX gives a globally well-defined separating continuous function. Thus, ^'^T*X 
is Hausdorff. 

By the continuity of n, ii K C ^'^TXi^X is compact, the same holds for 7t{K). 



But by ( [4.26D , i;(A) = 7r(Ubif,(Sfc(A))) (see (p.2lD), and tbi^bW) is compact for 
each 6, so S(A) is compact, hence closed in '^'^T*X. In particular C C S(A) is closed 
in S](A) if and only if it is closed in ^'^T*X, i.e. if and only if n^^{C) is closed in 
^''T*^X. Also, if C C t,{X) is closed, then TrQ^{C) is compact. 
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Note that for all C C t.{X) and for any R> with i? > A - inf Ai, 

(4.27) ^(Tr-'iC)) = n{n-\C) n U G ^^r^*^X : \^\' < R}). 

Indeed, if C e 7r-i(C), 7r(0 e "'^T^,^,,, then by C C S(A), |ia(7r(0)P < R. 
7r(ia(7r(^))) = 7r(^), so 7r(ia(7r(^))) is in the right hand side, showing one inclusion, 
and the other is clear. For any b then 

(4.28) n^\C) = TT,t{7T-\C) n U : l^l' < R} n ''T*c,X) n Eb(A), 

and TTob is continuous, so ■k'^ (C) is also compact (and similarly tt^^ (C) is closed). 
Thus, all the maps tt;,, TTb, are continuous. 

Also, fix ^ S ^'^T^iXa, write ^ = (VaAa), and choose a neighborhood C/q of 
y„ = p(^) in ax such that CTq n Cb = unless Ca C Ct- Let w^- = w : "'^T^^ ^ K 
be given by the following 7r-invariant function on ^'^T^X (also denoted by lu): 

(4.29) LoiO = (j/a - ya)' + \Za\^ + (Ca - ^a)' > 

in the coordinates {ya,Za,ia,^'^)- Suppose that U is open in S(A), f G f/. Thus, 
iiT = (S(A) n^'^Kt-X) \ 17 is compact, so, unless K is empty, ui assumes a minimum 
on it which thus has to be non-negative. But U!{(,) — implies Za = 0, so p(^) G Ca, 
and then ya = Va, ^a — Ca show that £,=£,■ Since ^ € [/, this shows that there 
exists S > such that cu > S on K. Replacing 5 > by possibly a smaller number, 
we can also assume that uj{£_) < S imphes p{£_) £ Uo- We thus conclude that if U is 
a neighborhood of f , then there exists 6 > such that 

(4.30) {e e S(A) : tj(0 < 5} C [/. 

Of course, we could have equally well used {Ta,Ha) in place of ^a- These sets are 
open since u is continuous, hence they form a basis for the topology as ^ and S vary; 
it is easy to see that if one restricts both of these to suitable countable sets, one 
still has a basis. Note that, separating the complement of U from ^ by a level set of 
w shows explicitly that S(A) is regular, and a simple compactness argument using 
these ujf (composed with cut-off functions on the reals as in the next paragraph) 
shows that S(A) is normal, hence it is a compact metrizable space. 

Composing ui with a C°° function on R supported near also shows that given 
any ^ S S(A) and any neighborhood [/ of ^ in S(A), one can construct a 7r-invariant 
C°° function / on ^^T*j^X for which UiO ¥" and supp /tt C U. This also shows 
the existence of smooth partitions of unity on S(A), smoothness understood as 
smoothness for the pull-back to ^'^T*X. 

If we define 

(4.31) i:(/) = UAe/E(A) 

for / C R compact, all the previous statements remain valid with trivial modifica- 
tions with S(/) in place of E(A). We write 

(4.32) Sb(/) = UAe/S6(A), ntj = 7rb|s,(/). 

Thus, S](/)=Ub^6(I]6(/)), 

(4.33) i e vr-;({f}) iff TTbii) = e and 3a € speCp^iH'') such that I^T +<JeI. 

We end this section by recalhng from |^, Section 5] the definition of the operator 
wave front set and the wave front set of distributions in the many-body setting. 
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Both of these are closed subsets of ^'^T*S^^ which were modelled on the fibred cusp 
wave front set introduced by Mazzeo and Melrose in [gOJ. Rather than using the 
invariant definition given there in terms of oscillatory testing, we recall the one 
given in terms of representation of ps.d.o.'s via quantization. 

So suppose that A € ^gc°°' (S+jC) is the left quantization of a symbol a £ 
p^gC°°{[E>l;C] X SI) which vanishes to infinite order at [E>l;C] x dS^. Then ^ ^ 
WF'sJa), ^ G ^''T*{Ca;X), peC'^,ii and only if there exists a neighborhood U of 
£, in «'=T*§» such that a vanishes at U' C (a[§'|;C]) x E" to infinite order where U' 
is the inverse image of U under the composite map 

(4.34) ^ : (c)[S:^;C]) x M" -B^^^^ (08%) x M" == ^'=rg*„_i§'^ — "^^ ^'"t*^,^. 

It follows immediately from the usual formulae relating quantizations that this 
definition is independent of such choices. For example, we could have equally well 
written A as the right quantization of a symbol with similar properties. 

The general definition for A g \I>^' {X,C), again following the paper |g^, in 
the explicit quantization form as in the previous paragraph, would also require 
the rapid decay of a in an open cone (conic in the cotangent variable, ^, i.e. in 
the second factor, R", in ( [4.34D ) that includes U' . For A € *g^°°''(X,C), a is 



rapidly decreasing in every direction as |^| —^ oo, so this statement is vacuous. 



and we recover the above definition. The main point is that if ^4 G ^tg^ {X,C), 
Aa{£,) is invertible, then there exists a microlocal parametrix for A, i.e. there exists 
G e *°f (X,C) such that Id = AG + Rr, Id = GA + Rl, with Rr,Rl € *°f (X,C), 
C ^ WFsc(i?fl), C ^ WFsc(i?L); see |^, Lemmas 14-15]. 

Werecall that the many-body wave front set WFsc(w) ofu e C^°°{X), defined in 
|34[ Section 5], is a subset oi^'^T*X, defined as follows. The definition is somewhat 
delicate due to the relation among the principal symbols crsc,o(^a(0) ^ ^ varies, 
see [p4| for a detailed discussion, but it becomes much simpler for approximate 
generalized eigenfunctions of H as the proposition following the definition shows. 
Although the definition is complicated, it ensures that all the desired properties 
hold: for aU u G C-°°{X), WFsc(u) is closed in "'=t*X, 

(4.35) 

Ul, U2 e C-°^\X) ^ WFsc(wi + U2) C WFsc(ui) U WFsc(w2), 

u e C-°^iX), A e r^/{X,C) ^ WFsdAu) c WFsc(u), 
A e *s;^''(X,C), u e C~°°{X), WFg^(A) compact, WFs^(A) n WFsc(u) = 
^Aue C°^\X). 
The central one of these means that '^gj,*(X,C) is microlocal'. 
Definition 4.1. ([|4[ Definition 5.2]) We say that 
(4.36) 

i i WFsc(u) n ^^T^, (Ca;X) iff 3A e *gf (X,C), KM) invertible in *gf (X",C") 

3i?, evI;g-."(X,C), e^WF^,(B,), 

3u, eC-°°(X), j-l,...,s, /eC°°(X), 

s 
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Similarly, the filtered version of the Sc-wave front set is given by 

(4.37) 

e i WF™''(u) n'^^Ta, (Ca;X) ifFBA e *gf (X,C), i„,o(0 invertible in *°f (X",^^), 

3u, ec-°°(x), j = i,...,s, / e i?3'^''(X), 

s 

The wave front set has the following property which is useful for approximate 
generalized eigenfunctions of H. 

Proposition 4.2. (^, Proposition 5.5]j Suppose that u € C~°°{X), A G M, and 
define W C ^'=r*X by 

(4.38) 
e iWn''T^,{Ca;X) tjf3^j e C^{R), V(A) - 1, 

3A e vi,^-.0(x,C), AaiO = Wl)aiO.Au e C°°(X). 

Then 

(4.39) WF5c(w) C WFsc((ff - A)w) U W. 

T/ie same conclusion holds with WF^^ replaced by WF™J and Au G C°°{X) by 

AueHZ'\X). 

Corollary 4.3. (Microlocal elliptic regularity.) Suppose that u G C^°°{X), A G M. 
Then WFsciu) C WFsdiH ~ \)u) U ±{X). 

Proof. If Ca ^ S(A), ^a e '''Tlj,Xa, One can choose V G Cf (R) such that ?A(A) = 

1 and i^{H)^{£,a) = ipiH"- + |faP) = by making the support of -0 sufficiently 
small. Then we can take A to be the zero operator, showing that f^ ^ W. Thus, 
Proposition |4.2| proves the corollary. D 

5. Broken bicharacteristics 

Since at C^ a particle can be in a bound state of any of the clusters Cb with 
Ca C Cb, it can be expected that the Hamilton vector fields associated to the b- 
external kinetic energy, i.e. to Ax^ , govern the propagation of singularities. If the 
particle is in a &-bound state of energy e^, its 6-external kinetic energy is A — e^. The 
symbol of Ax,, — (A — ep) at Cb — dXb (i.e. its sc-indicial operator) is |^f,p + €p — X. 
The Hamilton vector field of gb = \£,b\'^ (or \^b\^ + e/3 — A) on T*Xb is 

(5.1) Hl^2Cb-d^,. 

Following the general principle that it is more convenient to do analysis on compact 
spaces than analyzing uniform properties in a non-compact setting, we rescale H^ 
and regard the result as a smooth vector field on ^'^T*Xb, i-e. as an element of 
V(''^r*Xfc). Thus, the rescaled Hamihon vector field 

(5.2) -i/^^ = x-'H'g 
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of ffb = l^hPj introduced in |g3[, is 

(5.3) 

= 2n{xbd,, + fib ■ d^J - 2|/ib|'9r, + H|^,|2 + XbW, W e V^rT*Xb), 
so its restriction to Cb, also denoted by ^"^H^, is 

(5.4) ^'H'g - 2Tbfib ■ d^, - 2\fib\^dr, + i?|^,|.. 
Thus, for CaClCb, 

(5.5) ^''H'g = 2x^\^a ■ d^^ + U ■ a»„J = {xa/xb)rH;; + 2U ■ d^J. 
If fe = 0, then in agreement with p3| we write 

(5.6) ^""Hg = ^^=77°. 

Before defining generalized broken bicharacteristics, we discuss how the rescaled 
Hamilton vector fields are related for various b. First, note that ^"^TX X is the 
bundle direct sum of ^'^T^^Xb and the annihilator of ^'^Tc^{Cb]X). Let 

(5.7) lb : ^'T^Xb -- ^^ra^X 

be the inclusion map as the zero section in the second sunimand, so 

(5.8) ib{yb, n, fib) = ivb, 0, n, fib, 0), 
i.e. VbMO) = for all ^ e """T^^Xb- Then 

(5.9) M*''H'g = ^'Hg, 

since the corresponding statement for the non-rescaled Hamilton vector fields is 
easy to see and Xb/x is C°° near Cb, equal to 1 at Cb- 
Next, note that 

(5.10) T*K" = T*Xa X r*X". 

Thus, if / is a smooth function on '^'^T^ Xa, we can extend it to a smooth function 
on ^'^T*Xa using polar coordinates (i.e. so that it is independent of Xa, at least 
near Ca), and pull-back the result to T*R" using the projection to the first factor. 
The result is then smooth on >5cy*gn jjj g^ neighborhood of Ca (in fact, everywhere 
away from the closure of X"" in S" ). Thus, we have a pull-back, which we denote 
as F = {ttD* f. Note that F = (ttJ;)*/ is 7r-invariant near C^. Moreover, 

(5.11) """RgF ^ — «)*"'= iJ^/. 

One often drops the pull-back notation, and simply writes / in place of F. A similar 
formula also holds for ^'^TJ^/ by (5.9). 

It is worth calculating the derivatives of two functions in particular along the 
rescaled Hamilton vector fields. First, 

(5.12) ^^HgT = -2\fi\^, 
so by (IH), 

(5.13) -H^V = ^2\fLb\^. 
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Next, consider 

(5.14) rja = Za ■ Va ^ —r 



Then 

(5.15) -^^^^^(2r„,y, + 2|rn, 

X 

so 

(5.16) ^^if^V = —{2TaVa + 2|C?n - — (2t,77, + 2\iya. 

" X X 



m') 



We define generalized broken bicharacteristics of H — X partly following Lebeau 



19 1, but in such a way that all the analytic properties (compactness, applicability 
of positive commutator estimates) will be clear. However, the geometric properties 
will be less apparent, and we will devote some time to clarifying these. 

First, we say that a function / e C°°{'"'T*-^X) is 7r-invariant if for f, |' e """Tg^X, 
7r(^) = 7r(^') implies /(^) — f{C)- A tt- invariant function / naturally defines 
a function /^ on ^'^T*X by /7r(0 = fiO where ^ G ^'^Tg^X is chosen so that 

We also need to introduce notation for the upper and lower one-sided derivatives 
of functions / defined on an interval /: 

fit) - f{to) 



(i?+/)(to)=liminf 
(5.17) 

(I?-/)(to)-limsup 



(i^_/)(^o)=liminfM^, 

(i?+/)(to)=limsupM:i/M, 
t^to+ t — to 

fit) - fito) 



t-to 

A function / is thus differentiable at t^ if and only if these four assume the same 
(finite) value; then we write these as Df{to). Similarly, / is differentiable from the 
left at to if and only if (D_/)(io) = (-D~/)(io) is finite; we write the common value 
as Df{to-). 



Definition. (Definition p.l|) A generalized broken bicharacteristic of iJ — A is a 
continuous map 7 : / —> ^(-^)5 where / C M is an interval, such that for all to (z I 
and for each sign + and — the following holds. Let ^0 = 7(^0), suppose that 
Co e ''=T'5,Xa. Then for all 7r-invariant functions / e C°°("'^T^^X), 

(5.18) D±iU o j)(to) > Mi'^H'ji^o) ■■ Co e T^b^H^o), ^ C a}. 

In this section we usually drop the word 'generalized' for the sake of brevity. 

Remark 5.1. Considering — / in place of /, the definition immediately gives a sim- 
ilar estimate for D^{fT^ o 7) (to). In addition, the function / only has to be defined 
on '"'TIjX where [/ is a neighborhood of p in dX, 7(^0) e ^''T*{Ca;X). Indeed, 
otherwise one can consider 0/, where cj) is (the pull-back of) a cut-off function sup- 
ported in U that is identically 1 near p, without affecting any of the statements 
above. 
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Remark 5.2. The local coordinates j/a and ^a are actually differentiable along 7 at 
to since the Hamilton vector fields ^'^H^ applied to them give 7r-invariant results 



which agree with each other by (5^) and ( 5.11 ). In addition, over Cq, tt is one-to- 
one, so if J is an interval and 7|j lies in ^"^T^iX then 7|j is a bicharacteristic of 
-if,. 

Moreover, since Uf,'iT^'^(S(A)) is compact, it follows that for any 7r-invariant func- 
tion /, the derivatives D±{f.,^ o 7), D^{fTj o 7) are bounded (independently of to, 
and even of 7), hence /^r ° 7 is Lipschitz; in fact, uniformly so (i.e. the Lipschitz 
constant depends on /, but not on 7). This at once implies the equicontinuity of 
the set of broken bicharacteristics as in (l9| . 

Lemma 5.3. (cf. Lebeau's proof, |19| , Proposition 6]) For any compact interval 
I — [Ti,T2\, the set TZ of broken bicharacteristics "f : I —> S(A) is equicontinuous. 

Proof. Let d be a metric defining the topology on S(A). Suppose TZ is not equicon- 
tinuous. That is, suppose that there exists eq > 0, sequences Sk, s'f, in /, 7^ e TZ, 
such that \sk — s'f.\ < 1/k but 'i(7fc(sfe),7fc(Sfc)) > eo- Since S(A) is compact, 
one can pass to subsequences (which we do not show in the notation) such that 
7fe(sfe) and 7fe(s'j.) converge to points <^ and <^' respectively. Note that £, ^ C 
since d(^,C') ^ eo- For any 7r-invariant function /, fj^ o 7^. is uniformly Lips- 
chitz (i.e. the Lipschitz constant is independent of k), so \sk — s'f.\ < 1/k shows 
that |/^(7j.(sj,)) — fnilk{s'k))\ < M/k; M independent of k. But /^ is continu- 
ous, so limfc^oo/77(7fc(sfc)) = fniO^ limfc_,oo/7r(7fe(Sfe)) = UiC), so we conclude 
/^(^) = /^(^')- But let a be such that ^ S '"'T*Xa, p € C'^. All functions on 
X, pulled back by the bundle projection to ^'^T*X, are 7r-invariant, so we see that 
^1 g s'=T*Xa as well. But {£,a)j is 7r-invariant near ^ for all j, so this shows £, ~ £,', 
a contradiction. D 

Next, we note that the uniform limit of broken bicharacteristics is a broken 
bicharacteristic . 

Proposition 5.4. Suppose that I is a compact interval I ~ [Ti, T2\, "/k '■ I ^ "^W 
are broken bicharacteristics, 7 : / — > S(A), and 7^ ^ 7 uniformly on I. Then 7 is 
a broken bicharacteristic. 

Proof. Let / be a 7r-invariant function, to G /, ^0 = 7(^0), and let 

(5.19) CO = inf{-if,V(lo) : Co e 7r-i(eo) : Ca C a}. 

We need to show that D±{f.^ o 7)(io) > cq. Wc only consider D+ for the sake of 
definiteness. So let e > 0; we need to prove that there exists S > such that for all 
te {to,to + S), f^o-f{t)-f^oj{to) > (co-e)(t-tn). 

But by the continuity of ^'^Hgf on ^'^Tq X\j, there exists a neighborhood U of 

Co in E(A) such that """H^gfii) > cq - e/3 for aU b, i £ K^iO^ C e C/. Next, by 
the uniform convergence of the 7^, there exist 6 > and M E N such that for 
tE{to-6,to + S),k> M, 7fc(i) G U. Let 

(5.20) Fk-.I^R, Fk{t) = f^o^k{t)-{co-e/3)t. 

Thus, for t G {to,to + S), k > M, Fk satisfies D+Fk > 0, hence Fk is non-decreasing, 
so for t G (to, to + S), k> M, 

(5.21) /, o 7fc(i) - /, o 7fe(to) > (co - e/3)(i - to). 



MANY-BODY SCATTERING 31 

Now given t e {to, to + 6), simply choose k > M such that \f-^ o jk{t) — /tt ° 7(01 < 
e{t - io)/3, \U ° 7fc(io) - fn ° 7(io)| < e(i - io)/3, which is possible since /^ is 
continuous, so /^ o 7fe(i') -^ fn ° l{i') for all i'. Thus, using this particular value of 
k and ( |3.2l[ ), we conclude that 

(5.22) /^ o 7(t) - /, o 7(io) > (co - e)(t - to)- 

This holds for every i e (to, ^o + ^)j completing the proof. D 

We thus deduce the compactness of the set of broken bicharacteristics. 

Proposition 5.5. For any compact interval I = \ri,T2], the set TZ of broken 
bicharacteristics 7 : / ^ ^(-^) 's compact in the topology of uniform convergence. 

Proof. Since TZ is equicontinuous and S(A) is compact, any sequence of broken 
bicharacteristics in TZ has a convergent subsequence, converging uniformly over /, 
by the theorem of Ascoli-Arzela. But the uniform limit of broken bicharacteristics 
is a broken bicharacteristic, proving the proposition. D 

Corollary 5.6. (cf. Lebeau, fig, CoroUaire 7]) If ^ : (Ti,T2) —> S(A) is a general- 
ized broken bicharacteristic then 7 extends to a generalized broken bicharacteristic 
on [Ti,T2]. 

Proof Let 7„ = [Ti, Ti + S] ^ t{X) be given by 7„(i) = -f{t + 1/n), where (5 > is 
chosen sufficiently small. By the previous proposition, there exists a subsequence 
of {7n} that converges to a broken bicharacteristic 7 : [Ti,Ti + S] -^ 5](A). But 
then 7|(Ti,Ti+i5] = 7, so 7 gives the desired extension of 7 to [Ti,T2). The other 
endpoint can be dealt with similarly. D 

Remark 5.7. This corollary shows that every generalized broken bicharacteristic 
7o : J — > E(A) can be extended to another one, 7, defined over R (meaning that 
7 : M ^ 'S{X) and "/\j = 70). Indeed, by the corollary we may assume that J 
is closed. It suffices to show that, say, 7 can be extended to an interval that is 
unbounded from below; the extension in the other direction is similar. So suppose 
the lower endpoint of J is Ti G R (if it is — cx), we are done), and suppose that 
7o(ri) e ^''T^.Xa n Trba{^b{Xj)- such a b exists by ( |425| ). Then choose ^ € Eb(A) 



such that TTbaiO = 7o(Ti), and define 7|(-oo.Ti]uj by ^(t) = 7r(exp((i-ri)-i?g^)(0) 



for t < Ti, 7(i) — 70 (t) for t £ J . Directly from Definition (2T), 7 is a generalized 
broken bicharacteristic with the desired properties. 

Since Lipschitz functions are differentiable almost everywhere, and they are equal 
to the integral of their a.e. defined derivative, we can analyze 7r-invariant functions 
such as r in more detail. 

In fact, the Hamilton vector field ^'^H^ applied to — r gives 2|/ifcp > 0, so we 
deduce that — r is monotone increasing, i.e. r is monotone decreasing along broken 
bicharacteristics. In fact, writing F = — r o 7, note that ^ = 7r^^('C) being in Sb(A) 
gives 

(5.23) A - llfcp = A - r2 - j^^p e speCpp(ij''). 

Thus, A — r^(7(to)) ^ spcCpp(iJ'') for any b with Cb D Ca imphes that 

(5.24) D±F{to) > Co > 0, 



32 ANDRAS VASY 

SO there exist c > and 5 > such that 

(5.25) |r(7(t)) - T(7(io))| > c\t - to\ if \t - io| < 5, 

so the bicharacteristic 7 is not constant. Conversely, if there exists b with Cb D Ca 
such that A — r^ (7 (to)) G spec (iJ''), then the constant curve 



(5.26) 70 : K -> ^'=r*X, 7o(i) = 7(^0) for all t e 

^^(7(^0)), we hav^ po - w, ow ^^g 



is a broken bicharacteristic since for ^ G vr^ ^(7(^0)), we have fib — 0, so ^'^iJ^(^) 



vanishes. 

Thus, we define the 6-radial sets 

(5.27) -R^(A) = {6 e nW ■■ Mfa = 0, ±T > 0}, 
and their image by 

(5.28) R±{X)=UbTTb{R±W)- 

We cannot expect a non-trivial propagation theorem at points in R±{X) (since one 
of the broken bicharacteristics through them is stationary) , but we can expect such 
results elsewhere. We emphasize that there can be (and usually there are, though 
not if there are no bound states in the subsystems) non-constant bicharacteristics 
through R±{\); it is because there exists a constant bicharacteristic that we cannot 
expect an interesting propagation result. 

Note that if over an interval J C /, the image of 7 is disjoint from i?+ (A) U R- ( A) , 
then 7|j can be reparameterized using r as a parameter instead of t; Lipschitz 
functions on the reparameterized curve remain Lipschitz. 

Another interesting 7r-invariant function is r/ = rja = Za-i'a defined near ^'^T*{Ca] X) 
for some p E C'^. Since ^'^H^rj over Sfc(A) is 

(5.29) '^H'gTta = —{2TaVa + 2Wab\^), 

where Hb = {fJ-a, i^ab) is the decomposition corresponding to that of Xb with respect 
to Xa and its orthocomplement, so 

(5.30) iMbl'^lAiaP + li^afcP, 

we can obtain a monotonicity result for F = r]a ° J, defined for t near to. First 
note that F(to) = (as z == at 7(^0) € '"'T^, {Ca;X)). Next, F(t) > for < > to. 
Indeed, if F{ti) < 0, let tp = sup{t : F{t) = 0, t < h}, so to < tf, < ti; since F 
is continuous, F(t) < for t e (t[,,ti]. But by (HH), D±F{t) > cF{t) for some 
c > whenever F < 0, hence on [to,ti], so D±{e^'^-F) > there, so e~'^F is non- 
decreasing, so F{ti) > F(to) = 0, a contradiction. Thus, F{t) > for t > to, so 
by (|^), D±F{t) > -c'F{t) for some c' > 0, so D±{e''-F){t) > 0, hence e'='*F(t) 
is non-decreasing. Thus, with ria{t) — ?7a(7(0) — -^(0: there exists S > such 
that either rja is identically on [to, to + 6), or ria{t) > for t e (to, to + 5). Since 
^'^Hg\za\'^ = 4:{xa/x)ria (as rja is TT-invariant), we deduce that there exists S > 
such that either Za — for all t in (to, to + S), or Za ^ for all t in (to, to + S). If 
the former holds, then "f\(^to,to+s) is a differentiable curve in ^'^T^, Xa with tangent 
vector given by ^"^H^, by Remark pT3. In other words, we have proved the following 
lemma. 

Lemma 5.8. Suppose j is a broken bicharacteristic. Then there exists 6 > such 
that either j{t) e "''T^, {Ca;X) for all t £ {to, to + S), or j{t) ^ "''T^, {Ca;X) 
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for all t e {to^to +6). If the former holds, then ^lu^j-^^gj is an integral curve of 
^'^H'^. A similar result holds for {to — (5, to)- Ii^ the former case we say that 7 is 
forward/backward tangential at to, in the latter that 7 is forward/backward normal 
at ip- 

Remark 5.9. Note that tangential/normal behavior is the property of 7, not just 
of its value at to, unlike in the situation when there are no bound states in any 
subsystems. Indeed, in that case we only need to consider ^'^Hgf{£_o) for ^0 G 7ro(fo) 



in Definition 2.1, so in particular ^"^Hgrjai^o) = Wa\'^ — A — r^ — |/iap, so 7 is 



'9' 

^2 , I,, 12 



tangential if and only if r„ + |/ia| = '^ at 7(^0) G ^'^Tq, Xa- This remark also shows 
that our definition of broken bicharacteristics agrees with that of |34| , which in turn 
corresponds to the analogous definition introduced by Lebeau in the study of the 
wave equation p9| . In particular, our general definition is equivalent to the one 
discussed in the introduction if Ai = {0}. 

We can give a more geometric description of the broken bicharacteristics, pro- 
vided that the set of thresholds, Ai, is discrete, or if i? is a four-body Hamiltonian, 
A ^ Ai. Since this is improvement is not important for the propagation estimates, 
we only state the result below, and give the proof in the Appendix. Note that if 
there are no bound states in any subsystems, then Ai = {0} is certainly discrete, so 
in particular the proof applies in the setting of [ p4[ (though it was already proved 
there, with the proof based on Lebeau's results). 

Theorem 5.10. Suppose that Ai is discrete and 7 : M — > S(A) is a continuous 
curve. Then ^ is a generalized broken bicharacteristic of H — X if and only if there 
exist to < ti < t2 < ■ ■ . < tk such that -f\it^,tj+i], as well as 7l(-oo,to] "•'^^ 7[tk,+oo); 
are the projections of integral curves of the Hamilton vector field '^'^H'i for some a. 
Similar results hold if the interval of definition, R, is replaced by any interval. 

We end this section by analyzing the behavior of generalized broken bicharac- 
teristics under the assumption that A ^ Ai; this will be useful when studying the 
resolvent. 

Lemma 5.11. Suppose H is a many-body Hamiltonian, A ^ Ai, and 7 : / — > S(A) 
is a generalized broken bicharacteristic where I is a closed unbounded interval. Then 



(5.31) 7(/)C7(^)Ui?+(A)Ui?_(A). 

Moreover, if I is not bounded above, r(+oo) = linif^+00 ''"(7(i)) exists. A— t(-|-oo)^ € 
Ai, and there exists ^ e i?+(A) U i?_(A) and a sequence tj —f +00 such that 
^ — limj^oo7(tj) in S(A) (hence in particular t(^) = r(+oo)j. Similar results 
hold as t — » —00 if I is not bounded below. 

Proof. For the sake of dcfiniteness we take / = [to, +00), all other cases are very 
similar. 

Since t o 7 is monotone decreasing, and r is a bounded function on S(A), 
t(-I-oo) = limt^+00 T{t) exists (we wrote T{t) = T{j{t))). 

Now, let tn — > +00 be any sequence such that t„ > to for all n. Let 7„ : [0, 1] -^ 



S(A) given by 7„(t) = 7(i + t„). By Proposition 5.5, {7n} has a subsequence, which 
we write as 7^ — "fk{n) 1 which converges uniformly to a generalized broken bichar- 
acteristic 7' : [0,1] -^ T,{X). Thus, for t G [0,1], T{-f'{t)) = lim„^oo t(7;(<)) = 
lim„^oo j{t + tk(n)) = ^(+00), so T is constant on 7'. Since r o 7' is Lipschitz, it 
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is equal to the integral of its a.e. defined derivative, — 2|/i(t)p, so we conclude that 
|/z(i)p = a.e. along 7'. But then the arclength of the projection of 7' to §"~^, 
which is the integral of 1/^(^1, is also zero, so 7' is a constant curve, and hence, as 
shown above, there exists ^ G ^+(A) U i?_(A) such that 7'(i) = ^ for all t G [0, 1]. 
Thus, 



(5.32) 7([io,+oo))n(i?+(A)Ui?_(A))^0. 

Also, t{+oo) = t(^), and ^ e i?+(A) U i?-(A), so A - T{+oof e Ai. 

We still need to show that if {t„} is a sequence in [to, +00), and {j{tn)} converges 
in S(A), say ^0 = hm„^oo l{tn), then either ^ £ l{[to, +00)), or ^ G i?+(A)Ui?_(A). 
Now, if tn has a subsequence converging to some T G [to, +00), then 7(i„) converges 
to 7(r), so we are done. We may thus assume that i„ —> +00 as n ^ 00. By the 
argument of the previous paragraph, there is a subsequence ifc(„) such that 7^ = 
7(i + tfe(„)) converges uniformly (over [0, 1]) to a generalized broken bicharacteristic 
7', which is a constant curve, 7'(i) = ^ for i G [0,1], in R+{X) U i?_(A). In 
particular, 7(tfe(„)) converges to ^. But on the other hand, ^0 = lim„^oo 7(in), so 
Co = C e R+W U R-{X). This proves that 



(5.33) 7([io, +00)) c -/{[to, +(^)) U i?+(A) U i?_(A). 



D 



Remark 5.12. Note that under the assumptions of Theorem 5.10|, 7|[t,+ 



integral curve of '^'^H^ for a sufficiently large T, so 1^ = limt_>-|.oo l{t) exists in S(A). 
Moreover, if ^ G i?-|-(A) (i.e. it is incoming), then 7|[t,+oo) is a constant curve, so 

eG7([to,+(^)). 

6. Positive commutators 

Our approach to prove propagation of singularities along generalized broken 
bicharacteristics relies on positive commutator estimates. Just as in Q and [Q, 
we need to commute operators arising by quantization of 7r-invariant functions q G 
^ooj^sc J-* ^-j^ A = §!Ji as usual. That q is 7r-invariant means that at C'^, the regular 
part of Ca, given by Xa = 0, Zq = 0, qijja, Za,^a,O[za^=0 is independent of ^°, i.e. 
liVa, 0, £,a, C") is independent of ^a- This means that (except in some special cases) q 
is not symbolic at cotangent infinity, i.e. as ^ — > 00, since typically ^q derivatives will 
not give any decay in the ^° direction. Recall here that a (non-polyhomogeneous) 
scattering symbol r of, say, multiorder (0,0), r G .A(°'°^(S" x S"), is an element 
of C°°(]RJ^ X K?) that satisfies (product-type) symbol estimates in both sets of 
variables, i.e. for all a, /3 G N" there exists Cap > such that 

(6.1) [DZDlriw,O[<Ca0{w)-\^\{O-^^\. 

However, in our arguments the failure of q to be in yl(°'") (§" x §" ) will never cause 
a problem, only an inconvenience, since we always localize in the spectrum of the 
elliptic operator H by composing our operators with ipoiH) G ^g,,""' (^i^) where 
V'o G C^(I^)- Such operators tpo{H) are smoothing - their amplitudes decay rapidly 
in ^. 

So fix ijjo G C^(R; [0, 1]) which is identically 1 in a neighborhood of a fixed A. 
Thus, TpoiH) G ^S(°°' (A,C), so it is smoothing. At the symbol level, ipoiH) is 
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locally the right quantization of some 

(6.2) peC°°{[§l;C]x§l) 

which vanishes to infinite order at [§";C] x 9S", i.e. it is Schwartz in ^. 

We are then interested in the following class of symbols q. We assume that 
q € C°°(M.'^ X Wt) and that for every multiindex a, /3 e N" there exist constants 
Ca^jS and rria^p such that 

(6.3) \{DZDlq){w,0\ < a,^(u')-l"l(e)'"-'^- 
This implies, in particular, that 

(6.4) q G A^i^l X M"), 

i.e. that g is a 0th order symbol in w, though it may blow up polynomially in ^. 
Indeed, in the compactified notation, ( |6.3|) becomes that for every P £ Diffb(§"), 
acting in the base (w) variables, and for every /5 e N" there exist Cp^p and mp^p 
such that 

(6.5) |(P7^f)g|<Cp,^(0"-'^ 

It is convenient to require in addition that q be poly homogeneous on §" x M": 

(6.6) geC°°(S!; xR"); 

this stronger statement (implying ( |6.4[ )) automatically holds for the 7r-invariant 
symbols we are interested in. 

We next introduce the product symbol 

(6.7) a{w, w', = q{w, £,)p{w',£.)^ 

where ipo(H) is given locally by the right quantization of p. The main point is 



Lemma 6.1. The symbol a defined by (|6J) is in C°°(§'^ x [S'^;C] x W^) and it 
vanishes with all derivatives at §" x [§" ; C] x 9S" . Hence, it defines an operator 
A G ^5^' (§",C) by the oscillatory integral 

(6.8) Au{w) = (27r)-" / e'^'"-'"'^<a{w, w',^ uiw') dw' d^. 

If in addition q is tt -invariant, then for all S, G '^'^T* {Ca', X) Aa{S,) = q{S,)''Po{H)^{^). 
In particular, [A,H]^q vanishes for alia, so [A, H] G ^5^' (§",C). 

The last statements can be seen directly, but they follow particularly easily from 
the forthcoming discussion in which we quantize q itself (without the factor p) to act 
on oscillatory functions, and from the oscillatory testing definition of the indicial 
operators that we have adopted. 

Our immediate goal is now to prove that under certain conditions, the most 
important of which is that in a part of phase space, say U C S(A), ^'^H^q is negative 
for all 6, 

(6.9) i^{H)[A* A, H]iP{H) > B* B + E + F 

where B,E £ ^g^°^'° (§!;:, C) (so they are bounded on Ll^{^\)), F G '^^^^'^{^1,0) 
(so F is compact on i^c(S+))> ^^'sc{E) disjoint from C/ and V € C^(R; [0, 1]) has 
sufficiently small support near A. In view of B* B, this says that the commutator 
i[A*A,II\ is microlocally positive at U. The precise version of this statement is 



Proposition 6.5 at the end of this section. 
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Due to a square root construction in the ps.d.o. algebra, see |^, Corollary 9.7], 
this result follows immediately if we prove that all indicial operators oiiilj{H)[A* A, H]'i/j{H) 
over U are positive. There are two closely related ways to proceed prove the pos- 
itivity of the indicial operators. A general approach is to follow the proof of the 
Mourre estimate due to Froese and Herbst , but replace their global statements 
by the appropriate indicial operator ones. This approach is necessary in more gen- 
eral geometric settings. However, in the Euclidean setting it turns out that all one 
needs is the Mourre estimate in all the proper subsystems. We follow this path. 
Thus, we will need to evaluate some commutators quite explicitly, and use the 
Mourre estimate in the normal variables and the standard Hamilton vector field 
commutator estimate in the tangential variables. 

Although the evaluation of the commutator appears (and is) rather obvious if 
we write A = Qipo{H), think of Q as a scattering ps.d.o. obtained by quantizing 
the symbol q, and use 

(6.10) [A, H] = [QMH), H] = [Q, H]MH), 

this argument is rather formal since q is not a scattering symbol. Here we follow 
|37[ to make precise sense of this. 



So let q be as in Lemma 6T, and suppose that / e C°^{X;R), u — e^^l^v^ 



V € C°°{[X;C]). Thus, u e n^o^''''(IR") for s < ~n/2. Using our indicial operator 



-oo.O 



definition, A £ $g^ ' {X,C), the above calculation, ( 6.1C ), makes sense and gives 



the correct result if Q, defined by 

(6.11) Qu^{2tt)-^ I e'^<q{w,OHOd^, 



acts on such oscillatory functions. But u e n^oiJ'^'"(R"), so u € n^o-H"'''''(]R"), 
so the integral makes sense as a distributional pairing. The same holds if we re- 
place q{w,^) by gC^' = g(u',^)(^)~^'', and in addition, by choosing k large enough, 
g*^'') satisfies any fixed number of scattering symbol estimates as a scattering sym- 
bol of multiorder, say, (0,0). In addition, Qu = Q('=)(A + l)''u, and (A + 1)*^ : 
eif/^C°°{[X;C]) -^ e^//^C°°([X;C]). Since any fixed symbol estimate of e'^^'^'Pu, 
u £ e'-''/^C°°([X;C]), requires only a finite number of scattering symbol estimates 
on P, we conclude that 

(6.12) Q : e'-^/="C°°([X;C]) ^ e*//=^C°°([X;C]). 

In addition, to calculate Ru, R £ **c°(^), modulo a;'e*^/^C°°([X;C]) requires 
only a finite number of scattering symbol estimates on the amplitude r of R, we 
conclude that in all indicial operator calculations, as well as calculations of lower 
order terms, we can work as if Q were in 'I'*|;°(X). We also note that if i? G 'I'*^°(v?) 
has sc-principal symbol r{ya, Za,S,a,£,°'), its a-indicial operator, Rafiiya,S,a) is the 
translation-invariant ps.d.o. /(r) £ ^lf{X'^) on X"' given by 

(6.13) liriya,0,^a,-))v ^ Txlriya,0,U,n^x'^v, « G 5'(X'^); 

see [BTI |3. Here JFjfa denotes the Fourier transform on X"'. Correspondingly, when 
applied to w € e^^^^v, v £ C°°{[X;C]), e^^^/^Qu is given by l{q{ya,0,^a, ■))v mod- 
ulo xC°°{[X;C]), at C^. But q is 7r-invariant, so q{ya,0,^a, ■) is independent of ^", 
so I{q{ya, 0, ^a, ■)) is just multiplication by q.n{ya, 0, ^a), which we will simply write 
as q{tla,^,S,a) as in Section |[ Simlarly, ^^Hgq\z=Q is also independent of ^°, since 
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(zU=o is and ^^H'^g = 2\wa\-^^a-d^^ is tangent to -T^^S';, so I{{-'H;;q){ya, 0,^a, •)) 
is just multiplication by ("'=i7^q)^(yQ, 0,Cq, ■) = i^''H^q){ya,0,S.a)- 

While the preceeding discussion was quite general, we will be interested in rather 
special choices of q, which will make the indicial operator calculation simpler. In 
particular, these choices will allow us to use the standard Mourre estimate on X°', 
rather than essentially repeating its proof, as presented by Froese and Herbst g, 
see also |^, to obtain a corresponding estimate for more general operators in the 
normal variables. Thus, we will make a first order assumption on q at the collision 
planes at infinity, i.e. at the Cq, as well, in addition to its 7r-invariance. Namely, 
we assume that for all clusters a, q has the structure 

(6.14) dzq\z=o — Qi£.°' ■ dz, qi tt — invariant, 
so for all a, 

(6.15) g = go + qiVa + y^ ZjZjqij, 

with go,gi e C^{^''T*Xa), q^J € C°°{''T*§1). 

We use this as follows. For u € e^^^''C°°{[X;C]) we can write 

(6.16) 

[A, H]u = AHu - HAu = Q%Po{H)Hu - HQtljo{H)u = {QH ~ HQ)iJo{H)u. 

As indicated above, during calculation the indicial operator of [A, H], we can thus 
work as if Q were a scattering ps.d.o. Recall that our many-body Hamiltonians are 
of the form 

(6.17) i/ = i/« + Ax„ + la, la e xBiSliU), 



U a n eighborhood of C^, disjoint from Ct such that Cb 7^ Cq. Using ( 3.17 ) and 



(xl£), we write 



(6.18) [A,H]u=[Q,AxJMH)u+[Q,Ia]MH)u+[Q,H'']MH)u. 

Since the principal symbol in the scattering calculus is given by the Poisson bracket, 

for u = e'f/'^v, ij € C°°([X;C]), p G ff ([X; Ca]) with p[X; Ca]ip) e C^, 

(6.19) 

zx-'e-^f^-[Q,AxMp) = I{{-'''H^qmMp) = -(^=i?,"go)(0«b), df{p) = C, 
here we used that as ^'^H^q — ^"^H^qo is 7r-invariant, I{^'^Hgq) is a multiplica- 
tion operator. In addition, as la G xDiSg^{X,C), we see that a priori QIa,IaQ '■ 
e^^^^C°°{[X;C]) -^ xe^f^^C°°{[X;C]) with indicial operator given by the composi- 
tion of the two indicial operators. Since that of Q is I{q) = g Id, a multiple of the 
identity, we conclude that the (a, l)-indicial operator of [Q, la] vanishes, so in fact 
[QJa] ■■ e'f/''C°°{[X;C]) -^ x^e'f/''C°°{[X;C]). 

The only remaining term is [Q, H°-]'iIjo{H) . Note that Zj = \wa\^^{'w°-)j as usual, 
so ija — Iwal^^w"" • ^°. Thus, with Qi, Qij denoting the quantization of qi and qij 
respectively, the same way as Q is the quantization of q, 






(6.20) g = Oo + Iw^l-'Qiw^ ■D^a+^ z,ZjQ,j. 

i,3 

Here Qq, Qi depend on the ^^T*Xa variables only, so they commute with iJ". On 
the other hand, [zi, H""] e \E'gj. (§" ,C) (since all indicial operators oi Zi are multiples 
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of the indentity), and Q,j : e'f^='C°°{[X;C]) -^ e*^/"'C°°([X;C]), so the same holds 
for [Q,j,H'']. Thus, 

(6.21) [z,z,Qy,HT : e'f/^C°°{[X;C]) ^ |zpe'//-C°°([X;C]). 

This proves that for osciUatory functions u, 

(6.22) 

[Q,H'']MH)u^\wa\-'qi[w''-D^.,H'']MH)u + u', m' e |z| V-^'/^C°°([X;C]), 



where qi is evaluated at df as in (6.19). Thus, we conclude that 

where 
(6.24) 

Biya,U = MH"" + \U^)qqi Y,i^")jDi.o^)MH'' + l^aH e ^7r'~\X-) 

j 

is a ((j/a,Ca)-dependent) multiple of the generator of dilations on X'^. Although we 
worked at C^ for the sake of convenience, i[A*A, H]^ ^ is a well-defined element of 
C°°{'"'T*{Ca;X);^^^'"{X'',C'')y, in particular it is continuous on '"'T*{Ca;X), in 
which ^''T^,^{Ca; X) is dense, so ( |6^ ) in fact gives i[l*2^]^^ on ^''T*{Ca; X). 

We also°remark that if V G C;?°(M; [0, 1]), -0o = 1 on supp^, then tlj{H) = 
'iJjoiH)ip{H), so Ai;{H) = Qij}Q{H)i:{H) = Q^p{H) on oscillatory functions. Thus, 

(6 25) ^i^iH)[A^H]i,iH)^,{ya,U = - 2q^^H-qoiya,0,Ui^{H'' + U^f 

+ *[Sv,(2/a,Ca)+5v(2/a,Ca)*,i^1 

where 
(6.26) 

i 
This also proves the following lemma. 



Lemma 6.2. (see also [p4l Lemma 9.4]j Let q and A be as in Lemma 6.1, q be- 
ing TT-invariant. Suppose also that q is of the form (6.15). For every seminorm 
in ^5^' {X°',C°') and for every I G N there exist C > and m (^ N such that 
for every a and every ^ £ '^'^T*{Ca', X), p G C'a, the seminorm of [A, H]^ ^{^) in 
^-oo,0(.^a^(^a) -g i)ounded by 

(6.27) C{\qm + \q,{0\). 



Remark 6.3. Similar conclusions hold for every seminorm in C°°{'"=T*{Ca;X), ^'^'"{X^X")), 
which can be seen directly from our calculations in the preceeding proof. 



hi view of ( 6.24 ), the Mourre estimate is immediately applicable; we first recall 
it. We present it essentially as stated in [||, statement H^^a), Section 6.4. Let 

d°(cr) = infjcr - cr' : a' <a, a' e A' } > 0, if cr > inf A' 

^ ■ ' d°((T) = 0, if (T< inf a;,. 
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(6.29) d''^'^{a)=mi{d''{a'): a' e [cr - k,(t + k]}. 

Then given e > 0, CTq € K, there exists 5 > such that for aU cr < ctq and for aU 

ipi £ C^{R; [0, 1]) supported in {a - S,a + S), 



(6.30) 



Note that this definition of (i"(cr) differs sHghtly fro m tha t of g, namely our d°-{a) 
is bounded from above by theirs if cr > inf Aa, and ( 6.30 ) is trivial if cr < inf Aq. 

We apply this with a = A— |CaP, fo = A + 1. Thus, given e > there exists 
(5 > such that for all ip e C^{R; [0, 1]) supported in {\- S,X + S), 

> ^id'^'^X - ICaP) - e^H^ + \UY- 
Note that if cr > inf A'^, 
(6.32) 

2d'^(a) = 2inf{|4"|2 : a - l^' e speCpp(il^), Ct D C,} 

= inf{^^iJg^,(y,,0,Ca,4") : ^ - 141' e speCpp(ij''), ^ D Ca} > 0; 

here ^'^Hgija = |CaP is independent of j/a and ^a, so the appearance of these on the 
last line is irrelevant. For the sake of convenience, in the next paragarphs, we write 
inf = 0, then ( |6.32 ) holds for all cr e R. We remark that ^" • w" is the principal 
symbol of {w°' -D^a +D^a • w")/2 in the scattering calculus (here the part at infinity 
is the one that matters), which explains how it arises in the commutator estimates. 
Thus, with/ = [A-(5,A + (5], 

(6.33) 

2d''''(A-|Can 

^inirH'g^aii) 
= inf{-i/g^40 
= inirH'rjaiO 



A' - 141' - 141' e specpp(i/''), Cb D Ca, A' e [A - 5, A + S]} 
TTbiO =4 CbD Ca, 3<J e spec„„(H^), |A - leT - fil < <5} 



In our estimates, in the part of phase space where we wish to prove that the 
commutator is positive, we will have 

(6.34) q>0, qi<0. 

Note that due to ( |6.14D , for b with CbD Ca, 



(6.35) 



'H'q\. 



H^qo + qi"'m7]a 



sc Tja 
^9 



The Mourre estimate thus shows that 

(6.36) 

i^PiH^ + l^a^lB + B*,H^]i:{H^ + \ia\^) > -iqq^id^^X - \^a\') - e)7P{H- + \(a\y 



Substituting this into ( |6.23 ) yields 



(6.37) 



t^iH^ + \unA*A,H]Jo^{H'^ + ic.n 



> -2qrH^qo + 2q^d-^\X ~ \U') " 2egi)^(ff'' + |4 



|2\2 
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with 



where we used ( 3.33 ), ( 6.35| ). In ( |6.38 ), inf must be understood as —^'^H'^qQ, but 



due to the factors of il){H°' + |^q| ) in ( 6.23 ), that equation holds for any value 



replacing ^'=iJ^go + 2gid''''5(A-|CaP)-2egi if A-|CaP < inf A;,-(5 (which is exactly 
when inf arises above) . 

In order to estimate the commutator on some set U C *''^T*§", it is convenient 
to rewrite this as follows. Suppose that b G C°° {^'^T*'§,'^) is 7r-invariant, 6q > 0, 
/ = [A — (5o , A + (5o] , for all clusters 6, 

(6.39) i e n-jiU) ^ -F^^d) < -6(eT, 
and there exists C > such that for all ^ £ U 

(6.40) |gi(0| < CbiO'. 

Then for all e' > there exists S' > such that if supp ip C {X — S' , X + 6') then 

(6.41) ^^l;{H'^ + W)[A^Hl,i;{H'^ + W) > (2 - e')qhMH- + \U)\ 

We have thus proved the following proposition. 

Proposition 6.4. (cf. [p3. Proposition 9.6] j Suppose that H is a many-body Hamil- 
tonian and A G M. Suppose also that q-,b € C°°('"^T*§!j;;IR) are iT-invariant, for all 
clusters a 

(6.42) q = qo + qiVa + 0(|zq|^), go,9i TT - invariant, 

they satisfy the bounds ( |6.5| ), q,b,—qi > 0, and that there exist S > 0, C > 0, 
C" > 0, U <Z T,{I) , where I ~ [X — 5,X + S\, such that for every index a and for all 

(6.43) i e 7T-}iu) => ^'H^qio < -m^ 

and 

(6.44) C G vl(C^) ^ 9(0 < C^KO' «^d 191(01 < C'biO'. 

Let A G ^5c^' (-^iC) 6e as m Lemma \6.^ . For any e' > there exists 5' > such 
that if t/j e C^O^) IS supported m (A -S',X + S') and ^ e '"'T*{Ca;X) for some a 
with 7r(^) G U , then 

(6.45) i{rl:{H)[A^H]i^{H))^,{0 > (2 - e')6'gV(^a,o(0)'- 

We can now apply |34| Corollary 9.7] to obtain estimates on the original oper- 
ators. Note that the proof of this corollary consists of a square root construction, 
|34| Proposition 8.3], for which the existence of L^ eigenvalues in subsystems is 
irrelevant, and an indicial operator estimate, which was given in Proposition 9.6 in 



l34|; here its place is taken by the preceeding result. Proposition 6.4. Recall that the 
square root construction is via a ps.d.o. version of the functional calculus, and that 
(3.4(;) ensures that this can be done within the ps.d.o. calculus. Also, (6.46|) might 



appear weaker than its analogue in [p4|, but due to (6.42), the other estimates of 
Corollary 9.7 of that paper follow automatically. Note that 7r-invariant (continu- 
ous) functions on >5cy*gn ^^^^ |-,e regarded as (continuous) functions on ^'^T*§", so 
it makes sense to talk about the support of these functions in ^'^T*S\. 
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Proposition 6.5. (cf. p4, Corollary 9.7]j Suppose that the assumptions of Propo- 
sition 6.4 are satisfied and let C be as in (pA4), and U open. Suppose in addition 



that for any cluster a and any differential operator Q on '"^T*{Ca\X) there exists 
a constant Cq such that 

^(0 e U and 6(0 ^ 
f6 46) 

^ \Q{b-^qm\ < Cq and \Q{b-'qim\ < Cq. 

Then for any e' > 0, M > 0, and for any K C '"^T*X compact with K <Z U there 
exists (5' > 0, tA G C^(R) is supported in (A — (5', A + 5'), ■(/; = 1 near \, B,E Cz 
*sr'°(X,C), F e ^-sT''\x,C) with ^¥'sAB),^¥'sAE).^Y's,{F) c t{I), 

(6.47) 

"NY's.iE) nK^9, WF's,{F) C suppg, Ba,o{0 = KOqiO'^^iKHaAO), HO e K, 

.such that 

(6.48) 

i^P{H)x'^/^[A* A, H]x-^/^i^{H) - Mil){H)A* Ail){H) > (2 - e' - MC)B*B + E + F. 



Remark 6.6. Multiplying both sides of ( 6.4S| ) by -01 (^) such that V' = 1 on supp^'i 



shows that "0 can be replaced by any such 0i. 

7. Propagation of singularities 

For our positive commutator estimates it is convenient to consider two scenarios 
separately, though we present them parallel to each other. Recall that we are 
interested in proving propagation estimates in 5 = S;(A) \ (i?+(A) U i?_(A)). Note 
that f e 5 n '''^T^, Xa, ^ = {ya,fa,fia), mcaus that for aU b with Ca C Cb, 6 G 
^b^(0: l^bi^b) 7^ 0. The two possibilities are: 

(i) p,a 7^ 0. Then ^"^Hg does not vanish at ^, so there is tangential propagation 

along Ca. 
(ii) /2o = 0. Then there is no tangential propagation along Ca. (However, note 
that as ^ G S this automatically gives ^ ^ Sa(A), so certainly the particles 



cannot be in an a-bound state.) In this case, by (5.16) 



(7.1) "'K(Za ■ Va) = {Xalx){2Ta{Za ■ Va) + '^WhV) > 



(near tt^ {^)) since Hb ^ there, so Vab ^ by (5.30). Thus, we can use 
Za • t'a as a parameter along broken bicharacteristics. 
The main technical result on the propagation of singularities is the following. 

Proposition 7.1. Suppose that H is a many-body Hamiltonian. Let u G C^°^{X), 
A ^ Ai. LeU e 5 n "'T^,^Xa, C = iya,fa,pa), and suppose that ^ i V^Fsc{{H - 
X)u). Let rja = Za- Va- Then there exists a neighborhood Uq of ^ with the following 
properties. 

(i) If p-a 7^ and there exists a neighborhood U C Uq of ^ in S(A) and if there 
exists So > such that exp(— so'"^i?^)(0 ^ WF5c(^*) and for all s satisfying 
< s < sq, exp(— s^^i?^)(0 e U (i.e. the integral curve segment of'^'^H^ is 
completely inside U ), and in addition 

(7.2) eeC/ anrfr?a(O<0^e^WF5c(u), 

theni(^^¥sc{u). 
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(ii) If p,a — and there exists a neighborhood U d Uq of ^ in S(A) such that 
(7.3) ^eU and7jaiO<0^H'^Fsciu), 

then^^WFsc{u). 

Remark 7.2. Note that r]a{0 < implies Za ^ 0, so ^ ^ ''"T^, Xa- 

Proof. The main step in the proof is the construction of an operator which has 
a microlocally positive commutator with H near ^o- In fact, we construct the 
symbol of this operator. This symbol will not be a scattering symbol, i.e. it will 
not be in C°°(S" x §"), only due to its behavior as ^ — > oo corresponding to its 
TT-invariance. This will be accommodated by composing its quantization with a 
cutoff in the spectrum of H, ip{H), ip S C^iM) supported near A, as discussed in 



Lemma 3.1. This approach simply extends the one taken in |37|-[pJ|, though the 
actual construction is different due to the more complicated geometry. 
We introduce some notation and then fix Uq. We define cq > by 

(7.4) CO - i inf{--i7^V(6) : 6 e ^H^), ^ D Ca} > 0. 
Note that at Ca^ where |u;q|/|i(;| — x/xa — 1, 

(7.5) T = Ta- ria, SlUCC T = — - , Tq = j — , r]a = "i T' 

\w\ \Wa\ \Wa\ 

and if pa = 0, then ^'^HgTa{£,b) — when 7r(^(,) — ^. Thus, if pa — 0, then 

(7.6) CO = i mirn'^iza ■ i^aM) ■■ a e n^'iO} > 0. 

Now fix / = [A — (5o, A + (5o] and Uq to a neighborhood of ^ in S(/) such that 
Uq c UI^^TB'^b ■ ^O' ^ C-fc}, Uo n WFsc((ii^ - X)u) = (this is possible since 
WFsc((-ff — A)u) is closed), Uq is inside a fixed system of local coordinates, and 

if /2, ^ then e e C/o, C e KJiO, CbDCa^ -''H'^ri^t) > 3co/2 

iifia^Q then ^ e Uo, e& e KJiO, CbZ>Ca^ ""H^iza ■ Va){ib) > 3co/2. 

If jia ^ 0, we make an additional definition. Namely, we let (5 > be given by 
(7.8) -5 = f-r(exp(-so^^i?^")(0). 

We will use e £ (0, 1) as a small parameter that microlocalizes in a neighborhood 
of ^ ii p.a = 0, and in a neighborhood of {exp(— s*''^i?g)(f) : s 6 [0, sq]} if p-a 7^ 0; 
ii pa ^ we will always take e < S. 

Employing an iterative argument as usual, we may assume that ^ ^ WFg^ (u) 
and we need to show that f ^ WFg^ (u) . (We can start the induction with an 
I such that u e H*^\X).) 

Our positive commutator estimates at a point ^ = (jja,f,pa) will arise by 
considering functions 

(f> = f-T + e-Wza\^+Lj), if Ma 7^0, 

(7-9) /3 

0= ^^a • i^a + -(l^aP +w), if Ma = 0, 

e 
where cj localizes in the tangential variables i/a, Ta, fj-a. In fact, if Ma 7^ 0, then 
^'^Hg{^a) 7^ 0, so we construct w to be a 'quadratic distance' (in ^'^T^^Xa) from the 
^'^H"i integral curve through ^q constructed so that ^^H^uj ~ 0. That is, we define 
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w on a hypersurface through ^ that is transversal to ^'^H^, e.g. t^ = fa, to be given 
by a positive definite quadratic form in some local coordinates centered at £,, e.g. 
iVa — VaY + (Mo ~ Mo)^, and extend it to ^'^T^^Xa to be constant along the ^'^Hg 
integral curves; cf. |3J, Section 7]. We can then extend a; to a function on ^'^T*^X 
defined near Ca as discussed before ( 5.11 ), and then 

(7.10) """Hgoj = 0, so "'^-ffgcj = for aU h. 
On the other hand, if jia — 0, we can take 

(7.11) U;^ = {^a-Lf + {.ya-Va?\ 

now 

(7.12) -ij;c.o = 2(y, - va) ■ ""^H^y, = 4^., • (y, - yj, 
so {^"^HgUJol — C'ujQ. Thus, for C sufficiently large, we see that for 

(7.13) uj = e-^^ujo, ''-ff> > 

since —^'^HgT > on S', hence near ^. Since uj depends on the tangential variables 
only, we conclude in either case that for all h, 

(7.14) ""^hIlu > 0. 

Moreover, ^'^Hg\za\'^ — ^Zab-^ab under the decomposition ^^ = (^a, i^ab), so ^'^H^g\za\'^ < 

Cl\Za\. 

Now suppose that /Iq ^ 0, 

(7.15) 0< 2e, f-T>-S-e. 

Since u; > Q and |zap > 0, the first of these inequalities implies that f — t < 2e, 
and the combination of these two gives </> — (f — r) < 6 + 3e < 4:6 a.s e < 6. Then 
we conclude that 

(7.16) -d-e<f-T<2e, \za\ < 2{eSy^^, w^/^ < 2(<5e)i/2. 
Note that for all b with Ct D Ca, Zat ■ i^ab ^ Za- Va at ^'^T^ Xf,. Thus, if 

(7.17) Za-Va> -Cseco 
we see that 

(7.18) ""Hl\za\^ = 4Za ■ Va > -iCseCQ. 

Thus, using (^J) as well, we deduce that there exists C4 > (independent of e > 0) 
such that 

(7.19) if za-i^a> -C4e, then ^^iJ^V > -'^H^gT - co/2 > cq. 

Note that here e G (0, S) is arbitrary. 
Similarly, suppose that p,a = 0, 

(7.20) (/) < 2e, Za-iya> -2e. 
Then we conclude that 

(7.21) \za ■ i^al < 2e, |z„| < 2e/y^, lu^^^ < 2e/^. 



Thus, from (|7.7D , 

(7.22) ^'^H'c^ > ^''H'iza ■ va) - 2y^Ci > co if /? = (co/4Ci)^ 
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The positive commutator estim.ate then arises by considering the following tt- 
invariant symbol q and quantizing it as in Lemma |6.l| . Let xo G C°° (M) be equal to 
on (-00,0] and xo(0 = exp(-l/t) for t > 0. Thus, Xo(0 = ^"^Xo(i), t > 0, and 
XoW = 0, t < 0. Let xi G C°°(K) be on (-oo,0], 1 on [l,oo), with xi > and 
Xi(i) = exp(— 1/i) on some small interval (0, to), to > 0. Furthermore, for Aq > 
large, to be determined, let 

(7.23) q = Xo{Ao\2 - 0/e))xi((T - r + S)/e + 1) 
if jla 7^ 0, and let 

(7.24) q = xoiAo\2 ~ ^/e))xi{iza ■ i^a)/e + 2) 
if fta = 0. Thus, q{(_) = xo(2/Ao) > 0, and on suppg we have 

(7.25) (j) < 2e and f — t > —6 — e in case (i), 
and 

(7.26) < 2e and Za ■ i^a> — 2e in case (ii), 

which are exactly ( [7.15 ) and ( 7.20 ), so suppg is a subset of ( 7.16 ) and ( 7.21 ) 



respectively. We also see that as e decreases, so does suppg — suppg^, in fact, 
if < e' < e then q^ > on suppg^'. Note that in case (i), by reducing e, we 
can make q supported in an arbitrary small neighborhood of a compact backward 
bicharacteristic segment through ^, and in case (ii), by reducing e, we can make q 
supported in an arbitrary small neighborhood of ^. 

We at once obtain positivity estimates for ^'^H^q. The following argument works 
similarly for both /1q = and jla ^ Q', we consider the slightly more complicated 
case p,a ^ 0. Thus, if /!„ 7^ 0, then 

^'H^q = -A-,^e-^x:,{A^\2 4>/^))Xi{{f - r + 5)/e + l)-ff^^0 

-e-\,{A-\2-c^/e))^,{{f-T + 6)/e+irHlT. 

We b reak up the first term by using a cutoff that ensures that the hypothesis in 
( 7.19| ) is satisfied. Thus, let 



(7.28) x = Xi(4^a-i^a/(C4e) + 2), 
so 

(7.29) on suppx, Za-Va> -C4e/2, 

(7.30) and on supp(l — x), Za- ^a < —Cie/A. 
Then 

(7.31) '^H^q = -bl + e, 
with 

(7.32) bl = A,h-'xoiAo\2 - <^/e))xi((f - r + S)/e + l)rH^^<t>. 
Hence, with 

(7.33) b^ = coA^'e-\'„{A^\2 - 0/e))xi((r - r + S)/e + l)x, 
using ( [7.19 ) and (7.29), we deduce that 



(7.34) """H^^q < -b^ + ee- 
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Moreover, 

(7.35) b^ > {coAo/16)q 

since (p > f — t > —2d on supp q, so 

3 Xo(A-'(2 - 0/e)) - Ali2 - q^/e)-\oiA^\2 - cj^/e)) 

>(A§/16)xo(V(2-<^/6)). 

On the other hand, Cc is supported where either 

(7.37) -5-e<f-T<-5, w^l^ < 2{eSy^^, \za\ < 2{eSy^^ 

so near the backward direction along ^"^H^ bicharacteristic through f, or 

(7.38) 



Za 



Va < -C4e/4, and - J - e < f - t < 2e, lu^/^ < 2(e(5)^/^ \za\ < 2(e(5)^/^ 



But by our assumption exp{—so^''Hg){^) ^ WFgc(M), Sq > 0, so the same holds 
for a sufficiently small neighborhood of exp{—so^'^Hg){^) as WFsc(u) is closed. By 
choos ing e > sufficiently small, we can thus make sure that the region defined by 



( 7.37 ) is disjoint from WFsc(u). Moreover, by further reducing e > if nece ssary 



and using our second assumption, we can also make sure that the region ( [7. 38]) 



is also disjoint from WFsc(m), so that suppcc is disjoint from WFsc(m) for all c. 



Moreover, by ( 7.16 ), for e > sufficiently small, we deduce from the inductive 



hypothesis that suppq (hence supp 6) is disjoint from WFg^(u) n S(/). 

Moreover, with d denoting a partial derivative with respect to one of (y, z, r, /u, i'), 

dq = -A^'e-'xoiAo'C^ - '^A))xi((t - r + S)/e + 1)50 

-e-'xo{A^\2-4>/e))x[{{f-T + S)/e+l)dT. 

At any Ct with p € Cb, defined by x = 0, z" = 0, as above, (p is independent of v" 
at z" = so outside supp Cc 

(7.40) \d^"dq\ < Cf^b^ at z" = 0. 

In fact, outside supp Cc, but in the set where b is positive, 

(7.41) b-^dq = CQ^dcf), 



so the uniform bounds of (3.46) also follow. In addition, at any cluster b, \za\'^ = 
|za6p + |ZqP, ?7a = ^afc ■ i^ab + Zq ' i^a ^ Zgb ■ i^ab + fib-, and Zafc, Vab are 6-tangential 
variables, so </), hence q has the form (6.42) around each Cf,. 



Let ip £ C^(]R) be identically 1 near and supported close to 0. We also define 
(7.42) q = ^(x)g. 



Thus, q G Qoo ^scrp* j^--^ jg ^ TT-invariant function satisfying (6.5). Let A be the 
operator given by Lemma |6.l| with q in place of q, so in particular its indicial 
operators are <?(C)'/'o(-ff&,o(0)- Note that ( |6.44| ) holds with C = Wc^^A^^. So 



suppose that M > and e' > 0. Choose Aq so large that MC < e', and let U be 
the complement of Uc7rc(supp Cc) in E(/), and let K = WFsc(w) r\T,{I); so K C U 
by our choice of e > 0. By Proposition p^ we deduce that there exists 5' > 0, 
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ijj e Cf (R) is supported in (A ~ S' , X + S'), tp = 1 near X, B,E e ^s;^'°(X,C), 
F e vl/--.i(X,C) with WF's,{B),WF's,iE),WF's,iF) C E(/), 

(7.43) 

WF's,{E) n if - 0, WF^,(F) C SUPP9, B^AO - ^(^^(O'^Vli^a.olO), ^(0 G ^, 
such that 
(7.44) 
i'il!{H)x-^/'^[A*A,H]x-^^^'ilj{H)~Mij{H)A*A'ilj{H) > {2 - 2e')B* B + E + F. 

Let 

(7.45) Ar = x-'-^/^{l + r/x)-\ r e (0, 1), 

so Ar e ^s;"'+^/^(X, C) for r e (0, 1) and it is uniformly bounded in -^I'^^^^^^^X, C). 
The last statement follows from {l + r/x}^^ being uniformly bounded as a 0th order 
symbol, i.e. from {xdx)''{l + r/x)^^ < Ck uniformly {Ck independent of r). We 
also define 

(7.46) Ar = AArX-^/'^ij{H), B,. = BAr, Er = A^^A^. 
Then, with V'o G C^(M; [0, 1]) identically 1 near suppi/', 

(7.47) 
ix'+^/^[AlAr,H]x'+^'^ 

= i{l + r/x)-^ilj{H)x-^/^[A*A, H]x-^/^iIj(H)(1 + r/x)-^ 

+ fiP{H)A*x^+^/'^[ArX-^/^ , H]{1 + r/x)-^x-^^^iPoiH)A^{H) 

+ iilj{H)A*4>Q{H)x-^/^{l + r/x)-^[ArX-^/^,H]x^+^/^A'4>{H) + Hr, 

where Hr is uniformly bounded in Vl/g(.°°' (X, C). Note that H^ arises by commuting 
A, powers of x and A^ through other operators, but as the indicial operators of A 
and X are a multiple of the identity, A^ x and A^ commute with these operators to 
top order, and in case of A^, the commutator is uniformly bounded as an operator 
of one lower order. Then, multiplying ( 7.44 ) by {l + r/x)~^ from the left and right 
and rearranging the terms we obtain the following estimate of bounded self-adjoint 
operators on Ll^{X): 

ix^+^''^[A*.Ar,H]x^+^l'^ - ij{H)A*{G; + Gr)Ai){H) - Mij{H)A*Ai;{H) 

^^■^^^ > a;'+i/2((2 - e')B:.Br + Er + Fr)x'+'/^ 

where 

(7.49) Gr = iMH)x^^^^il + r/x)-^[ArX'^/^,H]x'+^^^, 

and Fr G *s^°°'"^'+^(X,C) is uniformly bounded in •i'g^'^'^\x,C) as r -> 0. Now, 
Gr G ^Sc°°' (-^j C) is uniformly bounded in ^g^ {X, C), hence as a bounded operator 
on Ls(.(X). Thus, if M > is chosen sufficiently large, then G^ + G* > -M for aU 
r e (0,1), so 

(7.50) ij{H)A* {Gr + G*r+ M)At/j{H) > 0. 



Adding this to ( 7.48| ) shows that 

(7.51) zx'+'^'[A*rAr,H]x'+'^' > x'+'/^{{2 - e')B:Br + Er + Fr)x'+^'\ 
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The point of the commutator calculation is that in Ll^{X) 

{u,[A;Ar,H]u) 
(7.52) = {u, A;Ar{H - X)u) - {u, {H - X)A*.Aru) 

= 2i Im(u, a; Ar (iJ - A)u) ; 



the pairing makes sense for r > since Ar E '^g^' {X,C). Now apply (7.51) to 



-Sc 

^-;-i/2y and pair it with x^^^^l'^u in L\J^X). Then for r > 
(7.53) ||Srw|P < |(u,^rw>| + |(M,i^rW>| + 2|(u,A;A^(i7- A)m)|. 



Letting r — > now keeps the right hand side of (7.53) bounded. In fact, Ar{ii — 
X)u e C°°[X) remains bounded in C°°(X) as r ^ 0. Similarly, by ( [7.43[ ), E,.u 
remains bounded in C°°(X) as r -^ since WFsc(u) n WFg^(£') = Al so, F^ is 
bounded in S(iJ™''(X), iJ-™'-'(X)), so (m, i^^w) stays bounded by (fzl^ ) as well 



since suppq n WFgj, (w) n S(/) = 0. These estimates show that BrU is uniformly 
bounded in LlJ,X). Since (1 + r/x)-i ^ Id strongly on S(i7™''''(X), i/™''''(X)), 



we conclude that Bx ' ^/^m G Lg^(X). By ( 7.43 ) and Proposition 4^ this implies 
that for every to, 

(7.54) e^WF™''+i/2(«); 

in fact that t, i WF";'+^/^(w) for aU ^ e X for which g > 0. This is exactly the 
iterative step we wanted to prove. In the next step we decrease e > slightly to 
ensure that WFg(,(F) C suppg is disjoint from WF^' (u), just as Hormander 

decreases t in the proof of 



12|, Proposition 24.5.1]. D 

Based on this proposition, we proceed inductively, assuming that propagation 
has been proved at C^ with Ca C C;,. Noting that r]a = Za ■ Va < ^ implies that a 
backward broken bicharacteristic through ^ must stay away from Ca, we can use 
Lebeau's argument, as presented in ||34| , Theorem 10.6], to prove propagation along 
broken bicharacteristics. The basic idea is that S, G WFsc(u) H ^"^T^iXa implies 
that either every point on the non-constant backward '^'^_ff°-bicharacteristic is in 
WFsc('u) (in a neighborhood of ^), or there are points in WFsc(u) \ ^"^T^iXa ar- 
bitrarily close to ^ such that there is a backward broken bicharacteristic through 
these points which is completely in WFgc(M) (by the inductive hypothesis). Using 
the compactness of the set of broken bicharacteristics, we can extract a sequence 
of backward broken bicharacteristics which converge to a backward broken bichar- 
acteristic through ^, which will prove the propagation statement. 



Theorem. (Theorem 2.i) Letue C"°°(X), A ^ Ai. Then 



(7.55) WFsciu)\WFsciiH-X)u) 

is a union of maximally extended generalized broken bicharacteristics of H — X in 
±{X)\WFsc{(H-X)u). 

Proof. As usual, broken bicharacteristic means generalized broken bicharacteristic 
in this proof. 

We only need to prove that for every cluster a, if 

(7.56) C G WFsc(w) \ WFsc((i/ - X)u) and ^ e ''^T^, Xa 
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then 



(7.57) 



3 broken bicharacteristic 7 : [— eo,eo] — *■ S(A), eq > 0, such that 

7(0) = ^" 7(0 G WFsc(^) \ WFsc((i/ - X)u) for i e [-eo, £0]- 

In fact, if (|7^)^([r57|) holds for aU a with C^ C Ca, let 

(7.58) 

TZ ^{broken bicharacteristics 

7 : («,/3) ^ (WFsc(u) \ WFsc((if - A)u)) n [J{'''T*„^Xa ■ C, c CJ, 

a < < /?, 7(0) - a, 

and put the natural partial order on 7?,, so 7 < 7' if the domains satisfy (a, /?) C 
(a',/3') and 7 = 7'|(a,/3)- Then 7?. is not empty (due to (7.57)) and every non- 
empty totally ordered subset of TZ has an upper bound, so an application of Zorn's 
lemma gives a maximal broken bicharacteristic of 7? — A in the intersection of 
WFsc(u)\WFsc((7i7 — A)u) with Uc^cCa^'^T^' ^a which passes through^. A similar 
maximal statement holds if we replace Cc C Ca by Cc C. Ca- 
Indeed, it suffices to show that for any a, if 



(7.59) 
then 

(7.60) 



^ e WFsc(^i) \ WFsc((i? - X)u) and ^ G ''T^> [Ca] X) 



there exists a broken bicharacteristic 7 : [— eo, 0] 
7(0) = C" lit) G WFse(^) \ WFsc((i? - 



-^ i](A), eo > 0, 
A)u), te [-eo,0], 



for the existence of a broken bicharacteristic on [0, eg] can be demonstrated similarly 
by replacing the for ward propagation estimates by backward ones, and, directly 
from Definition 2.1, piecing together the two broken bicharacteristics gives one 
defined on [— eo,eo]. 



by Proposition 7.1 



We proceed to prove that ( 7.59 ) implies ( 7.60 ) by induction on a. This is certainly 
true for a 



there are no normal variables Za, I'a, so rja = 
in the notation of that Proposition, showing that a segment of the backward 
bicharacteristic through 1^ must be in WFsc(w). Of course, this is simply Melrose's 
propagation theorem p3[ Proposition 7]. 

In addition, if a is arbitrary and S, G -R+fA) U _R _(A) , then the constant curve 7 
through ^ is a broke n bich ar acteri stic, so ( 7. 59 )^( 7. 6C ) holds with this 7. 

So suppose that ( |7.59| )^( 7.60 ) has been proved for all b with Ca C Cf, and 
that f S S(A) n ^'^T^i {Ca',X) satisfies ( 7.591 ). As noted in the first paragraph, we 
thus know that the intersection of WFsc(w)\ WFsc((-ff- A)m) with Uc^cCt'^'Tc'Xb 
is a union of maximally extended broken bicharacteristics oi H — X. We use the 
notation of the proof of Proposition 7.1 below. Let Uq be a neighborhood of ^ = 
{ya,fa,pa) in J^(A) as in Proposition 7T; we may assume that Uq n WFsc((-ff — 
X)u) — %. By Proposition 7.1, either every point on the non-constant backward 



^i?"-bicharacteristic segment through ^ is in WFsc(u) (in the neighborhood C/o of 



^), in which case we have proved (7.60), so we are done, or there is a sequence 
of points ^„ e [/o such that ^„ e WFsc(w), ^„ ^ ^ as n — > 00, and ?7a(Cn) < 
for all n. Since ria(£,n) < 0, ^„ ^ ^'^T^iXa- Thus, by the inductive hypothesis 

there exist broken bicharacteristics 7„ : (— e„,0] ^ S(A), e„ > 0, with 7„(0) = ^„, 
7„(t) e WFsc(w) for all t G (— e„,0], and 7„ is maximal with this property in 
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UofMiCaCCb^'^T^' XbCiWYsciu) . That is, if 7„ is an extension of 7„, then 7„(— e„) ^ 
Uc„cCb'"'T'p,Xb, so 7„(-e„) € ''T,!., Xq. By Remark 57 we can fix some ep > and 



extend (or restrict) each 7„ to a broken bicharacteristic defined on [— eo,0], which 
we keep denoting by 7„; by the previous remark 7„(— £«) S ^'^TQ,Xa- But e'~'*ria is 
non-decreasing (for sufficiently large C > 0) along broken bicharacteristics by the 



argument that preceeds Lemma 5.8, so we conclude that 

(7.61) VahrM < e^(*-*'''77a(7n(0)) < 

fori e [-eo,0], hence ZQ(7„(t)) 7^ 0, so e„ > eo, 7„(t) S WFsc(u) for alH e [-eo,0]. 



By Proposition 5^, there is a subsequence of 7„ converging uniformly to a broken 
bicharacteristic 7 : [— eo,0] — > S(A). Since WFsc(u) is closed, 7 : [— eo,0] -^ 
WFsc(w)- In particular, 7(0) — ^ and 7(i) G WFsc(w) for aU t e [— cq, 0], providing 
the inductive step. D 

Remark 7.3. As expected, the theorem does not provide any interesting information 
at the radial sets i?+(A) U i?_(A). 

This theorem essentially says that the propagation of quantum particles, modulo 
'trivial' (i.e. Schwartz) terms can be understood as a series of collisions in which 
the total energy and the external momentum are preserved. 

8. The resolvent 



Before we can turn Theorem 2.2 into a result on the wave front relation of the 
S-matrix, we need to analyze the resolvent. More precisely, we need to understand 
the boundary values 

(8.1) R{X±iO) = {H-{X±iO))-^ 

of the resolvent at the real axis in a microlocal sense. To do so, we also need 
estimates at the radial sets i?±(A). Since some of the Hamilton vector fields ^'^Hg 
of the metric g vanish at i?+(A) U i?_(A), the estimates must utilize the weights 
a;~'~^ themselves. In this sense they are delicate, but on the other hand they 
only involve x and its sc-microlocal dual variable r, so they do not need to reflect 
the geometry of C. The best known positive commutator estimate is the Mourre 
estimate, originally proved by Perry, Sigal and Simon in Euclidean many-body 
scattering [ p7[ , in which one takes q = x~^t with the notation of Section |6|. Since 
it is easy to analyze the commutator of powers of x with H (in particular, they 
commute with V), the functional calculus allows one to obtain microlocal estimates 
from these, as was done by Gerard, Isozaki and Skibsted ||, 0]. Thus, nearly all the 
technical results in this section can be proved, for example, by using the Mourre 



estimate and Theorem 2.2. In particular, apart from the propagation statements, 
they are well-known in our Euclidean many-body setting. 
We first recall the limiting absorption principle. 

Theorem 8.1. Suppose that H is a many-body Hamiltonian, and A ^ A. Let 
f e C°°{X), uf = R{\ ± it)f, t > 0. Then uf has a limit u± = R{X ± iO)f in 
H^/{X), I < -1/2, ast^O. In addition, 

(8.2) WFsc(w±)c$±(i?^(A)). 
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Remark 8.2. The first statement in the theorem also holds if we merely assume 
/ e iJ™'''(X) with /' > 1/2, but then WFsc(w±) has to be replaced by the filtered 
wave front set WF^' ^ ("±)- Moreover, R{X±iO) give continuous operators from 

Proof. The existence of u± in 7J™''(X), I < —1/2 follows from the Mourre estimate, 
as presented in ||2^. That ^ € WFsc(m+) imphes t(0 < (in fact, r(^) < -^a(A) 
where a{\) = inf{A — cr : cr S A, a < A}), follows from the work of Gerard, 
Isozaki and Skibsted ; see pil to see how the proof would proceed with our nota- 



tion. By Theorem p.2|, WFsc(m+) \ WFsc(/) = WFsc(u+) is a union of maximally 
extended generahzed broken bicharacteristics. So suppose that S, € WFsc(w+), 
and let 7 : K ^ ^W be a gener alized broken bicharacteristic in WFsc('u+) with 
7(^0) = C- Then, by Lemma |5.11 , r(— 00) = limt^_oo r^jit)) exists. If t(— cxj) > 



then r(7(t)) > — ^a(A) for large negative t, contradicting that r < — ■\/a(A) 



WFsc(u+). Thus, r(— 00) < 0, and hence by Lemma p.ll| , 7((— cx), to])ni?_(A) 7^ 



so ^ € $+(i?_(A)). We thus conclude that (|8^ holds. D 



A pairing argument immediately shows R{X ± iO)v also exists for distributions 
V e C^°°{X) with wave front set disjoint from the incoming and outgoing radial 



set respectively. Combining it with the propagation theorem, Theorem 2.2, we can 
deduce the following result. 



Theorem (Theorem 2A). Suppose that H is a many-body Hamiltonian, and A ^ 
A. Suppose also that v e C-°°(X) and W¥sc{v) n $-(i?+(A)) = 0. Let uf = 
_R(A + it)v, t > 0. Then uj' has a limit u+ — R{X + iO)v in C^°^{X) as t ^ and 
WFsc(u+)n$_(i?+(A)) =0. Moreover, WFsciu+) C $+(i?_(A))U $+(WFsc(w)). 
The result also holds with R+{X) and i?_(A) interchanged, R{X + it) replaced by 
R{X — it), <&+ and <&_ interchanged. 

Proof. As mentioned above, the first part follows from the self-adjointness of H, so 
that for t > 0, w e C-°°(X), / e C°°{X), we have v{R{X + it)f) = R{X + it)v{f); 
recall that the distributional pairing is the real pairing, not the complex (i.e. L^) 
one. The wave front statement of Theorem ^T| and the assumption on v show the 
existence of the limit u+ = i?(A + zO)ti in C~°°{X) and that in addition WFg"' {u+)C] 
$_(i?+(A)) = for every I < -1/2. The statement WFsc(u+) n $_(i?+(A)) = 
follows from a microlocalized version of the Mourre estimate due to Gerard, Isozaki 
and Skibsted H; see UM or |34] for a detailed argument. The final part of the 



conclusion follows from Theorem 2.2, much as in the previous proof. D 



9. POISSON OPERATORS 

The propagation of singularities theorem, especially in its form for the resolvent. 



Theorem 2.4, has immediate consequences for the wave front relation of all S- 



matrices. To see this, recall the definition of the Poisson operator from (2.29). 
First, let 

(9.1) la = Y. ^b, ia = {Ia/Xa)\c^ & C^{Ca). 

CalCb 

SO la is C°° near C^ with simple vanishing at C^, so la/xa is C°° there. Then, 
for A e (€0,00) \ A and g G C'^{C'^), there is a unique u G C^°^{X) such that 
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{H-X)u 



0, 


and u 


has the form 












u =e~ 


i\/A-e„/3;^dimC„/2^j/„ 


/2v^ 


-Eo 


IJ 


- + i?(A 




u_ 


= ((^'^)*Va)«.9) + 


xv'_, 


/ 


e 


c-(x), 



(9.2) ^ '•'' 

and w^ is conornial to the boundary of [X;Ca], with infinite order vanishing off 
the front face (the inverse image of Ca under the blow-down map), and has a full 
asymptotic expansion at the front face of the form 

oo 2k+2 

(9.3) 5]5]a;'=(logx)^a,(y,Z), 

fc=0 s=0 

Gj smooth in y and Schwartz in Z where we used the notation of ( |3.8|) -( p^ for 
the coordinates in the interior of the front face of [X;Ca]- (Uniqueness holds if 
we merely require that / € Hsi ^ (X), e' > 0, due to Isozaki's results M, [l5| 



applied to the difference of two solutions of the form (9.2), see pq; existence is a 
direct consequence of the following parametrix construction.) The Poisson operator 

Pa.+ W is the map 

(9.4) F„,+ (A) : C^{C',) ^ C-^\X) defined by Pa^+Wg = u. 

The scattering matrix is then given by the following formula, which is related to 
the representation formula of Isozaki and Kitada |l7| : 

(9.5) Sp^iX) - i ((g - A)P^,_(A))*P„,+(A), 

2^^A- ep 

A > max(ea,e^), A ^ A where P^,_(A) is a microlocalized version of P^^_(A), 
microlocalized near the outgoing radial set (see the following paragraphs for details). 

This formula reduces the understanding of the structure of the S-matrices to that 
of the Poisson operators. These in turn can be described in two steps. First, one 
constructs the Poisson operators approximately in the incoming region of phase 
space, and then one uses the propagation results for the resolvent to obtain the 
wave front relation of the Poisson operators. 

A very good parametrix, which we here denote by Pq,,+ (A), for Pa.+{X) in the in- 
coming region of phase space, i.e. where r is close to -y/A — £«, has been constructed 
by Skibsted [^ in the short-range and by Bommier ||l| in the long-range setting 
using a WKB-type construction under the assumption that €a S spec^(iy), i.e. 
Cct is below the thresholds of iJ°. Thus, writing their construction explicitly for 



the class of potentials that we are considering and using the notation of (3.8)-( |3.9|) 
for the coordinates in the interior of the front face of [X;Ca], they construct an 
operator 

(9.6) P„.+ (A):C--(C:)^C— (X) 

by constructing its kernel 



(9.7) K^^^+ = ^-i^/^=7^cosdist(y,y')/x^^I^iy'}/2./^^7^^^^^^ y^ ^^ y>^^ 

a+ G C°°{[X;Ca] X Ca) vanishes to infinite order at the 'main face' (the lift of 
Co = dX to [X; Ca]; this means that a+ is Schwartz in the fibers of the blow-down 
map, which are X"", i.e. it is Schwartz in Z), and a+ is supported where y is near 
y', which in turn implies (by looking at the phase of the exponential) that r is 
near y/X — £„. Note that the phase is a multiple of cosdist(?/, y') — y ■ y' , where we 
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consider Ca as the unit sphere in Xa, y,y' G Ca- In terms of EucUdean variables 
Wa = y/x on X", 

(9.8) Kl^ = e-'^^^^^'--y'\war'-^y'^/^^^^!i+iwa,y',Z), 

with a_|_ poly homogeneous, of degree in Wa, smooth in y', Schwartz in Z. 

The construction of K^ _,_, i.e. finding an oscillatory function of the form ( |9.7D 
satisfying that {H — X)K^ _,_ is in C°°{X x C^) at least near y = y', where H ~ \ 
is applied to i^J!, _|_ in the left factor X, is essentially given by a non-commutative 
WKB-type procedure. Since one already has the correct phase function, one only 
needs to solve transport equations starting at the diagonal y = y' , which, however, 
are operator-valued. To remove the errors in these equations, if" — £„ has to 
be inverted on the orthocomplement of its L^ null space. Since iJ" — €„ has a 
parametrix in Vl/g^. ' {X°',C'^) (because £„ is below the set of thresholds), this is 
possible with the resulti ng (g eneralized) inverse being in ^g^ ' (X'',C"). This gives 



a kernel Kq of the form ( |9.7D , defined near y = y' , x = 0, and satisfying (H — X)Ko 



is in C°°(X X C^) in the region where the transport equations have been solved 
(see [^, |l| for more details). Multiplying it by a cut-off function (jj G C°°{C'^ x C^) 
supported near the diagonal, identically 1 in a smaller neighborhood of the diagonal. 



we obtain a kernel K which is still of the form (9.7), now globally defined, and 
{H - X)K is in C°^{X x C^) near y = y' . Finally, we take ip G C^i,^) identically 1 
near A, and apply i^{H) £ '^g^' i^^C) to K in the left factor X, i.e. we compose 
the operator P given by the kernel K with ip{H), so that Pa,+ {X) ~ 'ip{H)P. 
Writing out the composition explicitly, similarly to how one shows that ^sc(-^,C) 
preserves oscillatory functions, see |3J, Section 4], shows that K^^ + i^ indeed of 



the required form, (9.7). (Recall that this is based on conjugating an Sc-ps.d.o. 
by an oscillatory function, observing that the resulting kernel is that of an Sc- 
ps.d.o, and using that Sc-ps.d.o's map smooth functions on X to smooth functions, 
which is quite straightforward to see by an explicit integral representation of the 
kernel.) In addition, {H - X)Pa^+{X)g = ^AH){H - X)Pg for all g e C~°°{C'^), 
which will allow us to use that {H — X)K is in C°°{X x C^) near y — y' . Also, 
iP-P^,+ W)g = [Id -^{H))Pg=yiH)iH-X)Pg where iPis) = (l-^(s))/(s-A) 
so V G '^irhg(^)' V'(-H^) G *sc i^T^)- III view of the corresponding property of 
{H — X)K, K^ + ~ K vanishes to infinite order near y — y' , x — 0, i.e. K^ ^ and K 
are 'microlocally equal' there. 

We note that in fact the construction of the kernel can be continued up to Causing; 
we just need to cut off a+ shortly beforehand. In particular, if a is a 2-cluster, then 
the construction works until reaching the outgoing radial set. We also remark that 
if the Vb are Schwartz, then a product decomposition immediately gives Pa.+{X), 
without a need to remove errors, so in that case one can allow Eq, to be arbitrary. 

The main properties of Pc,+ (A) are that for any g e C-°°(C^), h e C^iC'J, 
e' > 0, 

(9.9) WFsc(Pa,+ (A)5) C ^^T5, (Ca;X), 

^g ^Q^ P«,+ (A)/i - e-^^^^/^x'''^''^/'x"'/'^^^ii{7rr^c.){{<h) + xv'), 

v as in (|9.3D, 
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(9.11) iH~X)P^,+ {X)heC^iX), 

(9.12) 

WFsc(.Pa,+ (A)g) C{{y,T,fi) : y e supp.g, r = y^X- ea, /^ = 0} 

U {e e S(A) \ (i?+(A) U i?_(A)) : 3C G WF(5), ^ e 7a.-(C)}, 



(9.13) 



WFsciiH - X)P^,+ {X)g) 



C {e e S(A) \ (i?+(A) U R-iX)) : 3C € WF(5), ^ e 7o.-(C)}- 

See [ p7| Appendix A] for a discussion of these mapping properties if a is the free 
cluster; the general case is similar since we are working in the regular part of Ca, 
and a-i- is rapidly decreasing in X". However, since the notion of wave front set is 
more complicated than at the free cluster, we briefly indicate how to prove these 
statements. Before doing this, we discuss the consequences of (|9.9D-( 0^ ) for the 

actual Poisson operator Pq,_+(A). 

The Poisson operator, first defined on C'^{C'^) by ( |9.4[ ), is given by 



(9.14) 



P„,+(A) = Po,+ (A) - i?(A + tO)iH - A)P„,+ (A). 



Note that for g e C^iC'^), (H - A)P„,+ (A)5 G C°°(X) by (|1|), so the resolvent 
can certainly be applied to it. For g G C^{C'^), the right hand side of (9.14) is of 
the form ( |9.2| ) by ( p. 10 ), so it indeed yields the Poisson operator by its definition, 
(P^). Similar results hold for the Poisson operator P^^_(A) with outgoing initial 
data. Theorem |2.4|, (9.13) and (9.14) immediately show the following proposition. 



Proposition. (Proposition 2.6) Suppose that H is a many-body Hamiltonian, X ^ 
A, and (2.34) holds. Suppose also that g £ C'^{C'g). Then 

(9.15) WFsc(Pa.+ (A)g) \ R+{X) c $+(P-(A)). 

In addition, Pa.+ {X) extends by continuity from C'^{C'^) to distributions g € 
C-^{C'^) with (R+iX) X WF(5)) n 7^<,,+ = (i.e. 7^„,+ (WF(5)) n R+{X) = %). If 
g is such a distribution, then 



(9.16) 



WFs,(P„^+(A).g) \ R+{X) C $+(i?_(A)) U7^„,+ (WF(.g)). 



We now return to the proof of ( |9.9D -( |9.13D . First note that (|9.9D follows from 
the fact that supp a+ is disjoint from Causing x _D x X° for any D C C'^ compact. 
Next, ( |9.10 )-( [).ll ) are direct consequences of the stationary phase lemma. The 
logo; factors appear in the lower order terms (i.e. in v') due to the derivatives 

|, Sections 11-12] for a detailed 



falling on a;'-fa(y')/2VA-e„ ^^g stated in (|9j); see 

description. Note that in addition to the stationary phase lemma, ( |9.11 ) uses that 
{II—X)Pa.+{X)h = ■4}{H){H—X)Ph, and the kernel of (7J—A)P is rapidly decreasing 
near y = y' , where the phase is stationary. The stationary phase lemma also shows 
( p.l2| ) with h G C'^{C'^) in place of g (with the second set on the right hand side 
empty) , since all terms in the stationary phase expansion arise by differentiating h 
and evaluating the result at y, so we obtain rapid decay away from supp ft,. 

Concerning ( 9.12 ) (for g G C~°°(C^)) , we are done if we show that for any 
^ that is not in the right hand side of ( 9.12 ), we can write u = Bv + f with 



B e *g<?°'''(X,C), V G C-°°{X), f G C°°{X), C ^ WF'sciB). For the sake of 
simplicity, and since we only need this weaker statement here, we only state the 
proof with {(y, ^/X — ea,0) ■ y G suppg} in the right hand side of (9.12) replaced 
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by {(?/, ^/X — €„, 0) : J/ e C'^}; the improved statement only requires a more careful 
consideration of supports in ( |9.19| ). 

To find such B, w, /, we consider q e C^("'=T*Xa) as in Section |, p e C°°(X) 
supported near Ca, identically 1 in a small neighborhood of Ca, and take 'i/'o G 
C;?°(R) is identically 1 on supp-0 (recall that Pa.+{X) = ■ip{H)P), and let A = 



-00,0 



pQtpo{H) g \l/gj. ' (^,C) as in Section g. Then, as discussed above (with t]:{H) in 



place of A), the kernel -fC' of APa^+{\) has the form (9.7), i.e. 
(9.17) K' = e-'V^=^a-a7:^2.»/a(a')/2VT=7r^,(^^ ^^ ^^ ^/^^ 

a' e C°°([X;Ca] X Ca), and in addition 



(9 18) (2/' d(y/x){-^J^ - Eq y • y7a^)l(o,y,r)) ^ suppg 

^ a' vanishes to infinite order near (y,y') & Ca >i C'^. 
Suppose also that 

(9.19) yeC;^ (2/, VA-ec,0)^supp<j, 

so by ( p. 18 ), the diagonal y = y' is disjoint from supp a'. This means that, with Tx^ 



denoting the Fourier transform in Xa, J^Xa^Pa,+ iX) is a Fourier integral operator 
in the usual sense (on X* x C'^ where X* is the dual of Xa , which we identify with 
Xa for the sake of convenience) , with values in Schwartz functions on X°, and with 
canonical relation in 

(9.20) iT*Xa \ 0) X (r*c; \ 0) 

due to the previous remark. We also take C G '^'•'{C'^) with compactly supported 
kernel and amplitude (full symbol) c, and suppose that 

(y, £,) £ supp q and (y', () G cone supp cC] S*Ca 
(9.21) 



=^ (y, e) ^ 7a,- (y', C) = 7a,- (y', C) U {(y', VA-ea, 0}. 

(N.B. (y,^) G suppq automatically implies tha t (y, ^) is not in the second set on 
the right hand side, i.e. it is not incoming, by ( 9.19| ).) We claim that under these 



conditions APa^+ WC has a smooth rapidly decreasing kernel Kapc S C°° (-^ ^ C^) . 
This can be seen from the standard FIO calculus (with values in Schwartz functions 
on X"^ if stated this way) if we perform Fourier transform in Xa- Namely, as 
remarked above J-'xa^Pa,+ W is a Fourier integral operator, with canonical relation 



satisfying <^M). Thus, J^x„A^a,+ (A)C is also an FIO. But by ( |9l8| ) and (|^, 
(y',C) e WF'(C') implies that for no (C,y') G T*Xa is (C,y',y, C) in the canonical 
relation oi Txa^-Pct.+ W, hence smoothness of the kernel of JFjjf^APo, +(A)C follows. 
Rapid decay at infinity of the kernel of JFx„^^q,+ (A)C being automatic, this shows 
that Kapc G C°"{X x C^). This can also be seen more directly, writing out 
the composition as an oscillatory integral, in the uncompactificd notation, with 
Wa = y/x, 



KAPc{wa,Z,y") = / e-'^^^^-y'e'(y'-y">i\war'-^y'^'^^^^ 
(9.22) J 

a'(wa,y',Z)c{y',C)dy' dC, 

where Wa ■ y' stands for the Xc-scalar product of Wa G Xa and y' G Ca, we see 
that the phase of the y' integral is only stationary if (y, ^a) G la,-{y\ C), but a'c is 
rapidly decreasing there. 
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Now, given any g e C^°°(C'^), and any U C S*C'^ open such that WF(g) C U, 
we can find C G *°(C^) with compactly supported kernel such that WF'(C) C f/, 
WF'(Id-C) n WF(g) = 0. Thus, (Id-C)^ € C^(C^). Choosing q as above, i.e. 



satisfying ( |9.19D and ( |9.2l| ), we see that 

(9.23) Au = ylP„^+(A)Cg + AP^^+{\){ld-C)g e C°°{X). 
As {Id —i/jq{H))u — 0, this shows that 

(9.24) u = ^o(-ff)w = QMH)u + [MH) - QMH))u = Au+ (MH) - A)u, 

so if g = 1 near {y,0, so (y,^ ^ WF^,(B), B = V^o(i?) ~Ae *s;°^°(X,C), then 
(2/1C) ^ WFsc(u). But U, hence suppc, can be chosen arbitrarily small around 
WF(g); correspondingly for any (y, £,) such that {y, ^ ^ 7q _ (y', C) for any (y', C) e 
WF{g), and (y,C) is not incoming, we deduce that (y,^) ^ WFsc(u). The estimate 
for WFsc((^ — A)u) works the same way; now (H — X)u = ipo{H){H — X)Pg, with 
{H — X)P of the same form as P but with smaller (cone) support. This gives the 
improved estimate. 

Note if we only required finitely many of the symbol estimates on, say, the symbol 



whose left quantization is 5, a conclusion analogous to (9.12) would easily follow 
for any operator with kernel given by an oscillatory function as in (9.7), i.e. we 
could allow a+ to be arbitrary, and there would be no need for the factor ip{H) 
in the definition of P„,+(A). (Of course, ( |9.13 ) could not hold then.) Indeed, v = 



(A + 1)*^^ is of the same form as u, so m = (A + l)^'''i;. Choosing q € C^{^'^T*Xa 
C e ^°(C^) as above, so Q € ^~°°'°(Xa), regarded as an operator on X, we 
deduce just as above that QPa,+ {X)C has a smooth rapidly decreasing kernel, so 
( j.23| ) holds with Q in place of A. Let p <E C°°{X) be a cutoff function supported 



near C'a as above. Then u = {ld-pQ){A + ly^v + f, f = pQu e C°"{X), so if 
B = (Id-p(5)(A + I)"*" were in ^^^'°{X,C), we would be able to deduce that 
i ^ WFg^(i3) imphes ^ ^ WFsc(u). The problem is that q is not a symbol in (Ca, C°) 
jointly since it is independent of the latter. However, note that k can be chosen 
arbitrarily large, so it can be arranged that B satisfies an arbitrary, but finite, 
number of the symbol estimates corresponding to ^S(°°' {X,C), while keeping v in 
a fixed weighted Sobolev space. While this is sufficient for any ps.d.o. argument 
based on weighted Sobolev space estimates, etc., to go through, it is not sufficient 
for ensuring that B is in 'I's(°°' {X,C); that was the reason for the somewhat more 
complicate argument of the previous paragraphs. 

10. Scattering matrices 



In view of the pairing formula (9.5) for the S-matrix, which we restate below. 
Proposition 2.6 allows us to describe the wave front relation of the S-matrices rather 
directly. Recall that the S-matrix from an incoming cluster a to an outgoing cluster 
(3 is given by 

(10.1) Sp^{X) = i {{H - A)P^,_(A))*P„,+(A), 

2i^X- ep 

X > max(eQ,e^), A ^ A. This formula defines 5'/3a( A) as a map C^(C^) -^ C~°°(C^). 
Namely if 5 G C^{C'^) then SpaWg is given by the formula 

(10.2) {S0^{X)9,h)= i (P^.+ (A).9,(g-A)Pg._(A)/i), heC^iC',), 

2iy^X-ep 
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where the pairing makes sense since, by (9.13) with Pf}.~ in place of Pa,+ , {H — 
\)Pp_^{X)h E 5(R"). (Note that (•, •) denotes the complex L^ pairing, linear in the 
first variable, not the real distributional pairing; the two are related by taking the 
complex conjugate of the second term.) 

Before stating the full results on the wave front relation of SpaW, we first 
address the question whether SpaW maps C^{C'^) to C°°{Cl) and extends to a 
map C~°°{C'^) -^ C~°°{C'ij). This is closely related to the question whether the wave 
front set of the kernel of S',3„(A) is a subset oiT*{C'f,xC'J\{{OxT*C'JU{T*C'i,xO)), 
as we discussed in Section 0. By a pairing argument we can answer the first question 
affirmatively if for all g S C^{C'^), h S C~°°(Cj), the two terms being paired 



in (10. 2|), that is Pa,+ {X)g and {H — X)Pp^^{X)h, have disjoint wave front sets. 



Indeed, suppose that Ki = WFsc(P„,+ (A)g) and K2 = WFsc((H-A)P/3,_(A)/i) are 
disjoint. Being wave front sets, they are closed, and as tp{H)Pa,+ {X)g = Pa,+ {X)g 
a ip = 1 near X, ip £ C^(R), we see that Ki is compact, and K2 is also compact by 



( P.13| ) applied for Pp-iX). As discussed in Section^, there exists q £ C°°{^'=T*X), 



q TT-invariant, such that q = 1 on Ki, q = on K2, and q satisfies estimates 



(|6.3|), so we can define an operator A £ '^g^' iX,C) as in Lemma ^j which 



is of the form Q'il){H), tp as above, and Q acts on oscillatory functions. Then 
WF's,(7/;(F) -A)r\Ki^%hy construction, so {t}}{H) - A)P„,+ (A)g £ C°°{X). In 
addition, WFsc(v4*) n X2 = 0, so A*{H - X)Pp-{X)h £ C°°{X) as weh. Then 

(10.3) 
(P«,+ (A).g,(7?-A)P;3,-(A)/i) 
= {{4,{H) - A)P«,+ (A).g, {H - X)Pp,-{X)h) + {P„,+ {X)g,A*{H - X)Pp.-{X)h), 

with the equality a priori valid for g, h smooth, for then the second factor is 



Schwartz, shows that the pairing (10.2) makes sense (i.e. extends by continuity 
from smooth data) if Ki and K2 are disjoint. 

This is certainly the case if a and (3 are either the free channel or 2-cluster 



channels satisfying ( 2.34 ) (including the possibility that one is free, the other is 



such a two cluster channel). In fact, if /3 = is the free channel, a is arbitrary 



(but satisfies (2.34)), A > (so with the notation of Theorem 3.1 a(A) = A) then 



on WFsc((i/ - X)Pp,^{X)h), -VX < T < VA by (g, while WFsc(Pa,+ (A)g) C 
<I>+(P_(A)) by Proposition |2.6| , so as t is a decreasing function along generalized 
broken bicharacteristics, r < — •\/a(A) = — vA on it, so these t wo se ts are indeed 



disjoint. On the other hand, if /3 is a 2-cluster channel satisfying ( 2.34 ), a arbitrary, 
then one can arrange (by continuing the kernel construction until close to the — 
outgoing, i.e. + incoming, set) that on WFsc{{H — X)Pp,-{X)h), t is arbitrarily 
close to -^A — ep, which is again disjoint from the set r < — -y/a(A), hence from 
WFsc(Pa,+ (A)5) C $+(P_(A)). This sufiices to prove that Sp^iX) : C°°iC'^) ^ 
C°^{C'f^) in this case. By a similar argument, or by duality, Spa{X) extends to a 
map C~°°{C'^) -^ C~°°{C'ij) if either a is a 2-cluster channel or the free channel. 
In other cases Proposition |2.6| and ( |9.13D do not allow us to conclude that the 
wave front sets of the two terms are disjoint, and correspondingly we cannot expect 
that Spa{X) either preserves smoothness or extends to distributional data. We 
summarize these results in the following corollary. 
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Corollary 10.1. Suppose that H is a many-body Hamiltonian, A ^ A, and (2.34) 
holds. If P is a two-cluster channel or the free cluster then SiSaW preserves smooth- 
ness: S'/5q(A) : C'^{C'^) — > C^iC'f^). If a is a two-cluster channel or the free cluster 
then 5/5q(A) extends to a map on distributions: Spa^X) : C~°^{C'^) -^ C^°°{C^) 

The pairing argument of course allows us to describe the wave front relation of 
the S-matrices in general. If g E C-°°{C'^) satisfies 7^a,+ (WF(ff))^i?+(A) = 0, then 
Pa.+ {X)g is defined, and its wave front set is given by Proposition |2.6| . Similarly 
if /i e C^°°{C'ij), then Pp^-{X)h is always defin ed a nd the wave front set of its 
'error as a generalized eigenfunction' is given by ( p.l3| ). Thus, under the additional 
assumption that these two wave front sets are disjoint, which holds if 

(10.4) {7^,+(C) : C e WF{h)} n (7^„,+ (WF(.g)) U $+(i?-(A))) - 0, 

we see that {Si3a{X)g, h) is defined via ( |10.2[ ). Since h is arbitrary, except that its 
wave front set satisfies ( |10.4D , this allows us to conclude that if i^ C S*C'^ is such 
that 

(10.5) {7/3,+ (C) --CeKjn (7^„,+ (WF(g)) U $+(i?-(A))) = 0, 



then WF(5'/3a(A)5) n i^ = 0. But by the definition of Ufja and Ufj-, ( |10.5D is 
equivalent to 

(10.6) 7^/3a(WF(.g)) n is: = and 7^;^^_(i?_(A)) n K ^ $. 

(See the proof of Theorem 12.4 in pm for a more detailed version of the pairing 



argument.) We have thus completed the proof of Theorem 2.8 



Theorem. (Theorem ^.g) Suppose that H is a many-body Hamiltonian, A ^ A, 



and (2.34) holds. Then Spa{X) extends by continuity from C^ (C^) to distributions 
g e C-°°{C'J with 7^a^+(WF(g)) n i?+(A) = 0. If g is such a distribution, then 

(10.7) WF{Sfi^iX)g) C 7^-l_(i?_(A)) U np^{WF{g)). 

Appendix A. Geometric description of broken bicharacteristics 

The purpose of the appendix is to give a more geometric description of the 
generalized broken bicharacteristics, provided that the set of thresholds, Ai, is 



discrete, that is to prove Theorem 5.10. We also analyze four-body scattering 



where we can also obtain a similar simple geometric description for generalized 
broken bicharacteristics. 

It is convenient to break up the analysis of the structure of the generalized 
broken bicharacteristics into a local and a global step. Although until now we 
used the expression 'broken bicharacteristic' as a synonym for the cumbersome 
term 'generalized broken bicharacteristic', we now give a stronger definition for the 
former which is only valid in this Appendix. 

Definition A.l. Suppose / is compact, 7 : / — > S(A) is a generalized broken 
bicharacteristic. We say that it is a broken bicharacteristic if there exist to < ti < 
. . . < ifc, / = [toytk], such that for each j, 7|[(^t , ^i is the projection of an integral 
curve of ^"^Hg to S(A). In general, if / is not compact, we say that 7 is a broken 
bicharacteristic if 7I j is a broken bicharacteristic for all compact intervals J <Z I. 

Definition A. 2. 
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(i) An A^-body Hamiltonian has a locally broken bicharacteristic relation if for 
all / compact, and for all generalized broken bicharacteristics 7 : / ^ ^(•^); 
there exist to < ti < . . . < tk, I = [to,i/c], such that for each j, j\[t-.t+i] is 
the projection of an integral curve of ^"^Hg to S(A). 

(ii) We say that an A^-body Hamiltonian has a globally broken bicharacteristic 
relation if for all generalized broken bicharacteristics 7 : M ^ '^Wi there 
exist to < ti < . . . < tk, such that 7|(-oo,to]> 7[tfc,+oo)j as well as for each j, 
l\[t-.t+i]: are the projections of integral curves of ^'^Hg to S(A). 

When we talk about the length £(7) of a generalized broken bicharacteristic 7, we 
mean the length of its projection to §"^^ (recall that the projection is Lipschitz). 

Lemma A. 3. Suppose that H is a many-body Hamiltonian. Then there exists I > 
such that for every generalized broken bicharacteristic 7 : [a,P] — > S(A) of length 
^(7) !i I there exists a cluster a such that the image of "f lies in UcaCCb'^'^l'c" -^b- 

Proof. For each point p e S"^^ there exists a unique a such that p E C'^. Also, 
there exists Ip > such that expp{Bi^{0)) hcs in (Jc^cCt^'^T^'^b, where Bi^{0) 
denotes the open ball of radius Ip in TpS"^^ with respect to the standard metric. 
Now, {exp (Bj /2(0)) : p G E>"'~^} is an open cover of the compact set §"~^, so it 
has a finite subcover corresponding to some points, say, Pj, j — 1, . . . ,m. Let I — 
min{Zp^./2 : 1 < j < to}. If 7 is as above, let j be such that 7(a) G expp_{Bi^^^{0)). 
For any t £ [a,f3], the distance of 7(i) and 7(a) is bounded by the length of 7, 
hence by I < lpj/2, so the image of 7 lies in expp . {Bi^ _ (0)) which in turn lies in 

Uc,cc,''=7B,Xh for some a. D 

We next give an upper bound for the length of a broken bicharacteristic if Ai is 
discrete. Note that Ai is bounded, so under this assumption it is finite. 

Lemma A. 4. Suppose that Ai is discrete, and let Ci be the number of its elements. 
Suppose that 7 : / ^ ^W is a broken bicharacteristic, I an interval. Then £(7) < 

Cm. 

Proof. We recall from pa the explicit arclength parametrization, s = s{t), of the 
integral curves of ^^Hg with kinetic energy cr > 0. In terms of this parameterization 
for 7L. j/i, t{s) — ^Joj cos(s — So) where s varies in a subintcrval of (sq, sq + tt), 
and \tj,t'A is such that 7|[f ^t'.j is an integral curve of ^'^Hg with kinetic energy aj. 
Since t o 7 (which we usually just write as r) is monotone decreasing, this shows 
that the total length of the segments of 7 which are integral curves of ^"^Hg with 
any given kinetic energy cr > 0, is at most tt. If cr = 0, the bicharacteristic segment 
is constant. Since a must be such that A — cr G Ai, there are Ci possible values of 
cr, which proves our estimate. D 

Next, we note that in three-body scattering the uniform upper bound tt holds 
without any assumptions on the structure of Ai, 

Lemma A. 5. Suppose that H is a three-body Hamiltonian. Suppose that 7 : / ^ M 
is a broken bicharacteristic, I an interval. Then £(7) < tt. 

Proof. In three-body scattering kinetic energy is constant along generalized broken 
bicharacteristics (essentially because there are no positive energy bound states), so 
the proof of the previous lemma applies and gives the desired conclusion. D 
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Our strategy to analyze generalized broken bicharacteristics is to divide them 
into pieces of length < I, and reduce the analysis to that in a subsystem, by virtue 



of Lemma A. 3 



Lemma A. 6. Suppose H is an N-body hamiltonian all of its proper subsystems 
tiave a globally broken bicharacteristic relation. Then every generalized broken 



bicharacteristic of H of length < I, I as in Lemma A.S, is a broken bicharac- 
teristic. In addition, if in all of the proper subsystems the maximum number of 
breaks in a broken bicharacteristic (defined over R and of arbitrary energy^ is at 
most M^-i, then every generalized broken bicharacteristic of H of length < I has 
at most 2Mn-i + 2 breaks. 

Proof. Let 7 : [a,j3\ -^ S(A) be a generalized broken bicharacteristic of length 
^(7) < h so the image of 7 lies in a region {JcadCb^'^T^'Xh for some a. Thus, 
we can use local coordinates around a for describing 7. By Lemma p.q and the 
argument preceeding it (showing that r]a cannot change sign more than once), 
there exist some points a',/?' such that on [a, a') and on (/?',/?] the image of 7 is 
disjoint from '^'^T^, Xa, and on [a',/?'], 7 is an integral curve of ^'^Hg, where some 
intervals may be empty or reduce to a point. Consider the interval (/3',/9] for the 
sake of definiteness, and assume that it is non-empty. Let S{t) be given by 

(A.l) Sit) = - y (1 + \Za{t')\y/'\Za{t')\-' dt', 

SO S{t) is the solution of the ODE 

(A.2) dS/dt={l + \za{t)\''f'y\za{t)\, te{P',f3], 5(/3) = 0, 

where we wrote |^a(i)| = ka(7(0)l- Thus, S is C^ and its derivative is positive, so 
the same holds for its inverse function, S^^, defined on an interval J. We denote 
by Sh» (A — ICaP) the characteristic variety of the (proper) subsystem Hamiltonian 
H" at energy A - |Cap. Now let 7 : J -> S]ffo(A - l^aP), given by 

(A.3) j{s) = iz,{S-\s))/\z,{S-\s))\,,.,{S-\s))), 

so in terms of Euclidean coordinates, 

(A.4) 7(5) = iw^iS-\s))/\w^{S-\s))\,eiS-\s))), 

It is straightforward to check that 7 is a generalized broken bicharacteristic of iJ"; 
the change of parameters accounts for the change in normalization of the rescaled 
Hamilton vector fields. In fact, ^"^Hg = {w)Hg, while its analog in the subsystem 
is {w°-)Hg^, and the quotient (w°)/(w) is |2:a|/(l + |zaP)^/^. Thus, by the hypoth- 
esis, 7 is a broken bicharacteristic. Since £,a is constant along generalized broken 
bicharacteristics, we conclude that on (/?',/?], j\[/3'^p] is a broken bicharacteristic if 
and only if 7 is, so by the hypothesis it is a broken bicharacteristic, and it has as 
many breaks as 7, so if we assume uniform bounds in the proper subsystems, at 
most Mjv-i- 

A similar estimate holds for [a, a'), so on [a,/3], 7 is a broken bicharacteristic, 
and if we also assume uniform bounds in the proper subsystems, it can have at 
most 2A/Ar_i + 2 breaks, proving the lemma. D 

Corollary A. 7. Suppose H is an N-body hamiltonian, and all of its proper subsys- 
tems have a globally broken bicharacteristic relation. Then H has a locally broken 
bicharacteristic relation. 
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Proof. If 7 : / ^ "^W is a generalized broken bicharacteristic, / compact, then its 
total arclength is finite, so dividing it up into segments of length < I and applying 
the previous lemma proves the proposition. D 

Proposition A. 8. Suppose H is an N -body hamiltonian with Ai discrete, all of 
its proper subsystems have a globally broken bicharacteristic relation, and in all 
of these subsystems the maximum number of breaks in a broken bicharacteristic 
(defined over M and of any energyj is at most Mjsi-i. Then H has a globally 
broken bicharacteristic relation and there is a constant Mj^ such that every broken 
bicharacteristic of H of any energy has at most M^ breaks. 

Proof. Due to the previous corollary, every generalized broken bicharacteristic of H 



is a broken bicharacteristic. By Lemma A. 4, its length £(7) satisfies the inequality 



^(7) ^ CiTT. Dividing it up into pieces of length < I, of which (we can arrange 
that) there are at most Ciir/l + 1, shows that the total number of breaks is at most 
{Ciir/l + 1){2Mn^i + 3) which is independent of 7, proving the proposition. D 



An inductive argument completes the proof of Theorem 5.1C, which we restate 



Theorem A. 9. Suppose that H is a many-body Hamiltonian and Ai is discrete. 
Then H has a globally broken bicharacteristic relation, and there is a constant M 
such that every generalized broken bicharacteristic of H has at most M breaks. 

We next discuss what happens if Ai is not discrete. 

In the following we let Co > be a constant such that for s,s' G (0, tt), 

(A.5) Icoss-coss'l >Co"^|s-s'p. 

Thus, Cq is essentially the Holder constant in the Holder of order i estimate for 
arccos. 

Lemma A. 10. Suppose that 7 : / ^ M is a broken bicharacteristic, I compact, 
[tjjt'A, j = 1,... ,k, are subintervals of I such that the open intervals {tj,t'j) are 
disjoint. Suppose that the kinetic energy over the interval [tj,t'A is aj, and let 

k 

(A.6) Ar = 5^r(7(i;))-r(7(t,)). 

Let ii^j) be the length of ^j = 7|[t ,t']. Then 

k k 

l^-'7 j = l j = l, a-j>a 

< Co(min{a, : a, ^ Q, j = 1, . . . , fc})-i/2 fci/2|Ar|. 

Proof. We again use the explicit arclength parametrization, s ~ s{t), of the integral 
curves of ^"^Hg with ctj > 0. Thus, t{s) = .VctJ cos(s — sq) for 7^ where s varies in 
a subinterval of (sq , sq + tt) . This gives 

(A.8) Coaj'/'kisit'^)) - r(s(i,))| 1/2 > ,(t^.) _ ,(i^.). 

If (Tj — then both T{s{t'A) — T{s{tj)) and s(i') — s{tj) vanish. Summing over 
j and applying the Cauchy-Schwartz inequality to the left hand side proves the 
lemma. D 
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Note that in Lemma |A.(j , we really only need to assume that the estimates A/jv_i 
in the subsystems are uniform over bounded sets of energy to obtain a uniform 
estimate Mm over bounded sets of energy for H . 

Corollary A. 11. Suppose H is an N-body hamiltonian all of its proper subsystems 
have a globally broken bicharacteristic relation, and in all of the proper subsystems 
the maximum number of breaks in a broken bicharacteristic (defined over R) is at 
most Mjv-i- Suppose that 7 : / — > S(A) is a generalized broken bicharacteristic 
of H (hence a broken bicharacteristic under these assumptions), and suppose that 
there exist cq > 0, m > 0, such that the kinetic energy assumes at most m values 
which are less than cq. Then there is a constant M^ that depends only on cq, m, 
and H , such that 7 has at most M^ breaks. 

Proof. We only need to show that there exists M^ as above such that for every 
compact interval J C /, 7|j has at most Afjv breaks. Let n denote the number of 
breaks in 7|j. Note that r o 7 is monotone decreasing and it is bounded, with a 



bound given by (A — infAi)^/^ (since A — r^ — |/^P G Ai in E(A)). By Lemma A. 10 , 
the total length of bicharacteristic segments with kinetic energy at least cq is at 
most 2CoCq ' n^/^(A— inf Ai)^/^, while by the proof of Lemma A. 4, the total length 



of bicharacteristic segments with kinetic energy less than co is at most tott. Thus, 

(A.9) ^(7|j) < 2Coc~^''^n^''\\ - inf Ai)^/^ + mn. 

On the other hand, dividing 7|j into segments of length < /, of which we can 



arrange that there are at most £(7! j)/Z + 1, and applying Lemma A. 6, shows that 

(A.IO) 

n<(2Mw-i + 2)(£(7|j)// + l) 

< (2MAr„i + 2)(1 + m7r/0 + 2{2Mn^i + 2)CoCo ^/^^^/^(A - inf Ai)i/V/- 

Dividing through by 71^ " gives the desired estimate for n, uniform in the energy A 
as long as A stays in a bounded set. D 



Proposition A. 12. Suppose that H is as in Corollary A. 11 , and A ^ Ai. Then 
there exists Mn depending only on H and inf{A ~ E : X ~ E > 0, E € Ai} such 
that every generalized broken bicharacteristic 7 : / — )■ S(A) has at most M^ breaks. 

Proof. Note first that Ai is closed, so for A <^ Ai, ini{\- E : X-E > 0, E e Ai} > 
0. The kinetic energy tr > on any bicharacteristic segment satisfies A — ct G Ai, 



so (T > inf{A - E : X - E > 0, E G Ai} > 0. Applying Corollary A. 11 with 



Co = inf{A — E : X — E > 0, E e Ai}, m — 0, completes the proof. D 



Combining this result with Lemma A. 5 yields the following theorem 



Theorem A. 13. Suppose that H is a A-body Hamiltonian, and A ^ Ai. Then there 
exists M depending only on H and inf{A — E : X — E > 0, E € Ai} such that 
every generalized broken bicharacteristic 7 : / ^ "^W ^^^.s 'J^ most M breaks. 

While in 4-body scattering one is mostly interested in A ^ Ai, making this result 
useful, in the inductive step, to analyze 5-body scattering, the 4-body thresholds 
must be understood as well, as indicated by the proof of Lemma |A.6| , in particular 
by 7 : J ^ TiH^iX— |CoP) (the notation is the same as in the Lemma), for A— |fop 
can be a threshold of the subsystem H"^. 
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Appendix B. Lagrangian structure of the forward broken 

bicharacteristic relation 

In this section we show that if the set of thresholds of H is discrete, then the 
forward bicharacteristic relation is the finite union of Lagrangian submanifolds, 
and it maps the outgoing radial set to a finite union of Lagrangian submanifolds 
as well. Together with Theorem E^, (a simple modification of) this immediately 



implies the corresponding result, stated in Theorem 2.10, about the wave front 
relation of the various scattering matrices with incoming and outgoing channels 
satisfying the assumptions ( p. 34 ). 



First, recall that if Ai is discrete, every generalized broken bicharacteristic 7 
consists of a finite number of bicharacteristic segments (with a uniform bound on 
the number of these segments) . Let aj be the cluster and aj the channel in which 
the jth segment propagates. Also, let Cj be the cluster at which 7 breaks between 
propagation along the clusters Uj and flj+i, so Cc^ C Ca^ H Ca^^i and the 'break 
point' is in ^"^T^, Ac , i.e. over the regular part of Cc ■ Thus, we can associate a 
string 



(B.l) ai,ai,ci,a2,a2,C2,. ■ ■ , c™, a^+i, a^+i 



to 7, where each term has the same meaning as above. Since there is a uniform 
bound on the finite number of breaks that 7 has, there are only finitely many such 
permissible strings. Moreover, since we are assuming that there is a break at Cj, 
j = 1, . . . ,TO, we can impose that either aj 7^ Cj or aj+i ^ Cj. (Note that the 
requirement that the 'break points' are in '^'^T^, X^. is only necessary to associate 

unique Cj's to 7; it plays no role in the following arguments.) 

It will be convenient to use the linear structure of the collision planes in this dis- 
cussion. Thus, we will think of the relations, just as of the wave front sets, as conic 
sets, invariant under a natural M^ action in the base variables (i.e. in configuration 
space). Correspondingly, we talk about conic Lagrangian rather than Legendre 
submanifolds. Since multiplying a vector field by a non-zero function simply repa- 
rameterizes the integral curves, the bicharacteristics of ^"^Hg and Hg^ agree after 
reparamaterization, so the former are given by the (R"*"-quotient induced) projec- 
tion of straight lines in T*Xa = Xa x A* whose projection to A* is constant, say 
^, such that 1^1^ = 9a{S,) = A — £„ where ea is the energy of the bound state the 
particle is propagating in. We sometimes write |^| for |^|q to simplify the notation. 

We show that the forward broken bicharacteristic relation is given by a finite 
union of Lagrangians, one for each string as in (|B.1| ). To sec this, we write the part 
of the forward broken bicharacteristic relation corresponding to such a string as the 
composite of elementary relations, each of which is Lagrangian, and which intersect 
transversally in the usual sense. Thus, we let A = Ac ,q, ,a ,c _i to be the relation on 
T*Xc X r*Ac _i corresponding to a bicharacteristic emanating from the collision 
plane Ac _i, propagating along aj in channel aj, and then hitting the collision 
plane Ac . That is, A is given by the projection of the end points of bicharacteristic 
segments of Hg^ whose endpoints lie over Ac _i and Ac respectively, to T*Ac x 
T*Acj_i via the projections Traj,Cj_i and TTaj,Cj- We identify A* with Xa via the 
metric ga (induced by g); we do so in particular when talking about T*Xa = 
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Xa X X*. Thus, {w,^,w',^') £ A means that 
(B.2) e = 7r,^,c,(a, r = 7r,^,c,-,(l), ^ = 



w - 



In particular, this shows that A is smoothly parameterized by w G Xc-, w' e ATc _i, 
w 7^ w', so A is a graph over the base Xc x Xc _i of the cotangent bundle (away 
from w = w'). Note that if w = w' , the segment of 7 is trivial, hence 7 is described 
by a different string, namely one from which aj is missing. Thus, the tangent space 
of A at some point p g A, given by (w, w') as above, is spanned by 

V ■ dw + BjV ■ d^ + C-jV ■ d^r, V € X^^ , 
^^■^^ v' ■ d^, + B'^v' ■ d^> + Cy -d^, v' e Xc^_„ 

where Bj G End{Xcj), Bj e End(Xc^_J are given by 



B, = /^ ';^ i\w - W'\' Id^iw - TT.^.cK)) ®iw- ^a„c,{w'))), 



(B.4) ' ' 



Here, as in what follows, the center dot •, e.g. when writing v ■ d^i , etc., simply 
denotes the 9^', etc., component of the tangent vector. Note that A is a Lagrangian 
submanifold of T*Xcj x T*Xcj_i with respect to the usual twisted symplectic form 
Uc- — i^c _i, where tUc denotes the canonical symplectic form on T*Xc as well as its 
lift to the product (by an abuse of notation). As we discuss later, this implies C' = 



C* as well as Bj, B', self-adjoint. In addition, from (B^), \w~T:aj,cj{w')\ < \w — w'\ 
shows that Bj > 0, and indeed it is positive definite if \w — naj,cj{w')\ ^ \w — w'\, 
i.e. ii w' ^ Xc ■ If w' G Xc , then so is ^ = ^, and we see that the null-space of Bj 
is span{^}, its range the orthocomplement, span{^}^, of span{^} in Xc ■ Similarly, 
Bj < 0, and it is negative definite ii w ^ ^cj^i ■ 

With this notation, every element of the part of the forward generalized bro- 



ken bicharacteristic relation corresponding to the string (BJ) is in the composite 
relation 

yrj.Oj .^^a„i-|-i .Qtti + i ,a„i^i,CTTi ^ ^^Cjn^ajn^am,Cjn — i O . . . O J^C2,ct2,a2 ,ci ^ -'^ci ,ai ,ai ,ai 5 

and conversely, every element of the composite relation certainly corresponds to an 
element of the forward generalized broken bicharacteristic relation. For the sake of 



convenience we set cq = ai, Cm+i = cim+i- We now prove that (B.5) is a smooth 



Lagrangian. As indicated in the first paragraphs, a simple modification yields 



Theorem 2.10 if we also prove that the composite relation maps R- (A) n ^'^T^, Xa-^ 



to a finite union of smooth Lagrangians. We do not state this mapping property 



explicitly to avoid an overburdened notation, but we indicate in Remark B.2 how 
it is proved by a simple modification of the following argument. 



Proposition B.l. For each string as in (BT) the composite relation ( B.5 ) is a 



smooth Lagrangian submanifold of T*Xa„^^^ x T*Xai. 

The rest of this appendix is denoted to the proof of this proposition and to the 



indication of the minor changes required for the proof of Theorem 2.10 
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Proof. The proof proceeds by induction. Thus, let 

(tj.Dj Aj — Acj,aj.aj,Cj-i O . . . O Ac^ ,3,2.02, ci ° ■''■ci ,qi ,ai ,ai ! 

we need to show that for each j < m + 1, given that A'_]^ is smooth Lagrangian, 
so is a' = Acj^aj,aj.cj^i ° ^i-i' ^^ ^^"^ usual argument, this certainly follows if 
the intersection of Ac^.Q^.^ajCj-i ^ ^7-1 '^ith the partial diagonal diag in T* Xc^ 



X 



{T*Xc^_^ X T*Xcj_^) X T*Xca is transversal, i.e. for p e (A x t\.]_i) n diag, 
(B.7) rp(A X A«_i) + Tpdiag = Tj,{T*X,^ x r*X,^_, x r*X,^_, x T*X,,). 

Here 'partial' means that we take the diagonal in the central factor T*Xc_i x 
T*Xc_-^, and for the sake of simpler notation we wrote 

(B.8) A = A,^,„^.,^,e,_,. 

In this case the tangent space of the composite relation. A' is given by the projection 
of 

(B.9) rp(AxA«_i)nrpdiag 

to the first and last factors, i.e. to T*Xc- x T*Xcq. The transversality is equivalent 
to 

(B.IO) r(A X T*X,„) + T{T*Xc^ x A«_i) = r(T*X,^ x T*X,^_, x T*X,„) 
over (A x T*Xco) n {T*Xc^ x A5_^), and then TA^ is given by the projection of 

(B.ll) r(A X T*X,,) n T{T*X,^ X A«_i) 

to the first and third factors. Transversality certainly follows if we can find a 
Lagrangian subspace V of TT*Xc _i such that V is in the range of the differential 
of the projection A:_]^ — > T*Xcj_;^, and 

(B.12) TA + TT*Xc^ xV' = TT*Xc^ x TT*Xc^_^ , 

and in this case the Lagrangian subspace 

(B.13) V ^TAoV' d TT*Xc^ , 

given by the projection of TA n {TT*Xc^ x V) to TT*Xc^, is in the range of the 
differential of the projection A' -^ T* X^^ . Since the previous statements involving 
V' referred to tangent vectors, they are essentially equivalent to mapping properties 
of A on Lagrangian submanifolds A' of T*Xc_i ■ Thus, for the sake of convenience 
in notation, we will consider these mapping properties, i.e. we will take V' = TpA'; 
it is simple to reinterpret results in the desired form. 
To be concrete, for j = 1 we take 

(B.14) V' = VI ^ spanjw • d^> + BjV • 9^ : w G AT^J. 

Note that V is certainly in the range of the differential of the projection Aq = Aq ^ 
T*Xci by ( |B.3| ) (the projection is to the unprimed variables in the notation of this 
equation!). If 02 ^ ci, an equally good choice is V = V( = X^ C T*Xci, which 
simply corresponds to the Lagrangian given by plane waves Ci = const. 

Before continuing, we make some general remarks about Lagrangian submani- 
folds A of T*Xa. Suppose that A is (locally) a smooth graph over Xa', i.e. that 
for every p e A, the bundle projection T*Xa -^ Xa has a surjective differential 
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at p (which is thus an isomorphism). Then A is (locaUy) the image of a smooth 
bundle map F : Xa -^ T*Xa, i.e. the composite of F with the projection to the 
base is the identity map. Since T*Xa = Xa x X* naturally, F has the form 
F{w) = {w,F2{w)) where F2 : Xa -^ X*. Moreover, TpT*Xa can be naturally 
identified with the vector space T*Xa itself; so for each w G Xa, (i^2)*|ui = (-^"2)* 
is a linear map from Xa to X*. Let A — A^ be the induced endomorphism of 
Xa via the metric identification of Xa and X* . That A is Lagrangian means that 
for any two tangent vectors V, V € TpA, lu{V, V') = 0, where cj = cja is the stan- 
dard symplectic form given by d£_ A dw, so writing V, V £ TpT*Xa — Xa x X* 
asV = iv,v*), V = {v',{v')*), uj{V,V') = v*{v') - {v')*{v), so in our case, with 
V = F^v, V = F^v' , the Lagrangian condition becomes {F2)*v{v') — {F2)*v'{v) — 0. 
Since the metric identification means that {F2)a.v{v') = ga{Av,v') = Av ■ v' , the 
Lagrangian condition amounts to the statement that A is self-adjoint. 

A similar discussion also applies to T*Xc, x T*Xc-_i with the twisted symplectic 
form Wcj — Wcj_i , and shows that for Ac^^aj,aj,cj-i, which is a graph by ( |B.2| ), with 



the notation of (B^), Bj a nd B '^ are self-adjoint on Xc^ and ^cj_i respectively 



(this of course follows from (B^) as well), while C' = C* . 

Now we return to transversality. In fact, we show the following stronger state- 
ment. Let j be arbitrary, a = aj, d = Cj^i, c — Cj, a — aj, and suppose that 
A' is a Lagrangian submanifold of T*Xd which is (locally) a smooth graph over 
Xd of the form {w',S,') = {w' ,F'{w')) where for each w' € Xd, -F^'ltu' is given by a 
positive operator A' — A'^, on Xd (i.e. A' > 0; A' is automatically self-adjoint as 
discussed above), and suppose in addition that A' is positive definite if d = a. Then 
we show that the relation A — A^caad intersects T*Xc x A' transversally, and maps 
A' g T*Xd to a Lagrangian A C T*Xc with similar properties, i.e. A is (locally) a 
graph of the form {w, F{w)), A — A^, — -F*|«, is positive for every w G Xc, and A 
is positive definite if c ^ a. Note that the transversal intersection property auto- 
matically shows that A is smooth and Lagrangian. Once this is shown, the desired 
result follows by induction starting from j = 1. Indeed, to keep the induction going 
we need to check that A is positive definite if Cj — aj+i, but since this implies 
Cj ^ Qj, i.e. c ^ a, as discussed before, we see that our claim, that A is positive 
definite if c 7^ a, ensures that this holds. 

We remark that for j — I, the Lagrangian subspace Vi of Tp(^^'-fT*Xci defined 



in (B.14) certainly has the desired positivity; the operator A' = A'l is given by 
-Bi, hence it is positive definite unless ci = ai. If instead we take Vi — X^ 
(corresponding to plane waves), then A' = A'l = 0, so it is still positive (though 
not positive definite), which suffices if 02 7^ ci. 

We proceed to show the claimed mapping property of A. We need to show that 
for p = {w, C, w', O e {T*Xc X A') n A 

(B.15) TpiT*X, X A') + TpA = Tp{T*X, x T*Xd). 

But, due to the first summand, d^, d^ are in the left hand side, hence so is v' ■ dw' + 
B'jv' ■ d^', v' e Xd, due to (|B.3| ). Since the vectors wq • dw' + A'wo • d^' , wq G ^d, 



B] A' 



are also in the left hand side due to the first summand, (B.15) will follow if the 

, , , • II' 

block matrix 



is invertible, i.e. if i?' — A' is invertible. But B' < 0, in 
fact negative definite if a 7^ d, while A' > 0, in fact positive definite if a = d, we 



conclude that A' — B', is positive definite, hence invertible. This proves (B.15) 
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Next, we need to find 

(B.16) Tp{T*Xc X A') n TpA. 

Dimension counting shows that we will have found the whole intersection if for 
every v e X^ we find an element of the intersection which is of the form 

(B.17) v-dio+Avd^+vi-dio'+v2- 5^ . 

To do so, fix ?j e Xc, and consider elements of TpA which are of the form 

(B.18) v-d^o + BjV-d^ + Cjvd^' +v' -d^'+B'jv' -d^'+C'jv' -d^, v' e X^. 



For these to be in the intersection (B.16), we need that there exist fo G Xd such 
that the projection of ( B.18| ) to T(^t^i^^r-^T*Xd is equal to vq ■ d^' + A'vq ■ d^', i.e. 
such that 



(B.19) 



^'o, 



CjV + B'jv' = A'vq. 



The existence of such i/ and vq thus again follows from the invertibility of 



L^;- -A' 



which holds since A' — B', is positive definite, and we see that in particular v' = 
(A' — B'A^^CjV. The projection of the corresponding tangent vector in ( B.16| ) to 

T{vj,i)T*Xc is given by 

(B.20) vd^, + {Bjv + C'jv') ■d^=v-dy, + {Bj + C* {A' - B'j)-^Cj)v ■ d^ 

where we used C' — C* . Thus, the projection of the intersection {T*Xc x A') n A 
to T*Xc is a smooth Lagrangian A with tangent vectors of the form 

(B.21) vd^,+Avd^, A^Bj + C*{A' -B'j)-^Cj. 

Since A' — B', is positive definite, the same holds for its inverse, so C*{A' — 
B'A^^Cj > 0. Moreover, Bj > 0, so we conclude that ^ > 0. If in addition 
c ^ a, then Bj is positive definite, hence the same holds for A. This provides the 
necessary inductive step and proves Proposition B.l. D 



Remark B.2. The outgoing radial set, i?_(A) n^^'^T^, X^ gives another example of 

a Lagrangian A' such that the operator A' , defined above, is positive. Namely, a 
simple calculation shows that on fj, = 0, t 



--\/A — (7, cr e Ai, A' is given by 



(B.22) 



A' 



VT 



-(|wpld-i«® w). 



hence A' > 0. Thus, if 02 ^ ci, which holds in this case (otherwise the bicharacter- 
istic were constant!), the proof of Proposition B.l actually shows that the forward 



broken bicharacteristic relation maps it to a finite union of smooth Lagrangian 
submanifolds. 



We now indicate the changes required to prove the S-matrix result. Theorem ^.10. 



Proof. (Proof of Theorem 2.10.) With the notation of (B.l), if ttt, = 0, then 7 is 
an unbroken bicharacteristic, so the endpoints of 7 are related by the antipodal 
relation, which is certainly Lagrangian. So we may assume m > 1. In addition, due 



to the assumptions (2.34), the ^'^Hg bicharacteristics will only break upon hitting 
Cci with ci ^ a (i.e. a-tangential bicharacteristics cannot become normal to a 
without hitting Ca, sing)- Similarly, Cm 7^ b. 
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The only change from the proof of Proposition B.l is that we need to replace the 



first and last Lagrangians in the composition ( |B.5| ) by the Lagrangians Kp C r*C^ x 
T*Xc^ and A„ C T*Xc,xT*C'^ given by the (twisted) graphs oid{-^X- epuj-w') 
and d{—\/X — Eq ^' ■ w) respectively; here we wrote (w, C, w' , ^') and {w, ^, uj' , C') for 
the coordinates on T*C'ij x T*Xc„^ and T*Xc^ x T*C'^ respectively. Thus, on Aq, 
^ = — VA — ec tTq^ci ('^')j so the differential of the projection of Aq, to T*Xci is 
surjective since a ^ ci. In particular, we can choose a Lagrangian subspace T^' of 
T(^yj^^)T* Xc-^ which is in the range of the differential of this projection, and which 
is a graph over X^ of the form [w, Aw) with A positive definite (e.g. the identity). 



Then the proof of Proposition B.l shows that the composite relation 

(B.23) Ac„_Q^^a„,,c„_i O . . • O Ac2,a2,a2.ci ° A^ 

is a smooth Lagrangian. Since the differential of the projection of A^ to T*Xc^ is 



also surjective as 6 ^ Cm, the composition of A^ with (B.23) is transversal, so the 
result is a smooth Lagrangian submanifold of T* C^ x T* C^ . In view of Theorem 2 



this and Remark B.2 complete the proof. D 



Remark B.3. The only reason for the individual composite Lagrangians discussed 
in the above proof not to be cl osed in T*C'i^ x T*C^ \ is that in the definition 
of the elementary Lagrangians, (|B.2| ), we excluded w = w' . Instead, the boundary 
points corresponding to w = w' were included in a composite Lagrangian given by 



a different (shorter) string. As discussed after (B.l), we do not need to assume 



that the break points are over the regular part of the Cc ■ Each of the composite 



Lagrangians corresponding to a string (B.l) is smooth through the singular part of 



the various Cc ', they simply intersect each other there. 
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